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Supervisor’s Foreword 


A promising accelerator light source mechanism called steady-state microbunching 
(SSMB) has been actively studied in recent years. The idea of SSMB is to scale 
the longitudinal focusing of the electron beam in a storage ring from the conven- 
tional radiofrequency range to optical laser wavelengths. The combination of 
microbunching-enabled coherent radiation and the high repetition rate of electron 
beam circulating in a storage ring makes SSMB a high-average-power high-flux 
narrowband photon source, with wavelength extendable to soft X-ray. Such a high- 
flux narrowband light source allows sub-meV energy resolution in angle-resolved 
photoemission spectroscopy, which could provide new opportunities for fundamental 
physics research, for example, to probe the energy gap distribution and electronic 
states of superconducting materials. An SSMB-based kW-level extreme ultraviolet 
(EUV) source is also appealing for the semiconductor industry to be used in EUV 
lithography for high-volume chip manufacturing. 

From the accelerator physics perspective, the six orders of magnitude shortening 
of bunch length compared to that in a conventional storage ring provide challenges as 
well as tremendous opportunities for the accelerator physics research of SSMB. The 
Ph.D. thesis of Dr. Xiujie Deng is devoted to the theoretical and experimental studies 
of SSMB, with important results achieved. The contribution of his thesis can be 
summarized into three categories: first, answer the question of how to realize SSMB; 
second, reveal what radiation characteristics can be obtained from the formed SSMB; 
and third, experimentally demonstrate the working mechanism of SSMB in a real 
machine for the first time. All these achievements, in particular, the first proof-of- 
principle experiment, are of crucial and fundamental importance for the development 
of SSMB. In addition, I believe the efforts of Dr. Deng on precision microbunching 
dynamics will be of growing value to the general accelerator community, as the 
requirement for beam manipulation becomes more and more demanding. I whole- 
heartedly recommend this thesis to all accelerator scientists, especially those who 
work on advanced light sources, and also to the potential users of SSMB. 


Beijing, China Wenhui Huang 
May 2023 
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Chapter 1 A) 
Introduction E 


Particle accelerators as photon sources are advanced tools in investigating the struc- 
ture and dynamical properties of matter, and have enabled advances in science and 
technology for more than half a century. The present workhorses of these sources are 
storage ring-based synchrotron radiation facilities [1-3] and linear accelerator-based 
free-electron lasers (FELs) [4-7]. These two kinds of sources deliver light with high 
repetition rate and high peak brilliance and power, respectively. Some applications, 
however, do need high average power and high photon flux. Kilowatt extreme ultra- 
violet (EUV) light sources, for example, are urgently needed by the semiconductor 
industry for EUV lithography [8]. Another example is that to realize high energy res- 
olution in synchrotron-based angle-resolved photoemission spectroscopy (ARPES), 
which is highly desired by fundamental condensed matter physics research, we need 
the initial radiation photon flux before monochromator is high enough. To obtain 
high average power and high photon flux, a high peak power or a high repetition rate 
alone is not sufficient. We need both of them simultaneously. 

The key of the high peak power of FELs lies in microbunching, which means 
the electrons are bunched or sub-bunched to a longitudinal dimension smaller than 
the radiation wavelength so that the electrons radiate in phase and thus cohere [9— 
11]. The power of coherent radiation is proportional to the number of the radiating 
electrons squared, therefore can be orders of magnitude stronger than the equivalent 
incoherent radiation in which the power dependence on the electron number is linear. 
The Self-Amplified Spontaneous Emission (SASE) scheme [6, 7] of microbunching 
making the high-gain FELs so powerful, however, is actually a collective beam 
instability which degrades the electron beam parameters and the microbunching 
can only be exploited once. The repetition rate of the radiation is thus limited by 
the repetition rate of the driving source, i.e., the linear accelerator. There are now 
active efforts devoted to improve the repetition rate of FEL radiation, for example by 
implementing the superconducting technology. However, the realization of a high- 
average-power, continuous-wave (CW), narrowband, short-wavelength light source 
remains a challenge. 
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Fig. 1.1 A schematic layout of a conventional storage ring (left) and an SSMB storage ring (right) 


A mechanism called steady-state microbunching (SSMB) has been proposed [12, 
13] to resolve this issue. The idea of SSMB is that by a phase-space manipulation 
of an electron beam, microbunching forms and stays in a steady state each time 
going through a radiator in a storage ring. The steady state here means a balance of 
excitation and damping, a true equilibrium in the context of electron storage ring beam 
dynamics. The schematic layout of an SSMB storage ring and its operating principle 
in comparison to a conventional storage ring is shown in Fig. 1.1. SSMB replaces 
the conventional bunching system in a storage ring, namely the radiofrequency (RF) 
cavity, with a laser modulation system. As the wavelength of laser (~ um) is typically 
six orders of magnitude smaller than that of an RF wave (~ m), a much shorter bunch, 
i.e., microbunch, can thus be anticipated by invoking this replacement together with 
a dedicated storage ring magnetic lattice. 

The microbunching in SSMB is from the active longitudinal focusing provided 
by the laser modulator, just similar to the conventional RF bunching through phase 
stability principle [14, 15]. The radiation in SSMB, unlike that in an FEL, is a 
passive process and the radiator can be rather short, for example it can be a simple 
dipole magnet or a short undulator. The SSMB modulator is also much shorter than 
the radiator undulator in a high-gain FEL. Therefore, there is no FEL mechanism 
invoked in the bunching or radiation process in SSMB. If there is some unavoidable 
FEL effects, it needs to be controlled within a safe region to not destroy the steady 
state micobunches. 
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To provide adequate and stable longitudinal focusing such that microbunches can 
be formed and sustained, SSMB requires a powerful phase-locked laser to interact 
with electrons on a turn-by-turn basis. The realization of such a laser system usually 
demands an optical enhancement cavity. A laser cannot effectively interact with the 
co-propagating electrons if the electrons go through a straight line, as the electric 
field of a laser is perpendicular to the laser propagation direction. A modulator which 
bends the electron trajectory transversely is thus needed. The modulator is usually an 
undulator, which is a periodic structure of dipole magnets with oscillating polarity. 
Note that to avoid the head-on collisions, i.e., the Compton back-scattering, between 
the reflected laser and the electrons, a four-mirror optical cavity, instead of a two- 
mirror one, is chosen for the illustration in Fig. 1.1. 

Note that we have not presented explicitly the energy replenish system for SSMB 
in the illustration. The modulation laser in principle can be used to compensate the 
radiation energy loss of the electrons, just like the traditional RF, but this may not 
be a cost-effective choice. Besides, the electron beam current and output radiation 
power will also be limited by the incident laser power. Instead, one may just use 
a traditional RF cavity for the energy compensation. If a larger filling factor of 
the electron beam is desired, the energy supply system could also be one or several 
induction acceleration cavities. In the present envisioned high-average-power SSMB 
photon source, induction linac is tentatively used as the energy compensation system 
and the filling factor of the electron beam in the storage ring can be rather large, for 
example larger than 50%. 

Once realized, SSMB can combine the strong coherent radiation from microbunch- 
ing and the high repetition rate of beam circulating in a storage ring to provide high- 
average-power, high-repetition (MHz to CW) narrowband radiation, with the wave- 
length ranging from THz to soft X-ray. Such a novel photon source could provide 
unprecedented opportunities for accelerator photon science and technological appli- 
cations. For example, SSMB is promising for generating kW-level EUV radiation for 
EUV lithography [16]. Energy-tunable high-flux narrowband EUV photons are also 
highly desirable in condensed matter physics study, such as used in high-resolution 
ARPES to probe the energy gap distribution and electronic states of superconduct- 
ing materials. Ultrahigh-power deep ultraviolet and infrared sources are potential 
research tools in atomic and molecular physics. Moreover, new nonlinear phenom- 
ena and dynamical properties of materials can be driven and studied by high-peak and 
average-power THz sources. Besides high power, SSMB can also produce ultrashort 
(sub-femtosecond to attosecond) photon pulse trains with definite phase relations, 
which could be useful in attosecond physics investigations. 

This dissertation is devoted to the theoretical and experimental studies of SSMB, 
with important results achieved. The contribution of this dissertation can be sum- 
marized as: first, answer the question of how to realize SSMB; second, reveal what 
radiation characteristics can we obtain from the formed SSMB; and third, experi- 
mentally demonstrate the working mechanism of SSMB in a real machine for the first 
time. More specifically, in Chaps. 2 and 3, we have conducted in-depth theoretical 
and experimental studies on single-particle effects vital for the formation and trans- 
portation of microbunching in a storage ring. Chapter 2 is on longitudinal dynamics, 
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while Chap. 3 is devoted to transverse-longitudinal coupling dynamics. Chapter 4 
is the theoretical and numerical investigation on the average and statistical charac- 
teristics of the radiation generated from the formed microbunching. In Chap. 5, we 
report our work on the first successful demonstration of the mechanism of SSMB, 
performed at the Metrology Light Source in Berlin. Finally, in Chap. 6 we present a 
short summary of the dissertation, together with some useful formulas and example 
parameters of SSMB storage rings aimed for kW-level infrared, EUV and soft X-ray 
radiation, respectively. Summarizing, the highlights of this dissertation are: 


e Presents the first proof-of-principle experiment of a promising accelerator light 
source mechanism. 

e Covers precision longitudinal and transverse-longitudinal coupling dynamics in a 
storage ring. 

e Provides useful formulas and example parameters for high-power infrared, EUV 
and soft X-ray light source design. 


The work presented in this dissertation is of fundamental importance for the devel- 
opment of an SSMB-based high-power photon source. 
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Chapter 2 A) 
SSMB Longitudinal Dynamics ciecie; 


In this chapter, we study the single-particle longitudinal dynamics of SSMB. The 
motivation is to answer the question: how to realize the short bunch length and 
small longitudinal emittance in an electron storage ring, as required by SSMB? 
Note that the curvilinear (Frenet-Serret) coordinate system and the state vector 
X = (x, x’, y, y’, z, 8)", with T representing the transpose, are used throughout this 
dissertation. For the longitudinal dynamics without coupling from the transverse 
dimension, what we can play are the momentum compaction and RF systems, for 
SSMB the laser modulators. The momentum compaction is a measure of particle 
energy dependence of the recirculation path length 


_ AC/Cy 1 f Duis) 9, B 


a= = 
AE/Eo Coj) p(s) 


where Co is the ring circumference, Ep is the particle energy, D, is the horizontal 
dispersion which is a measure of the energy dependence of particle horizontal posi- 
tion, p is the bending radius. Considering the energy-dependent velocity, the particle 
energy dependence of the revolution time can be quantified by a parameter named 


phase slippage factor 
AT/To 1 

= =a i (2.2) 
AE/Eo y? 


n 


with y the Lorentz factor. For linear dynamics, the phase slippage to longitudinal 
dimension is like the drift space to transverse dimension, while the RF kick in lin- 
ear approximation to longitudinal dimension is like the quadrupole to transverse 
dimension. The difference is that the sign of phase slippage can either be positive or 
negative, while the drift space can only have a positive physical length. To account 
for the impact of local or partial phase slippage on the evolution of longitudinal optics 
around the ring, Courant-Snyder analysis can be invoked for linear dynamics study 
beyond adiabatic approximation. Such new derivations are necessary to accurately 
describe the dynamics of the SSMB mechanism. Usually there is only one RF cavity 
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in a storage ring, the longitudinal optics can be manipulated with more freedom with 
multiple RFs. For example, the strong focusing principle can be implemented in the 
longitudinal dimension to realize ultrashort bunch length, not unlike its transverse 
counterpart. For nonlinear dynamics, both the nonlinearity of the phase slippage and 
the sinusoidal modulation waveform can lead to subtle and rich beam dynamics. In 
the following we will investigate along this brief review. Parts of the work presented 
in this chapter have been published in Refs. [1—4]. 


2.1 Linear Longitudinal Dynamics 


2.1.1 Longitudinal Courant-Snyder Formalism 


SSMB means an ultrashort electron bunch in an equilibrium state. One successful 
method of realizing short bunches in an electron storage ring is the implementation of 
a quasi-isochronous lattice, which means particles with different energies complete 
one revolution using almost the same time. The reason behind is the well-known 
oz X 4/|n| scaling law of the “zero-current” bunch length given by Sands [5], in which 
n =q — = is the global phase slippage factor of the ring as introduced just now. 
However, from single-particle dynamics perspective, there is a fundamental effect 
limiting the lowest bunch length realizable in an electron storage ring originating 
from the stochasticity of photon emission time or location. This stochasticity results 
in a diffusion of the electron longitudinal coordinate z even if the global phase 
slippage of the ring is zero as we cannot make all the local or partial phase slippages 
zero simultaneously. The partial phase slippage factor from sı to s is defined as 


nls “oe -{ (=e =) as (2.3) 
ee Csi pl) P 


The physical picture of the partial phase slippage and quantum excitation in both the 
particle energy and longitudinal coordinate, therefore the longitudinal emittance, is 
shown in Fig. 2.1. 

Due to this quantum diffusion, there exists a lower bunch length limit and the 
energy spread diverges when the bunch length is pushed close to the limit. This 
effect is of vital importance for SSMB and other ideas invoking ultrashort electron 
bunches or ultrasmall longitudinal emittance in storage rings. It is first theoretically 
investigated by Shoji et al. [6, 7], and recently more accurately analyzed by us using 
the longitudinal Courant-Snyder formalism [2, 3, 8]. The key to understanding the 
effect is to change from the global viewpoint to a local one, i.e., the quantum excitation 
at different places around the ring actually contribute to the longitudinal emittance 
with different strengths, just like its transverse counterpart. 

For an accurate analysis of this effect, here we invoke Chao’s solution by linear 
matrices (SLIM) formalism [9]. SLIM is an early effort to generalize the classical 
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Fig. 2.1 Physical picture of Bliceon emission 
the partial phase slippage a 
factors and quantum H(S2, So) 
excitation. Particles undergo 
diffusion in both the particle 
energy and the longitudinal 
coordinate in each turn, 
giving rise to longitudinal 
emittance growth 
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(S150) 
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Courant-Snyder theory [10] from 1D (2D phase space) to higher dimensions. It 
invokes 6 x 6 transport matrices and applies to 3D (6D phase space) general coupled 
lattice without the assumption of a small synchrotron tune. Concerning the evaluation 
of equilibrium beam parameters in an electron storage ring, SLIM can be viewed 
as a method of solving linear Fokker-Planck equation [11, 12] without adopting 
the adiabatic approximation. Therefore, it can account for the variation of one-turn 
map around the ring. In other words, the contribution of diffusion and damping, 
for example the quantum excitation and radiation damping, to the eigen emittances 
depends on the local one-turn map. 

The three eigen emittances €z of a particle beam, with k = J, I, IIT, are defined 
as the positive eigenvalues of i ZS, where i is the imaginary unit, © = (XX7) are 
the second moments of the beam and 


S= (2.4) 


0 
0 
0 
0 


The eigen emittances are invariants with respect to a linear symplectic transportation 
T, as 
iLnewS = iT Eou TS = TCi LoaS)T |, (2.5) 


in which the last step has invoked the symplecticity of T, i.e., TST = S. There- 
fore, i LpewS is related to i XoigS by a similarity transform, thus having the same 
eigenvalues. 

In an electron storage ring, the equilibrium state is a balance between quantum 
excitation and radiation damping. According to SLIM [9], the equilibrium eigen 
emittances are given by 


_ cry f JExs(s)/? 


E ds, 2.6 
nor i 


CK 
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and the second moments of the beam are 


Sy=2 J` eRe[EvE;\], (2.7) 
k=1,11,111 


where a, are the damping constants of the three eigen modes, Cy = 55r.h/ 
(48./3m_.), with fi the reduced Planck’s constant, re the electron classical radius, 
Re[] means taking the real part of a complex number, * means complex conjugate, 
and E; are eigenvectors of the 6 x 6 symplectic one-turn map, satisfying the follow- 
ing normalization condition 


; i, k= 1,11, III, 
ESE =] ”. (2.8) 
fp Sarr TIT 


and EÍ SE; = O fork Æ j, in which * means complex conjugate transpose. Ex; is the 
i-th component of the eigenvector Ex. 

To simplify the discussion, here we only consider the horizontal and longitudinal 
dimensions and use the state vector X = (x, HZS 5)". Under the assumptions that 
the ring is planar x-y uncoupled and the RFs are placed at dispersion-free locations, 
which is the typical setup for present synchrotron radiation sources, the betatron coor- 
dinate Xz = BX can be introduced to parametrize the transfer matrix in a diagonal 
form, with the dispersion matrix given by 


1 0 0O-D, 
0 1 O0-D, 
sala D. -D,1 0 G 
0 0 0 1 
The one-turn map M of X is related to the one-turn map Mg of Xx by 
M = B 'M,B, (2.10) 
with 
_ [Mae 0 
E ( 0 Ma) l 
(2.11) 


M _ (608 ®, ,+a,,, sin ®,, By,z sin By, 

428 —/x,z sin By, cos yz — ay, sin Oyz)? 
in which ®, = 27v, and ®, = 2x v, are the betatron and synchrotron phase advance 
per turn. The eigenvectors of Mg can be expressed using the Courant-Snyder func- 
tions as 
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J 0 
1 Ar 7 7 1 0 


"= Ka e€? v, B| VE elem, (2.12) 
0 J: 


where ®; and ®77; are phase factors which do not affect the normalization of eigen- 
vector and the calculation of physical quantities. Therefore, the eigenvectors of M 
are 


i-a, 
a ve 
1 = : 1 Poe p! : 
E; = Blv, = — VPs e?ii, E =B lv, = =] vB * | eén, 

I x V2 —/Bx Di. + at Dx III Z V2 JE 
0 i=; 

VB: 

(2.13) 


According to SLIM, the equilibrium horizontal and longitudinal emittance are then 


55 arty? [f Hy(s) 


e = (J) = S, 
Mea a a 
saig 55 ary’ f B(s) l 
Z = z = S, 
96/3 ar lo(s)’ 
in which 
7 = EZD + [ara -DD + Be’ -DD 
r ies 2 2 $ 
: f > (2.15) 
PPE (z -Dx -— Dx") + [a-(z — Dix — Dyx')+ 6-5] 
s 2B, 


are the horizontal and longitudinal action of a particle, and () here means particle 
ensemble average, œ y and a, are the horizontal and longitudinal damping constants, 


Uo Uo 
= —(1- D), = — (2+ 9D), 2.16 
aH Ey ), OL DE +D) (2.16) 
$ (1-21) Dx ds 
where Up is the radiation energy loss of a particle per turn, D = a 
D 2 S 
with n = = the field gradient index, ar = a is the fine-structure constant, 


Ae = àe /27 = 386 fm is the reduced Compton wavelength of electron and Hy = 


ve + 2a, Dy D’. + BD? — DDD) is the horizontal chromatic func- 
tion. Therefore, it is the longitudinal beta function B, at the bending magnets that 
matters in determining the contribution of quantum excitation to the longitudinal 
emittance €,. A physical picture is given in Fig. 2.2 to help better understand this 
argument. 
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Fig. 2.2 A physical picture 5 Phase slippage 
to explain why a larger e. 
longitudinal beta function 8; | 

means a larger contribution Quantum | © — 7 =e ji ----------- 
to longitudinal emittance €z, excitation | $ ---------Y------- 
with a given strength of - 
quantum excitation 


Small 8; 


After getting the equilibrium eigen emittances, we can obtain the second moments 
of the beam according to Eqs. (2.7) and (2.13), more specifically, 


€x By —€ ay O 0 
—ExAy Ex Vx 0 0 


= TX 
Xp = (XpXz ) = 0 0 eb. <e,0, | (2.17) 
(0) 0O —€,0, €:y; 
and 
Ey È 
E = (kX) = B! Z; (Boy = ( ZE HEL), 2.18 
( ) a ( ) Xi; Ez ( ) 
where 
Z= Exx + ey: D? =6,0, + Ez Yz Dx D! 
4 (ertr + Ey DD) eye + ey D? J’ 
Noe — Ex (œx Dx +r B,D’) = €z Dx EzYz Dx 
MAL = ( Ex (Yx Dx + a, D!) = €,a,D’, ey: D! , a 


— Ezz Ez Yz 


D= (07% + €,B; Ps i 


The distribution of a Gaussian beam is related to the second moments matrix of the 
beam according to 


Y(X) = 2 exp ( XTE x) = : exp ( Js £) 
~ ayva \ 2 aree ee)’ 


2.1.2 Classical o, « ./\n| Scaling 


Now we first reproduce the classical o, « ./|n| scaling using this longitudinal 
Courant-Snyder parameterization. To simplify the discussion further, in this section 
and the following we focus on the longitudinal dimension only and the state vector 
X = (z, 5)’ is used. We treat first the case where there is only one RF placed at a 
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dispersion-free location. In this case, the linear longitudinal one-turn map observed 
in the middle of the RF cavity is 


10\ (1 —nCo\ (1 0 1—4nCy  —nCo ) 
M=(, = 2.21 
(: ) ({ 1 he ) Care i= tG ean 


with h = e Vrrkrr cos ġs/ Eo quantifying the RF acceleration gradient, where e is the 
elementary charge, Veg is the RF voltage, krr = 27r /Àrp is the RF wavenumber, ¢, 
is the synchronous phase and Eg = ymec” is the electron energy. The Rs, = —nCo, 
a measure for the dependence of z on ô, of the ring and the RF kick h can be viewed 
as the longitudinal drift space and quadrupole, in correspondence to their transverse 
counterparts, respectively. Note however that as mentioned before, the R56 can be 
either positive or negative, while the physical length of a drift space is always positive. 
The linear stability requires that 


(= He 
a0 


<130<hnCy <4. (2.22) 


Actually if the sinousidual modulation waveform is taken into account, the longitudi- 
nal dynamics is more accurately modeled by a standard kick map [13]. We then need 
hnCo be small enough to avoid strong chaotic dynamics. An empirical safe criterion 
is thatO < AnCo < 0.1. For rings working in the longitudinal weak focusing regime, 
|vs| < 1, we then have 


h p? —JhnCo if n > 0, 
-fazen Eon] ee (2.23) 


2 VhnCo if n < 0. 
Therefore the longitudinal beta function £, at the RF center is 


Mi. —7Co nCo 
ae b, d, h a 


In this dissertation, we use the subscript 5 to denote results which are the same with 
that obtained in Sands’ classical analysis [5], although the method used here to get 
these results is different from that of Sands. As |v;| « 1, therefore 


bzs > | — Col. (2.25) 


We will see later in Sect. 2.1.6 that in a longitudinal strong focusing ring, |v,| can be 
close to or even larger than 1, and 8, can then be the same level of or smaller than 
| — Col. 

Using this zs to represent B, of the whole ring, we then get the longitudinal 
emittance obtained in Sands’ analysis 
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55 apa2y? 
és = APAY Pes t -d (2.26) 
96/3 a |p(s)| 


For a ring consisting of isomagnets, p is a positive constant and 


Uo QnAaary? 
— Js = J; 5 2.27 
QL 2Eo ( ) 


with J, = 2 + D the longitudinal damping partition number [5] and nominally J, ~ 
2, we have 


B:s = 1Co 9 /nCo 
Ss sin ®, h’ 
1 2 n VhnCo 
vs = — arcsin x , 
In| 2x 


yr | C 
ozs = y EzsÊzs © me l x oss pzs, (2.28) 


J EzS | Ca y? 
Oss = yV EzSYzs © Aa Tr, 
Bes Js p 


Es © na 1 G05 © oe B 
zS J, p h zS965S s PzS> 
where C} = ae = 3.8319 x 107! m. Therefore, to generate short bunches in an 


electron storage ring, we need to implement a quasi-isochronous lattice, i.e., a small 
n, and a high RF acceleration gradient, i.e., a large h. We also note that the energy 
spread of an electron beam in the classical analysis is dominantly determined by the 
beam energy and bending radius of the bending magnets, and has little dependence 
on the bunch length or global phase slippage of the ring. 


2.1.3 Beyond the Classical o; « y |n| Scaling 


2.1.3.1 Analysis 


Using a single 6,5 to represent that of the whole ring is valid in usual rings where the 
relative variation of £; is negligible and the electron distribution in the longitudinal 
phase space is always upright. But when the global phase slippage is small, the partial 
phase slippage can be significantly larger than the global one and the variation of 
B, and beam orientation in the longitudinal phase space around the ring can be 
significant, thus the classical o; x ./|n| scaling fails. Now we present an accurate 
analysis of this effect using the longitudinal Courant-Snyder formalism. 
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If there is only a single RF cavity placed at a dispersion-free location in the ring, 


then at a specific position s;, the ring can be divided into three parts, with their 
longitudinal transfer matrices given by 


T (spr, sj) = € ae aa) l 
PE (; : i (2.29) 


1 —7(s;, G 
T(s;, SRF) = ¢ ms; ~ ’) : 


where Ñ (srr, Sj) + 7(S;, Spr) = n, In the analysis, the RF cavity is assumed to be a 
zero-length one. The one-turn map at s; is then 


M(s;) = T(s;, Spe) T (Spe, Spr) T (rr, sj) 


_ (1=A(sj, sREVACo —NCo + Ñ (sj, SREDA(SRE, SACS (2.30) 
= h 1 — 7 (srr, Sj)hCo ` 


Therefore, 


Mı2(s;)  —nCo + Ñ (sj, SRE)ACSRE, sj)hCE 


(sj = = 2.31 
B.(s;) sin ®, sin ®, ( ) 
Note that 6, is always positive, and 
dp.(s;) _ [i(srr, si) — 1(5;, srr) ] hCo (=e 1 ) 

dsj sin ®, pls) 7y? (2.32) 

D,(s;) 1 

= 2a;(s;) =z he 

p(s;) y 

which is different from the conventional relation dbr = —2qa,,, in transverse dimen- 


sions [4]. 

The first term in the numerator of Eq. (2.31) is the conventional global phase 
slippage. The second term reflects the impact of the partial phase slippage on 6z. 
In usual rings, the second term is much smaller than the first term, therefore 6, 
is almost a constant value around the ring. As mentioned, the classical formulas of 
bunch length ozs, energy spread oss, and longitudinal emittance €,s in last section are 
actually obtained with such approximation. Now with both terms in the numerator 
of Eq. (2.31) considered, the more accurate formula of the longitudinal emittance is 
then 


(Bz) p 


(2.33) 


— (: y rao ESI SRP D0 = MACS, a | 


n 
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Note that (), here means the radiation-weighted average around the ring, defined as 


P 
_ f pops 
= hor, 
f ors 


i.e., the average is actually conducted at places with nonzero bending field. After 
getting the longitudinal emittance and Courant-Snyder functions, the bunch length 
and energy spread at a specific location s; are then 


(P) p (2.34) 


E€ - 7 hC 
olsi) = y Ezp: (si) © ozs, T if — (Si, SRF) | (SRE, si) 

E 
o5 (si) = y €Y: (si) © oss.) 7 : 

zS 


(2.35) 


We remind the readers that the energy spread and y, are unchanged outside the RF 
cavity. In addition, if the contribution of + is negligible in the definition of n, then 
a, and f; will vary notably only inside the bending magnets. Actually, the chromatic 
H, function, a parameter quantifying the coupling of horizontal emittance to bunch 
length as can be seen from Eq. (2.19), also changes only inside the bending magnets. 
Both arguments reveal the fact that in ultrarelativistic cases, bunch length changes 
only inside the bending magnets. We will see this clearly in Fig. 2.3. 
By investigating the bunch length at the RF cavity 


Ez n E 2 
Oz (SRF) © ozs, | — = osae + (Ñ (sj, SRE))p — n(7(Sj,SRE))pCo, (2.36) 
EzS hCo 
we observe that there exists a lower bunch length limit when 7 approaches zero 


Oz limit = Oss4/ (Ñ (Sj, SRF)) oCo. (2.37) 


This limit is the main consequence of the unavoidable quantum diffusion of longitu- 
dinal coordinate in a storage ring. It has little dependence on the global phase slippage 
and RF voltage, once the beam energy and dispersion function pattern around the 
ring is given. Since ozs « ./|n], the above bunch length limit means £ will diverge 
as ņ approaches zero. The energy spread will thus diverge in this process. 

While the bunch length at the RF cavity will saturate at the limit given by Eq. (2.37) 
with the decrease of 7, the bunch length at other places, from which the partial 
phase slippage to the RF cavity is large, may first decrease and then increase. The 
reason is that the increased energy spread will lead to bunch lengthening through the 
partial phase slippage from the RF cavity to the specific location. In other words, the 
longitudinal beta function ratio between that at the RF cavity and that at the specific 
location may increase with the lowering of 77. 
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2.1.3.2 Experimental Verification 


The above analysis has been confirmed by numerical simulation as presented in 
Ref. [2]. Now we introduce our experimental work on this quantum diffusion effect. 
The experiment was conducted at the Metrology Light Source (MLS) [14—16] of the 
Physikalisch-Technische Bundesanstalt in Berlin. For usual rings, the bunch length 
limit given by Eq. (2.37) is a couple of 10 fs to about 100 fs, while the typical bunch 
length in operation is in 10 ps level. So this effect is negligible in almost all existing 
rings. However, with the accelerator physics and technologies continue to advance, 
more ambitious goals of bunch length can be envisioned and realized in the future 
to benefit more from the electron beam. For example, in an SSMB storage ring, the 
desired bunch length is sub-micron or even nanometer, which corresponds to sub- 
fs in unit of time. The quantum diffusion investigated here then becomes the first 
fundamental issue that needs to be resolved. With such motivation to develop an 
SSMB light source, and considering that it is a fundamental physical effect by itself, 
we believe it is important to experimentally verify this effect. 

To observe the influence of this effect, we need the second term in the bracket 
of Eq. (2.33) to be comparable or larger than 1, which is non-trivial for many of 
the existing storage rings. Other collective and single-particle effects stand in the 
way before arriving at such a small value of n. However, due to the dedicated quasi- 
isochronous lattice design and the individually independent magnet power supplies 
of the MLS storage ring, there is great flexibility in tailoring the lattice optics to obtain 
a locally large and globally small phase slippage simultaneously, thus opening the 
possibility to see this effect in an existing machine. Another characteristic making 
the MLS an ideal test bed of single-particle beam dynamical effects is that it can 
operate with a beam current ranging from 1 pA (a single electron) to 200 mA. 

We have prepared two quasi-isochronous lattice optics at the MLS, named lattice 
A and B, respectively. Lattice A is the standard quasi-isochronous lattice, while 
lattice B is developed and dedicated for this experiment. The optical functions of 
the two lattices are shown in Fig. 2.3. Other related parameters of the two lattices 
are given in Table 2.1. The key difference of these two lattices is that lattice B has a 
much larger partial phase slippage and average value of (8-.),. Therefore, the bunch 
length limit in lattice B (469 fs at 630 MeV) due to this quantum diffusion is larger 
than that in lattice A (115 fs at 630 MeV). Note that with the given parameters set, 
B, in lattice A is almost a constant value around the ring, while £, in lattice B varies 
significantly and at many places is much larger than that in lattice A. 

As can be seen in Fig. 2.3, the magnitudes of horizontal dispersion function Dx 
of lattice B are large at some of the bending magnets, which according to Eq. (2.3) 
means the local phase slippage increases or decreases sharply within them, leading to 
a large variation of local phase slippage 7 and 6,. The small global phase slippage n is 
realized by canceling the contribution of positive and negative D, at different bending 
magnets. We remind the readers that this lattice can also be used for the delayed 
alpha-buckets study in which the momentum differences of particles in different 
alpha-buckets can be translated into large arrival time differences through the large 
partial phase slippage [16], which might be useful for some user experiments. 
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Table 2.1 Parameters of the two lattices of the MLS storage ring used in the experiment 


Parameter Value Description 

Co 48 m Ring circumference 

Eo 630 MeV Beam energy 

Uo 9.14 keV Radiation energy loss 
FRE 500 MHz RF frequency 

VRF 600 kV RF voltage 

hRF 0.01 m7! RF acceleration gradient 
Oss 4.4 x 1074 Classical energy spread 
Ex 197.3 nm Lattice A 

Js 1.95 Lattice A 

(7(8j, SRE)) o 2.5 x 1075 Lattice A 

y (FP (Sj, 5RF)) 1.6 x 1073 Lattice A 

Oz limit 34 wm (115 fs) Lattice A 

Ex 219.4 nm Lattice B 

Js 1.95 Lattice B 

(7(Sj, SRE)) o —5.5 x 1073 Lattice B 

y (FP (Sj, SRE)), 6.7 x 1073 Lattice B 

Oz limit 142 um (469 fs) Lattice B 


To evaluate the possibility of verifying this effect experimentally, the bunch length 
and energy spread evolution around the ring in these two lattices have also been 
presented in Fig. 2.3. Note that the bunch length formula in Eq. (2.35) contains only 
the contribution from longitudinal emittance. Considering the bunch lengthening 
by horizontal emittance at dispersive locations, according to Eq. (2.19), the more 
accurate formula of bunch length is [1, 18] 


0z = Vy €,B, + EHy. (2.38) 


Strictly speaking, Courant-Snyder and dispersion functions are only well-defined in 
a planar uncoupled lattice and only when the RF cavity is placed at a dispersion-free 
location. For a general coupled lattice, the more accurate SLIM formalism should be 
referred, i.e., 


o= [2 >) «lEs, 
kaL ILII 


os= [2 J` exlExel?. 
k=1,11,11T 


(2.39) 
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Fig. 2.3 Two lattices used in 
the experiment. Evolution of 
a x,y, b Dx and Bz, € oz, d 
as, around the ring. In this 
plot, the RF cavity is placed 
at srF = 0 m and Ver = 600 
kV is applied, the global 
phase slippage used is 

n = 1 x 1075. The dipole 
magnets are shown at the top 
as blue rectangles. Each 
dipole has a length of 1.2 m 
and bends the electron 
trajectory for an angle of 
7/4. Px, y and Dy are 
obtained by fitting a model to 
the BPM-corrector response 
matrix (LOCO) [17]. The 
bunch length and energy 
spread evolution are 
calculated based on the 
longitudinal Courant-Snyder 
formalism and SLIM 
formalism 
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On the other hand, although the RF cavity is placed at a dispersive location in lattice 
B, we have confirmed that the Courant-Snyder parametrization for beam dynamics 
analysis in this case is still largely valid, since the difference of result between that 
given by the longitudinal Courant-Snyder formalism and the more accurate SLIM 


formalism is very small. 


As can be seen in Fig. 2.3, in which the global phase slippage 7 is lowered to 
be 1 x 1075, which corresponds to a synchrotron frequency of f, = 2.2 kHz with 
Ver = 600 kV, the energy spread grows to be os = 7.9 x 1074, while the classical 
energy spread is oss = 4.4 x 1074. Such an amount of energy spread growth is 
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detectable by measuring the spectra of Compton-backscattered (CBS) photons from 
the head-on collision between a CO) laser with the electron beam at the MLS [14, 
19]. In addition, the bunch length difference in these two lattices are large enough 
to be observable by evaluating the spectra and power of coherent THz radiation, and 
invoking streak camera measurement. 

To exclude the influence of collective effects, the beam current is lowered to around 
6 wA/bunch in a multi-bunch filling mode in the experiment. There is no indication 
of microwave or other collective instabilities. The beam is stable (no fluctuation of 
radiation source point observed) and its width and energy spread are independent 
of the beam current when the single-bunch current is as low as the value applied in 
the experiment. The horizontal chromaticity has been carefully corrected close to 
zero (about 0.05) to minimize the beam energy widening arising from the betatron 
motion of particles as will be reported in Sect. 3.2. The longitudinal chromaticity has 
also been corrected to a small value to mitigate longitudinal nonlinear dynamics. We 
note that a large quantum diffusion of longitudinal coordinate (a root-mean-square 
value of 0.54 um or 1.8 fs per turn in lattice B at 630 MeV) actually helps suppress 
collective beam instability of ultrahigh frequency, as it will disperse any fine time 
structure in an electron beam like density modulation and energy modulation [7]. 

To get an idea about the bunch length in the two lattices, first we measure the coher- 
ent THz radiation spectra and power as a function of the synchrotron tune in the two 
lattices. The shorter the electron bunch, the higher frequency range the coherent THz 
radiation spectra extends and the larger radiation power we can obtain. In the experi- 
ment, the synchrotron frequency fs, thus the global phase slippage n (fs x /|n]), is 
controlled by slightly changing the quadrupole currents while keeping the dispersion 
function pattern unchanged. The THz beamline has its source point at 7; bending 
angle (s = 38.775 m) at the 7-th dipole, counted from s = 0 m in Fig. 2.3 which is 
where the RF cavity is placed. To get the coherent synchrotron radiation emission 
spectra in the THz spectral range, a commercial, Michelson-type FTIR spectrome- 
ter (Vertex 80v) in combination with a 4K liquid helium cooled composite silicon 
bolometer was used for measuring interferograms. After fast Fourier transform of 
the data, the emitted spectrum can directly be accessed. For this experiment a series 
of 128 interferograms have been acquired and the average Fourier transformed. 

The measured coherent THz radiation power, integrated with wavenumber from 
1 to 20 cm™!, together with the theoretical bunch length at the THz observation 
calculated using Eq. (2.38), are shown in Fig. 2.4a. The measurement results agree 
with our expectation reasonably well. In particular, we notice that in lattice B, the THz 
power first increases and then decreases, with the lowering of the synchrotron tune, 
while the radiation power in lattice A monotonically decreases and then saturates 
in this process. This observation agrees well with our theoretical prediction of the 
bunch length evolution in these two lattices. Not presented here, we also notice that 
the frequency range of the spectra evolves consistently with the integrated power, 
i.e., a larger THz power corresponds to a higher frequency range coverage. To be 
more rigorous, we remind the readers that the bunch length in lattice A at the THz 
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radiation observation point in principle will also diverge, as explained in last section, 
if we push the phase slippage factor of the ring even closer to zero, which in practice 
is a demanding work. 

During the measurement of coherent THz radiation, we at the same time employed 
a streak camera to measure the electron bunch length directly. The streak camera at 
the MLS is installed at the undulator beamline (opposite the RF cavity, s = 24 m 
in Fig. 2.3). For the experiment, the undulator was closed from the “open” gap of 
180 mm to 45.7 mm to have the fundamental-mode undulator radiation at a visible 
wavelength available for the streak camera. The measurement results of bunch length 
and the comparison with theory is presented in Fig. 2.4b. Note that we have shifted 
the measured raw data downwards by 8 ps in the plot. The errorbars in the plot are 
the standard deviation of the fitted results for each single column of the recorded 
streak camera image. Again we observe the significant difference in the two lattices 
concerning the bunch length evolution as a function of the synchrotron tune, which 
agrees qualitatively with theory. However, quantitatively the measured raw data of 
bunch length deviates notably from the theoretic prediction. 
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Realizing that there will be unavoidable systematic errors concerning the streak 
camera measurement because we are close to its resolution limit, we try to use the 
model below to fit the data with the theory, 


AZft = AZZ easur — noise’ — offset, (2.40) 


where Az means the bunch length. Note that here we use the full width at half 
maximum (FWHM), instead of the root mean square, to quantify the bunch length, 
since in the real case, the bunch profile is unavoidable non-Gaussian to some extent, 
especially when the phase slippage factor of the ring is small. The noise in the above 
equation is used to model the square sum-type error, while the offset accounts for the 
systematic shift concerning the measurement results. The fitted data (noise = 7 ps 
and offset = 3 ps applied) agrees well with the theoretical curve as shown in Fig. 2.4b. 
We remind the readers that all the data points in the plot are modeled with the same 
noise and offset. 

Further, we have measured the electron beam energy spread in the two lattices, 
using the head-on CBS between a CO, laser with the electron beam. Note that the 
RF voltage applied in the above bunch length measurements is 500 kV, while now 
it is 600 kV when doing the energy spread measurement. The measurement of CBS 
photon spectra and the evaluation of electron beam energy spread based on it is a 
well-established method implemented at the MLS, and is used in this experiment 
to confirm the energy widening as we push the bunch length close to the limit, 
by lowering the global phase slippage 7. More details about this CBS method can 
be found in Refs. [14, 19]. Quantitative analysis revealing the energy spreads os 
normalized by the classical energy spread oss, and its comparison with the theoretical 
prediction from Eq. (2.35) for the two different lattices are shown in Fig. 2.4c. The 
error bars in Fig. 2.4c are the root-mean-square uncertainties of the measurements 
and are due to calibration errors and counting statistics. The data acquisition time 
of a photon spectrum is 15 min. It can be seen from Fig. 2.4 that in lattice B the 
energy spread grows significantly with the decrease of n, in the figure synchrotron 
frequency f,, to the level of 1 x 1075, while the energy spread stays almost constant 
in lattice A. Again the measurement agrees qualitatively with the theory. 

There is still some deviation of the measured energy widening and the theoretical 
prediction for lattice B. Candidate explanations are: first, there is some uncertainty 
in the determination of synchrotron frequency fs, especially when f, is lowered to 
2 ~ 3 kHz, considering the fact that the peak of the synchrotron frequency spectrum 
then can be as wide as 0.5 kHz; second, there could be some remaining higher-order 
phase slippages which may contribute to the energy spread growth when n is small 
due to its impact on the longitudinal phase space bucket, while the theory assumes a 
linear phase slippage. 

The above presented measurements of bunch length and energy spread are very 
demanding and are moving on the edge of the experimentally accessible parameter 
space. Nevertheless, we see a nice qualitative agreement with the theory presented 
in this section, proofing important experimental evidence to support the theoretical 
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analysis. As far as we know, this type of investigation can actually not be performed 
at any other operating storage ring. However, we recognize that the deviation of 
the quantitative numbers between the measurements and theory concerning both 
the bunch length and energy spread emphasizes the need for an even more improved 
model. Summarizing we state that our experimental work supports the existence of the 
analyzed quantum diffusion effect, and the argument that the quantum excitation on 
longitudinal emittance at a given location depends on the longitudinal beta function 
there. The evidence however is not strong enough to claim this is a fully consistent 
proof of the effect. 


2.1.4 Campbell’s Theorem 


We point out that quantifying the impact of variation of £; around the ring on €; using 
partial phase slippage variance (77) Pi (7); as that done in Ref. [6] and also our 
previous publication Ref. [1], is not generally correct. The reason is that while the 
photon emission process is stochastic, the evolution of partial phase slippage around 
the ring is deterministic. So the diffusion of z each turn d? due to quantum excitation 


is 


2 
d = (z2) = (z)? = CRZ N) (=). (2.41) 
0 
instead of z 
d? = (2°) — (z2) = Ca (7) 0 — (ñ) (N) (=) (2.42) 
0 


as that given in Ref. [6], where 7 is the partial alpha slippage calculated using the 
final observation location as the ending point, (N) is the expected number of emitted 
photons, u is the photon energy, (u?) and later also (u) mean the average is taken 
with respect to the photon energy spectrum. 

This result can be understood with the help of Campbell’s theorem [20]. From 
this theorem some expectation result for the Poisson point process follows. For 
example, for the application in synchrotron radiation, we have 6 = — 0; ae where 
the subscript ; means the i-th photon emission. Then according to Campbell’s theorem 
we have 


; (2.43) 


(8?) — (8)? = (N) (5) = Taya (25), 
E E 


where Ñ is the number of photons emitted per unit time in the dipoles and Tyipole 18 
the total time within dipoles. Equation (2.43) is why N (u?) appears so often in the 
calculation of energy spread, emittance, etc., in electron storage ring physics. Note 
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that the relation in Eq. (2.43) holds as long as the radiation is a Poisson point process. 
It is independent of whether (uv) = 0 or not, and is also independent of the detailed 
spectrum of the photon energy. In other words, the key of a Poisson point process 
is the randomness in whether there is a kick or not, i.e, the kick number, and not in 
the randomness of the size of the kicks. The importance of this theorem for electron 
dynamics was first pointed out by Sands [21]. A proof can be found in the article of 
Rice [22] and a less rigorous but simpler one in the lecture note of Jowett [23]. 

Now we can understand Eq. (2.41) as follows. Suppose that the RF is our observa- 
tion point. We divide the ring into many sections, and in each section 7(spr, $j) does 
not change much. Then the change of electron longitudinal coordinate in one turn is 
z= pe Zj, with z; = >>; Cor(srr, Sji) r the contribution due to photon emissions 
within the section j. According to Campbell’s theorem, the variance of z; is 


` . Ju? 
Var(zj) = C37 (sre, spin), (2.44) 
0 


where ¢; is the time within the dipoles in section j. As the photon emissions in 
different sections are uncorrelated, then the variance of z is the sum of variance of 
Zj 


mth 
| n° (SRF, Sp); ; 2 
a [ii RFs? j draid) 


(2) — (2) = CG 
9 Tiotal 


: 0 (2.45) 


2/2 u 
= C (7 (srr, $j) p(N) (=). 
Fo 
in which Tiota = >>; t j is the total time within the dipoles. So now we have obtained 
Eq. (2.41) following Campbell’s theorem. 

We can also view the above argument from another way. Given the same dispersion 
function pattern, which means the same (ñ?) a (ny? as it is independent of the 
observation point, a different longitudinal beta function pattern can be generated if 
the RF is placed at a different location, therefore resulting in a different longitudinal 
emittance according to Eq. (2.14). 

Changing the RF location means shifting 7(s;, srr) up or down as a whole. 
According to Eq. (2.33), the equilibrium emittance is a parabolic function of the 
shifted value. When the RF is placed at a location such that (7(s;, Spr)» = a we 
arrive at the minimum longitudinal emittance 


2 
€z,min = €zS (: + ci) : (2.46) 


The maximum longitudinal emittance is realized when the RF is placed at a place such 
that | (ins i> SRF) p — 7| reaches the maximum possible value. When the minimum 
longitudinal emittance is reached, the bunch length at the RF is 
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C 2 
Oz, min (SRF) = ou) + | ai (7)? = (2) | Cc. (2.47) 


In the case of ultrasmall 7, we have 


Oz min (SRF) X O55 (Ao _ (3Co, (2.48) 
and ( g 2) 7 
Es (7°) p — (ny5) Co 
z,min a ——— C h + 
eie aT ( a ViCo/h (2.49) 


> Lossy) (i?) — ()2.Co. 


The equality holds when £;s = ,/ eo = Co,/ (i?) — (71)2.. Therefore, the variance 


of partial phase slippage can be viewed as a parameter to quantify the lowest pos- 
sible contribution of this effect to the equilibrium bunch length at the RF and the 
longitudinal emittance with a dispersion function pattern given, if we can choose 
the location of the RF as we want. However, in a real machine, the RF location is 
fixed, and Eqs. (2.33) and (2.35) should be referred. This is why we state that using 
(7°) = y to quantify the impact of this effect is not generally correct. 


2.1.5 Minimizing Longitudinal Emittance 


It is clear that the quantum diffusion of z needs to be carefully treated for the realiza- 
tion and long-term maintenance of ultrashort bunch or small longitudinal emittance 
in either a multi-pass device or a single-pass transport line with bending magnets and 
large dispersion. A lower operating energy is preferred for suppressing the strength 
of quantum excitation. Note that the energy scaling laws of this effect are differ- 
ent in the one-turn or single-pass and steady-state cases; for the single-pass case, 
i.e., Eq. (2.41), the root-mean-square diffusion of longitudinal coordinate d; « y?°, 
while for the steady-state case, i.e., Eq. (2.48), Oz min X y, because the radiation 
damping time also depends on y. 

As can seen in Eq. (2.14), the longitudinal beta function £, with respect to the 
longitudinal dimension plays a role similar to that of the chromatic function Hy in the 
transverse dimension. As both the longitudinal and transverse emittances originate 
from quantum excitation, for a ring consisting of identical isochronous cells, the 
same scaling law of the theoretical minimum emittance (TME), .e., €:... TME & y?0? ‘ 
concerning the beam energy and bending angle of the dipole can be expected. Note 
that the TME is independent of the bending radius. But we will show soon that the 
bunch length limit does depend on the bending radius. According to the scaling, a 
ring consisting of a larger number of isochronous cells, each with a smaller bending 
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V 


Fig. 2.5 A symmetric dispersion function pattern which makes each half of the bending magnet 
isochronous, which is desired in minimizing the theoretical minimum bunch length and longitudinal 
emittance when there is only a single RF or laser modulator in the ring. The lattice design for the 
realization of the dispersion function pattern can be found in Ref. [25] 


angle, can better minimize the emittance than a ring consisting of fewer cells with 
larger bending. Generally, it is easier to realize small emittance in a larger ring. 

Equation (2.48) gives the lower bunch length limit by optimizing the location of 
RF cavity, with a given dispersion function pattern. To make this limit as small as 
possible, in addition to ensuring a small global phase slippage, the variation in the 
partial phase slippage should also be well confined by means of dedicated lattice 
design. More specifically, the strategy is to tailor the horizontal dispersion func- 
tion, thus to minimize the longitudinal beta function at the bending magnets. At the 
MLS, the small global phase slippage is achieved by means of an overall integration 
cancellation between the large positive and large negative horizontal dispersions at 
different dipoles [15, 16]. Therefore, the partial phase slippage varies sharply within 
the dipoles, leading to a large partial phase slippage variation and significant quan- 
tum diffusion of z. To obtain small global and partial phase slippages simultaneously, 
such cancellation should be done as locally as possible, and the magnitudes of the 
dispersion at the dipoles should also be minimized, thus making the partial phase 
slippage vary as gently as possible. In other words, each partial component of the 
ring should be made as isochronous as possible. In this sense, the dispersion function 
pattern in Fig. 2.5 is the most locally isochronous bending magnet [24], i.e., each 
half of the bending magnet is isochornous. 


2.1.5.1 Constant Bending Radius 


Now we present some quantitative analysis of the minimization of longitudinal emit- 
tance. In this section we use the partial R56, defined as 


i %/D) 1 
AEE E E f ( 5) ds, (2.50) 
i p(s) y 
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for the analysis. As can be seen Eqs. (2.48) and (2.49), it is \/(F*), — (F)? that 
determines the theoretical minimum bunch length and longitudinal emittance. Here 
we present an analysis of the above scaling by evaluating ,/(F?), — (F 2 of a single 


isochronous bending magnet with constant bending-radius, and whose dispersion 
function is symmetric with respect to the dipole middle point as shown in Fig. 2.5. 
To simplify the analysis, we start from the middle point of the dipole where the 
dispersion angle is zero D/,, = 0, then the dispersion as a function of angle ¢ is 


D(¢) = Dycos¢ + p(1 — cos ġ), (2.51) 


where Dm is the dispersion at the middle of the dipole. As we are mainly inter- 
ested in the relativistic case, ignoring the contribution of a on F, the condition of 
isochronicity of each half of the dipole is 


f pido, (2.52) 
0 


with 6 the bending angle of each dipole. Substituting Eq. (2.51) into Eq. (2.52), we 
get 


e 1 
Dm =p (1-7) * -= 0”, 
o( 5) 24? 


8 1 
D: =p |1- — | © =°, 
$ tan $ 2’ 


where D, is the dispersion at the entrance and exit of the dipole. Then we have 


(2.53) 


n 
n> 


$ f 
FÆ) = [ D(B)dB = p -g -2 sin s). 


1 fo 
F= | Fod =O, 


(F?) E 2 s+ (2) ag 7 p (8) 
P 0 =! 6 2 tan $ sin (£) ’ 


_ v210 
ee 2520” 


(F?)p 

(2.54) 

Note that the F in this section is defined with the middle point of the bending magnet 

as the starting point. Therefore, for a ring consisting of such isochronous isomagnets 

(note that the global phase slippage of the ring is non-zero for a stable beam motion), 
we have 
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J/210 
Oz min © O58, (F*) 9 — JE- F} = =S G Jer 6? x Spy}, (2.55) 


and 


/ V210 Cy 2 3 x 293 
€z min © 20A (F?)p F5 = 1260 J,” yo xy 0. (2.56) 


We can also derive the above emittance scaling using directly the longitudinal 
beta function evolution in the dipole 


B.(@) = Ban — 22m F (Ø) + Yan F” (0), (2.57) 


where 2m, Bzm, Yzm are the longitudinal Courant-Snyder functions at the middle of 
the dipole. Then 


1 
(Bz) = ala a (ddd 
P? (a2, + 1) [2 (6? — 24) + 30 (0 + 3 sin(@)) csc? ($)] 


2.58 
TAB... + Ban (2.98) 


n 28° (im + 1) 


30240B-m + Bam: 


The minimum average of 6z, (Bz) ,min, thus the minimum longitudinal emittance, 
€z min, 1S realize when 


(Bz) p,min N 210 3 


zm = 0, Ban = = 0”. 2. 
T= Pe 2 2520 ” paa 
The corresponding minimum longitudinal emittance is 

55 apy’ 27 (Bz)p min _ V210 C3 20°, (2.60) 


Se ea oat P —«:1260 J 


which is the same as that given in Eq. (2.56). If the longitudinal damping partition 
number J, = 2, then for practical use we have 


o; minium] ~ 4.932 [m] Eo [Ge V]6 [rad], 
Ez minlnm] ~ 8.44E5[GeV]6° [rad]. 


(2.61) 


This is the main result of our analysis of theoretical minimum bunch length and 
longitudinal emittance in a longitudinal weak focusing ring. A comprehensive study 
of minimizing longitudinal emittance can also be found in Ref. [4]. 
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Actually if we relax the condition of isochronousity for each dipole, and make 
Dm and Df, also variables that can be optimized, then we have 


1 fo 
(Bz) = F b: (ġ)do 


9 
2 


D20? 4D,,p0* — 10p76° + 504D29264 — 15 D2 p266 


= Bom + 126, * 240 Bim 16128087 4 
(2.62) 
When Dm = — ae Di, = 0, and Bum = Loo, we reach the minimal of (6z) p,min = 
M 963 and longitudinal emittance 
Ez,min = Si 203, (2.63) 


420 J,” 


This is also the result Eq. (29) in Ref. [4]. The emittance given in Eq. (2.63) is smaller 
than that given in Eq. (2.60) or equivalently Eq. (2.61). However, we should note that 
when isochronousity of each dipole is broken, to make the longitudinal beta function 
in different dipoles identical, we need RF kick between each two neighboring dipoles 
to adjust œ, there, or we need a very long drift space between them if we consider the 
contribution of +, on Rs¢. This means at least N RFs or laser modulators are needed 
if there are N dipoles in the ring. This is not very feasible in reality. In addition, 
placing RFs at dispersive locations will make the dynamics becomes transverse- 
longitudinal coupled. In the more-confronted or practical case of a single RF in the 
ring, Eq. (2.61) is a more self-consistent evaluation of the theoretical minimum bunch 
length and longitudinal emittance. 

Since high-power EUV radiation is of particular interest for EUV lithography [26], 
let us now do some evaluation based on our investigations to see if we can realize 
high-power EUV radiation in a longitudinal weak focusing SSMB storage ring. For 
coherent 13.5 nm EUV radiation generation, we need an electron bunch length around 
3 nm or shorter. The lower limit of bunch length o; min should be smaller than this 
desired bunch length to avoid significant energy widening. Here we assume o; min < 
2 nm. If Ey = 400 MeV and p = 4 m (B = 0.334 T), then oss = 1.7 x 1074. To 
realize oz min < 2 nm, according to Eq. (2.61), we need 6 < 0.0797 rad ~ ae 
which means 79 bending magnets are required in the ring. If the length of each 
isochronous cell with a single bending magnet can be designed to be around 2 m, 
then the circumference of the ring can be about 180 m, considering the sections of 
laser modulation, radiation generation and the energy supply system, etc. 

Applying the lowest phase slippage factor realizable in practice at present, whichis 


about n = 1 x 1076, to realize a3 nm bunch length in such a ring, the required energy 
os 


2 
chirp strength is then h = nCo (=) = 5.88 x 10° m™!. The effective modulation 


voltage of a laser modulator using a planar undulator is related to the laser and 
undulator parameters as [27] 
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JJ]K |4PLZoZ Lu 
V, = Ed LOOR tana! (2.64) 
y ÀL 2ZrR 


in which [JJ] = Jo) — J\(x) with J, the n-th order Bessel function of the 
first kind and x = oR ,K= = kT = 0.934 - Bo[T] - A, [cm] is the dimensionless 
undulator parameter, determined by the undulator period and magnetic flux density, 
P; is the modulation laser power, Zo = 376.73 Q is the impedance of free space, 
Zr is the Rayleigh length of the laser, L,, is the undulator length. The linear energy 


chirp strength around zero-crossing phase is therefore 


V JJIK [4P,ZZ L, 
ia pa A LOOFR tan! kt, (2.65) 
Eo y2mc? AL 2ZR 


where k; = 27 /àÀ;z is the wavenumber of the modulation laser. If àz = 270 nm, A, = 
4cm, Bo = 1.02 T, L, = 1.6m, Zr = Eu , then the required modulation laser power 
to get h = 5.88 x 10° m7! is P; = 2.75 GW. At present, 1 MW stored average laser 
power is the state-of-art value we can realize in practice for an optical enhancement 
cavity. So we can only operate the cavity in pulsed mode, which means the average 
radiation power will be limited. To put it another way, if a CW optical cavity and 
a practical global phase slippage is applied, a longitudinal weak focusing SSMB 
storage ring can only realize bunch length as low as tens of nanometer with a beam 
energy of several hundred MeV. Such an SSMB ring can provide high-power radiation 
with wavelength Àg Z 100 nm. We have presented in Table 6.1 of the final chapter 
an example parameters set of SSMB ring for high-power infrared radiation. We 
remind the readers again that considering the nonlinear modulation waveform, we 
actually need 0 < hnCo < 0.1 to avoid strong chaotic dynamics in a longitudinal 
weak focusing ring. 

The above limitation of longitudinal weak focusing scheme is the motivation for 
us to develop the longitudinal strong focusing SSMB and transverse-longitudinal 
coupling SSMB, or generalized longitudinal strong focusing SSMB, to compress 
the bunch length further for coherent EUV and soft X-ray radiation generation. We 
will present the details of these advanced scenarios in the following part of this 
dissertation. 


2.1.5.2 Transverse Gradient Bends 


The above analysis of theoretical minimum bunch length and emittance is for a 
constant bending radius. To minimize the longitudinal emittance further, transverse 
and longitudinal gradient bending magnets (TGB and LGB) can be invoked. Below 
we conduct some calculations based on the similar dispersion configuration as shown 
in Fig. 2.5, but this time using a TGB. The Hill’s equation for the dispersion is 
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d? D(s) 1 1 
r (= ko) D(s) =< (2.66) 


For simplicity, here we only investigate the case of a constant bending radius p(s) = 
p and a constant transverse gradient k(s) = k. To simplify the writing, we denote 


1 
Sa k. (2.67) 
p 


If g > 0, then the solution of Eq. (2.66) is 


D(s) = D; cos (./gs) + D; sin (V5) + : [1 — cos (/gs)], 
ve one (2.68) 
D' (s) = —D,,/g sin (ygs) + D; cos (ygs) + Ta sin (ygs) , 


where Dj; and D; are the initial dispersion and dispersion angle at the origin s = 0 
m, respectively. If g < 0, then the solution of Eq. (2.66) is 


D(s) =D h( igi eg IS) a o Te] 
s) = D; cos gls) i [ cos ( g Jl 
vig! Isle 
1 
D'(s) = -Di yg] sinh (/igis) + D{ cosh (/igis) + — sinh (/igis). 
viglo 

(2.69) 

Below, we present the derivations for the case of g > 0 and the results are similar 


when g < 0. Like our previous calculations, we set the origin at the middle of the 
dipole where D/, = 0, the dispersion as a function of angle ¢ is then 


: [1 — cos (./gp¢)]- (2.70) 


D = Dm + — 
($) cos (/gp¢) a 


Substitute Eq. (2.70) into the isochronicity condition Eq. (2.52), we get 


1 2 1 7 
Dn = f v8P2 l ~ —— p0? (: + aie). 


i a 
gp sin (/8p3 24 240 
(vVEr3) l (2.71) 
6 1 5 1 1 
p.= D( )- 1 a N pe? (14+ Zep?) 
2 gp tan (,/gp5) 12 60 


The ,/(F?), — (F)? in this case can be calculated as follows 
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8 
2 


sin (V803) 


$ 1 
rø=f Ddo’ = E sin (ve). 
0 &P 


(E= |, Pode =0 


6gp" tan (V82 (3)) sin (YB ($))° 
o v210 3 gp°07 
(Fa a=), ~ 77500 °° (1+ 40 ) 


(2.72) 
For example, to reduce ,/{F), — (F y by a factor of two compared to the case of 


no transverse gradient, we need 


gee 20 


Te cd se os 


(2.73) 


For the example shown in last section, pọ = 4m, 0 = m, then k = 63 m~?, which is 
a practical gradient. 


Besides the influence on ,/(F?), — (F i the transverse gradient may also affect 


the damping partition and hence has an impact on the bunch length and longitudinal 
emittance. For the specific case of a constant bending radius with a constant transverse 
gradient we are treating, we have 


1 27 
h = ¢ —ds = —, 
p p 


(2.74) 
Dy 3 
l= f al + 2p?k)ds = 0, 


where J and J, are the radiation integrals [5]. Then J, = 2 + a = 2. So a dipole 
with a constant bending radius and a constant transverse gradient is not very flexible 
in controlling the damping partition number, due to the constraint of isochronous 
condition. A varying transverse gradient may be helpful to minimize the longitudinal 
emittance, and optimization of the transverse gradient profile based on numerical 
method can be invoked. The application of TGB can also be analyzed following the 
same formalism, which we do not detail in this dissertation. 
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2.1.5.3 Transverse Emittance Scaling 


For completeness, now we present the horizontal emittance scaling in a longitudinal 
weak focusing SSMB storage ring. To realize the dispersion function pattern shown 
in Fig. 2.5, in thin-lens approximation, the horizontal optical functions at the dipole 
middle point are correspondingly! 


p0 sin By, 
xm = 0, (2.75) 


Pom = ~ 3 14cos 0,’ se 


with ®, the betatron phase advance per isochronous cell, which usually lies in 
(x, 277). We have assumed there is only a single dipole each isochronous cell. The 
normalized eigenvector corresponding to the horizontal plane at the dipole middle 


point is 
V Pxm 


i 


oe ee, id, 
E,(0) = Zl o dd, (2.76) 
Jian Dm 
0 


The transfer matrix of a sector dipole with no transverse gradient is 


cosa p sina 000 p(l —cosa@) 
ane cosa 00 0 sina 
0 0 1 pa 0 0 
— sing —p(1—cosa)0 0 |p (2 +a- sina) 
0 0 00 0 1 


Then 


H, (a) = 2|E75(a)|* = 2|S(@)E;5(0)|? 


2 
+ cosa) : 


dap 0 sin ®, ie a 3 1 + cos ®, E 
31+ cos ®, 6 sin ®, 24 
(2.78) 
Note that H,.(—a) = H, (œ). Putting in «y ~% = x A = iC, E, then according 


to Eq. (2.14), the equilibrium horizontal emittance is 


' Private communication with Zhilong Pan. 
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NIS 
z 
os, 
g 


55 apà?y’ 2r f 

€y = — = x2x 
96/3 1c, 4 0 0 

Z= p. 


55 Vy ®, 1 0, 
27 oe a tan ( ‘+ cot ( ) 63, 
24/3 CE L72 2 80 2 


Putting in the numbers, we have 


(2.79) 


2 3 1 ®, 1 ®, 
€,[nm] = —366.5 Eo[GeV]@° [rad] 9 tan 5) + 10 cot 7 ‘ (2.80) 


Note that the horizontal emittance diverges as ®, approaches x or 277. 


2.1.6 Longitudinal Strong Focusing 


2.1.6.1 Analysis 


The analysis in the above sections considers the case with only a single RF. When 
there are multiple RFs, for the longitudinal dynamics, it is similar to implement mul- 
tiple quadrupoles in the transverse dimension, and the beam dynamics can have more 
possibilities. Longitudinal strong focusing scheme for example can be invoked [8, 
28], not unlike its transverse counterpart which is the foundation of modern high- 
energy accelerators [29, 30]. The linear beam dynamics with multiple RFs can be 
treated using SLIM the same way as that with a single RF. When all the RFs are placed 
at dispersion-free locations, the Courant-Snyder parametrization can be applied as 
analyzed in previous sections. Here we use a setup with two RFs as an example to 
show the scheme of manipulating £, around the ring. The schematic layout of the 
ring is shown in Fig. 2.6. The treatment of cases with more RFs is similar. 
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Fig. 2.6 A schematic layout of a storage ring using two RF systems for longitudinal strong focusing 
and an example beam distribution evolution in the longitudinal phase space. Note that the tilted 
angles of the beam distribution and bunch length ratios at different places do not strictly correspond 
to the parameters in Table 2.2, but only to present the qualitative characteristics 
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We divide the ring into five sections, i.e., three longitudinal drifts (R56) and two 
longitudinal quadrupole kicks (A), with the transfer matrices given by 


(RẸ _ (10 ft Re 
To = (j l , Tre = hy 1 , Tp = 0 i ’ 


10 1 RO ey 
Tre = E D , Tp = f x | ‘ 
Then the one-turn map at the radiator center is 
Mr = Tps Trr2Tp2Trr: Tp1. (2.82) 


Linear stability requires that |Tr (Mgr) | < 2, where Tr() means the trace of. For 
the generation of coherent radiation, we usually want the bunch length to reach its 
minimum at the radiator, then we need œ, = 0 for Mr. 

With the primary purpose to present the principle, instead of a detailed design, 
here for simplicity we only discuss one special case: RY = RẸ, hı = h = h. The 
treatment of more general cases with different signs and magnitudes of R and RË 
and hı and hz is similar, but the same-signed RP and Re might be easier for a 
real lattice to fulfill. For example if RY, Re > 0, a possible realization of them are 
chicanes. 

For the special case of RY? = Re hı = hy = handdenotet; = 1+ Reh, h = 
2+ RPh, we then have 


togi 
Mp = (6182-1 “FE |, 2.83 
i i) eo 


The linear stability requires |¢)¢2 — 1| < 1, and the synchrotron tune is 


2r 


a — 1] if 82281 
e arccos [¢ı¢2 — 1] if ais (2.84) 
1 — + arccos [ġġ — 1] if += < 0. 


Here we give one example parameter set with a stable linear motion as shown in 
Table 2.2. According to the longitudinal Courant-Snyder functions given in Table 2.2 
(note the values of £, and the signs of œ), the evolution of electron distribution in the 
longitudinal phase space around the ring (note the bunch lengths and orientations) is 
qualitatively shown in Fig. 2.6. If we can realize €, < 5 pm in such a strong focusing 
ring, then we have 0; (Sraa) = /€zBz(Sraa) < 3 nm. We remind the readers that the 
contribution of modulators to longitudinal emittance should be carefully counted in a 
longitudinal strong focusing SSMB storage ring. Note that the energy chirp strength 
needed here is one order of magnitude smaller than the example of using a single RF 
or laser modulator to realize 3 nm as discussed just now. This benefit originates from a 
much compressed £, at the radiator in a strong focusing ring. However, we recognize 
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Table 2.2 An example 


. Parameter Value 
parameters set corresponding =m Eo ee 
to the setup shown in Fig. 2.6. Rsg 15 pm 
The subscripts —/+ means RO —100 um 
right in front and after the R p 15 um 
corresponding element h Sxi ne 

1 0.25 

o2 7 

Vs 0.115 
Bz (Sraa) 1.9 um 
Bz (SRF) 121 ym 
Oz (Stad) 0 

az (SRFI-) -7.9 

2 (SRF1+) =1.9 

tz (SRF2-) 1.9 

az (SRF2+) 7.9 


that the laser power needed (20 MW level if 270 nm-wavelength laser is applied) 
is still demanding, and here our primary goal is to present the principle based on 
which the interested readers can choose and optimize the parameters for their target 
applications. We will discuss in Chap. 3 the application of transverse-longitudinal 
coupling scheme to lower the requirement on the modulation laser power further, 
to make the optical cavity can be operated in CW mode, thus to improve the filling 
factor of electron beam in the ring and the average output radiation power. 


2.1.6.2 Discussions 


Here we make several observations from the above analysis and numerical example, 
which we believe are important. First, 6, in a longitudinal strong focusing ring 
can be at the same level of or even smaller than the ring | R55 = —nCo|, while in a 
longitudinal weak focusing ring 6; >> | — nCo|. Therefore, the bunch length can thus 
be much smaller than that in a longitudinal weak focusing ring. This is the reason 
behind the application of longitudinal strong focusing in SSMB to realize extreme 
short bunches [28, 31]. We remind the readers that the longitudinal emittance of 
electron beam in a longitudinal strong focusing ring still cannot be smaller than that 
given in Eq. (2.63), due to the intrinsic partial phase slippage, thus the evolution of 
longitudinal beta function, in a dipole. 

Second, £; changes significantly around the ring in the longitudinal strong focus- 
ing regime. Therefore, the bunch length and beam orientation in the longitudinal 
phase space varies greatly around the ring, as shown qualitatively in Fig. 2.6. This 
means the adiabatic approximation cannot be applied for the longitudinal dimension 
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anymore. Actually the adiabatic approximation also breaks down in the case corre- 
sponds to Fig. 2.3, where the change of £; around the ring is significant although the 
total synchrotron phase advance per turn is small. Therefore, the global synchrotron 
tune is not a general criterion in the classification of whether the adiabatic approx- 
imation fails. The evolution of 6, is more relevant. The argument is based on the 
fact that R56 can be either positive or negative, therefore the local synchrotron phase 
advance can be either positive or negative. While in the transverse dimension, the 
drift length and betatron phase advance are always positive. 

The breakdown of adiabatic approximation may have crucial impacts on the study 
of both the single-particle and collective effects. For linear single-particle dynamics, 
the longitudinal and transverse dimensions should be treated the same way on equal 
footing and SLIM formalism can be invoked. The treatment of nonlinear single- 
particle dynamics is more subtle as the longitudinal dynamics now is strongly chaotic. 
For the collective effects, many classical treatments should be re-evaluated and some 
new formalism needs to be developed. For example, the Haissinski equation [32] for 
calculating the equilibrium beam distortion cannot be applied directly then. Also, 
to our knowledge, there is no discussion on coherent synchrotron radiation (CSR)- 
induced microwave instability in a longitudinal strong focusing ring. The scaling law 
obtained in the longitudinal weak focusing [33] cannot be applied directly. 3D CSR 
effects and also the impact of bunch lengthening from transverse emittance on CSR 
needs more in-depth study. This is especially true for an SSMB ring, considering the 
fact that the beam width there is much larger than the microbunch length, while the 
contrary is true in a conventional ring. The contribution from horizontal emittance 
can easily dominate the bunch length at many places in an SSMB ring. This on the 
other hand, will be helpful to suppress unwanted CSR and may also be helpful in 
mitigating intrabeam scattering (IBS) [34, 35], as extreme short bunches occur only 
at limited locations. The IBS in a longitudinal strong focusing ring, and a general 
coupled lattice, also deserves special attention. To our knowledge, the IBS formalism 
of presented in Refs. [36, 37] can be applied for such purposes, as they are based 
on 6 x 6 general transport matrices. An IBS formalism can also be developed based 
on SLIM formalism [9], in which eigen analysis has been invoked and applies to 3D 
general coupled lattice with longitudinal strong focusing. 


2.1.7 Thick-Lens Maps of a Laser Modulator 


In the previous discussions, we have approximated the function of a laser modulator 
by a thin-lens RF-like kick. This means that we have ignored the phase slippage or 
R56 of the laser modulator itself. We need to know if this approximation is valid or 
under what circumstance we can use this approximation. 

Here we first derive the phase slippage factor of the undulator and then get the 
thick-lens transfer matrix of a laser modulator. The path length of an electron with a 
relative energy deviation of ô wiggling in a planar undulator is 
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Lu Lu 1/K 2 
L(6) = f V14+(4'dz 7% f 1+ 5 (= cos(kuz)) dz 
0 0 Y 


ee ee ee 
E ARa] 2y J 


in which k, = 27/2, is the undulator wavenumber, y, is the Lorentz factor corre- 
sponding to the resonant energy. The Rs¢ of an undulator is then 


(2.85) 


Lu ~ L(6) J Ly — L, F K?/2) 


R56 = = 
5 y2 y2 


= 2N,A0, (2.86) 


where N, is the number of undulator periods, Ay = EY : A, is the central wavelength 


of the on-axis fundamental spontaneous radiation. As can be seen from Eq. (2.86), 
the undulator R56 is twice the slippage length of the undulator radiation. 

As mentioned, the RF or laser modulator kick in linear approximation is like a 
longitudinal quadrupole and the R56 of the laser modulator is like the longitudinal 
drift space length of this longitudinal quadrupole. Assuming that the energy modu- 
lation is uniform along the undulator, then similar to the thick-lens quadrupole in the 
transverse dimension, we have the thick-lens transfer matrix of a laser modulator 


cos (v Rsoh) Sanh 


sin(./— h į if Rsh < 0, 
ED cos (V=Rssh) 
M = cosh (/Rs6h) Rss sinh( V Rs6h ) (2.87) 
hsi NED © Reh |, if Rsh > 0 
h sinl zeh ; : 
ar COSH (VRs6h) 


A laser modulator in linear approximation is therefore like a thick-lens quadrupole 
in the longitudinal dimension. A thin-lens approximation is applicable when | Rs6/| 
<i. 

After discussing the linear map, now we take into account the fact that the modula- 
tion waveform of a laser is actually sinusoidal. In principle, we can get an approximate 
analytical nonlinear thick-lens transfer map of the laser modulator using the tech- 
niques of drift and kick and Lie algebra [38, 39], by slicing the interaction into several 
smaller pieces and concatenate the maps of thin-lens kickes and drift spaces. Here 
we use a more straightforward method, i.e., to implement a symplectic kick map as 
below in a numerical code, to give the readers a picture. The kick map implemented 
in the code is as follows 


fori =1:1:N, 
z =Z + àoô 
ô = ô + A; sin (kzz) (2.88) 
z = Z + àoô 


end 
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Fig. 2.7 Impact of modulator undulator period number N, (thus R56) on the single-pass modulation 
process in a laser modulator. The injected beam is “ô = 0 and z € [0, Ag]”. (x, y) in the figures 
from ee right: (Zentrances Sexit)} (Zentrance, Zexit)3 (Zexit, Sexit). Parameters used: Ay, = 1 um, A= 
1x 107 


In other words, we have split the undulator into N, small “i dispersion + modulation 
+ 5 dispersion”. For a plane wave A; = x with A being the total energy modulation 
strength, while for a Gaussian laser beam A; is a function of i. 

To simplify the analysis, in the example numerical simulation, we consider the case 
of a plane wave, i.e., Aj = y First, we want to see how the single-pass modulation 
waveform is like based on Eq. (2.88). We choose parameters àz = Ap = l um, A = 
1 x 107°. The single-pass modulation of a line beam, with “ô = 0 and z € [0, Ao]”, 
as a function of N,,, namely Rs6, of the undulator is shown in Fig. 2.7. When changing 
N,,, we keep A unchanged. As can be seen, the waveform deviates from sine wave 
when N, increases. The beam distribution in phase space at the undulator exit (right 
sub-figure of Fig. 2.7) is similar to the beam de-coherence in an RF bucket. 

Now we consider the multi-pass cases, i.e., we consider the impact of modulator 
Rs56 on the phase space bucket. But here we do simulation only for the longitudinal 
weak focusing with a single laser modulator, as we only aim to give the readers 
a picture about such impact. We use parameters of A, = ào = 1 pm, A=1x 
107-3, Co = 100m, n = 5 x 107” inthe simulation, and choose to observe the beam 
opposite the modulator center where a, = 0 with N, = 0, 40, 320, respectively. Note 
that 7 is the phase slippage factor of the whole ring, including the modulator. When 
changing the undulator priod number N,, we keep n a constant. The results are shown 
in Fig. 2.8. It can be seen that the modulator R56 only distorts the bucket slightly 
when N, is 40. But when N, is as large as 320, it will have a profound effect on 
the longitudinal phase space topology. Its impact on longitudinal strong focusing is 
more subtle as the particle motion in a longitudinal strong focusing ring is strongly 
chaotic if the nonlinear modulation waveform is taken into account. The study of 
such effect can refer more straightforwardly to numerical simulations. Besides, the 
undulator R56 could also have an impact on the coherent radiation induced collective 
instability in the laser modulator [40, 41]. 
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Fig. 2.8 Impact of modulator undulator period number N, (thus R56) on the longitudinal 
phase space bucket in longitudinal weak focusing regime. Parameters used: Àz = 1 pm, A= 
1x 10-3, Co = 100m, n =5 x 1077 


2.2 Nonlinear Longitudinal Dynamics 


After resolving the issue of quantum diffusion of z by means of dedicated lattice 
design, we can then apply the quasi-isochronous or low-phase slippage method to 
realize short bunches in SSMB. However, the phase slippage is actually a function 
of the particle energy 

n(6) = no +md+ m5? +... (2.89) 


When no is sufficiently small, the higher-order terms in Eq. (2.89) may become 
relevant or even dominant, and the beam dynamics can be significantly different from 
those in a linear-phase slippage state. Proper application of dedicated sextupoles and 
octupoles may be needed to control these higher-order terms. 

The beam dynamics of the quasi-isochronous rings have been studied by many 
authors [15, 16, 42]. Here, we wish to emphasize two points that have not been 
well investigated before and might be important, for example, in the SSMB proof- 
of-principle experiment to be introduced in Chap. 5 and the longitudinal dynamic 
aperture optimization in SSMB. 


2.2.1 For High-Harmonic Bunching 


For seeding techniques such as coherent harmonic generation (CHG) [43] and high- 
gain harmonic generation (HGHG) [44, 45], it seems that to date, linear phase slip- 
page or Rs has been applied for microbunching formation. Here, we wish to point 
out that one can actually take advantage of the nonlinearity of the phase slippage for 
high harmonic generation. Intuitively, this is because a sinusoidal energy modulation 
followed by a nonlinear phase slippage can lead to a distorted current distribution, 
which, in some cases, can lead to large bunching at a specific harmonic number. 
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Fig. 2.9 The beam evolution in longitudinal phase space, final current distribution and bunching 
factor, when 7(5) = 716 is used for microbunching in CHG or HGHG, as modeled by Eq. (2.90) 


Figure 2.9 shows an example simulation of using 7(6) = 716 for microbunching. It 
can be seen that there is significant bunching in the second and fourth harmonics, 
while no bunching is produced in the fundamental and third harmonics. The reason 
can be found from the following derivation of the bunching factor. 

The microbunching process in the case of a single energy modulation followed 
by a dispersion section, as that in CHG and HGHG, can be modeled as 


8’ = 5 + Asin(kz2), 
z’ = z—n(8)Cod’, (2.90) 


where ky, = 27/4, is the wavenumber of the modulation laser, A is the elec- 
tron energy modulation strength induced by the laser. The bunching factor at the 
wavenumber of k is defined as 


b(k) = | de oi. (2.91) 


[oe] 


where p(z) is the normalized longitudinal density distribution of the electron beam 
satisfying ae dzp(z) = 1. According to Liouville’s theorem, we have dzdé = 
dz'dé'. Therefore, the bunching factor can be calculated in accordance with the 
initial distribution of the particles oo(z, ô) as 


[o6] [e6] es 
b(k) = i J dzd6 polz, de &), (2.92) 
—oo J—00 
Here we consider the simple case of n(6) = no + 716, then 


o. poo ik( noCoô+nı CoS? U0 
ba) = f J dzd8 po(z, ô)e (measmcosr+ ag) 
-0 J -o0 (2.93) 


Co A2 
2 


ik| -+0004 +2n: CoôA) sin(k,z)— 7 cos(2k1)] 
e y 


Adopting the notation Y = k (noCoA + 2nıCoôA), = —kn,CoA?/2, and using 


the mathematical identity e'* %9 = °°. e” J, [x], we have 


42 2 SSMB Longitudinal Dynamics 


ioe) foe) 
e—ikz+iY sin(kyz)+iZ cos(2kyz) — > 5 Ipl¥ gl Zlexe ( iftk (p — 2q)k1)z q5]): 


i i (2.94) 
If the initial beam is much longer than the laser wavelength, and considering that 
z ie e` dow = 8(t), where ô(t) is the Dirac delta function, the bunching factor 
will not vanish only if k = (p — 2q)k,. The bunching factor at the n-th harmonic of 
the modulation laser is then 


a ink, (| noCod+m1C 524 CoA? oo 
n=) dê e (w oô+n1 Co 2 ) 8) > Ins 2mlY | Jin[Z]. (2.95) 


02. m=—0o 


Here we consider the simple case of an initial Gaussian energy distribution p9(6) = 
JEP (- +), where o; is the initial RMS energy spread. If nı = 0, then Y = 


nkLnoCoA, Z = 0, and paws Jny2ml Y |Jin(Z] = J [nkt noCoA], and we have 


m=— 00 


B otwor] l 2.96) 


Da = Jn [nkrnoCoAlexp | 2 
which is a familiar result for HGHG [44] if we adopt the notation Rs6 = —noCo. If 
No = 0, then Y = 2nk,n,CodA, Z = —nkyn, CoA? /2, meaning that we have 


1 ce 1CoA? er\e 
bn = Tal dô exp [int (nco t 1 2 ) exp 202 D Jn+2m[Y]Jm[Z]. 


(2.97) 
The two exponential terms in the integral are even functions of ô, while 
Jn+2mlY¥ ]Jin[Z] is an odd function of 6 when n is odd; thus, b, is nonzero only 
for an even n. This is why bunching occurs only in the second and fourth harmon- 
ics but not in the fundamental and third harmonics when we use n(6) = 716 for 
microbunching, as shown in Fig. 2.9. 
Following the derivations and according to the relation 


n 


cos" (x) = 5 y (3) cos(2m — n)x, (2.98) 


m=(n+1)/2 


it can be seen that the energy modulation at the fundamental frequency can be cast into 
[i x (n — 2p) + j x (n — 2q)]-th harmonic bunching through the term n,,_ ;6” in the 
function of n(6). For an odd n, bunching at all harmonic numbers are possible, while 
for an even n, only bunching at the even harmonic numbers is possible. The optimal 
bunching condition for a specific harmonic requires the matching of (6) with the 
energy modulation strength. However, the analytical formula for the bunching factor 
will become increasingly involved with more higher-order terms of the phase slippage 
considered. Thus, it would be better to refer to numerical code to calculate and 
optimize the bunching factor directly for a specific application case. For storage rings, 
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another relevant point is that the distribution of the particle energy in the nonlinear 
phase slippage state may also have an impact on the high harmonic generation, and 
this phenomenon is also easier to be studied by means of numerical simulation. 

The approach of applying a nonlinear phase slippage for high harmonic bunching 
can be considered to share some similarity with echo-enabled harmonic generation 
(EEHG) [46, 47], in which the sinusoidal energy modulation and dispersion in the 
first stage can be viewed as the source of the distorted current distribution in the 
second stage of modulation and dispersion for microbunching. We have also noticed 
the work on optimizing the nonlinearity of the dispersion to increase the bunching 
factor for EEHG [48]. Based on similar considerations, tricks can also be applied on 
the energy-modulation waveform using different harmonics of the modulation laser, 
for example, forming a sawtooth waveform to boost bunching, as will be discussed 
in Chap. 3. 


2.2.2 For Longitudinal Dynamic Aperture 


Similar to the transverse dimension, there is a region in the longitudinal phase space 
outside of which particle motion is not bounded and can be lost in a ring. We refer 
this stable region as the longitudinal dynamic aperture. Here in this section, we want 
to show that, by properly tailoring the nonlinear phase slippage, the longitudinal 
dynamic aperture can be enlarged significantly compared to the case of a pure linear 
phase slippage. Only symplectic dynamics is considered in this discussion. 


2.2.2.1 Longitudinal Weak Focusing 


The longitudinal dynamics of a particle in a ring with a single RF can be modeled 
by the kick map 
Ôn+1 = Ôn F Al[sin(krrZn) — sin ds], 


(2.99) 
Zn+1 = Zn — n(ên+1)Coôn+1, 


where A sin ¢; = Uo/ Eo where Us is the radiation loss of a particle per turn. For the 
case of longitudinal weak focusing, the kicik map can be approximated by differenti- 
ation and Hamiltonian formalism can be invoked for the analysis. Denote ¢ = kez, 
then the equation of motion is 


aa (2.100) 


dọ _ kren n+ )Co 5 _ dH 
= T = 
A . . 
T7 neno — sin ġ;) = =p 


with Tọ being the revolution period of the particle in the ring. For n (8) = no + 16 + 
726°, the corresponding Hamiltonian is 
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Hd, 8) = onr ( 3m8? + Fm3? + znd!) + £ leos — cosg, + (Ø — gy sind. 
3 0 
(2.101) 

In writing down the closed-form Hamiltonian, we have implicitly assumed that the 
motion is integrable, i.e., there is no chaos. But we need to keep in mind that the 
dynamics dictated by Eq. (2.99) is actually chaotic even with a linear phase slip- 
page [13]. But here we ignore this subtle point as the chaotic layer is very thin in the 
longitudinal weak focusing regime. We remind the readers that the chaotic dynamics, 
for example the bucket bifurcation, can actually also be applied for ultrashort bunch 
generation [49]. 

To analyze the motion, we need to find the fixed points of the system 


| 2H = : | sing, — sing = 0, (2.102) 


5n(5) = 0. 


To determine whether a fixed point is stable or not, we need to check the trace of 
the Jacobian matrix around the fixed point. If n(ô) = no, there are two sets of fixed 
points: 


fe 2 (Øs, 0), (2.103) 


UFP: (x — ġ;, 0), 


in which SFP stands for stable fixed point while UFP stands for unstable fixed point. 
If n(6) = no + 716, there are four sets of fixed points: 


SFP: (dy, 0), (z - $,-2), 


nı 


UFP: (x — ¢,,0), ($ -2). (2.104) 


If n(8) = no + 16 + 26’, there are six sets of fixed points 


2 2 
se: 0, (1-2 + (at) B). (e -0 (a) 8). 


2 2 
UFP: (r — On: 0), (o m (#4) 7 , (+ m (#4) 3 
(2.105) 


To see the impact of 7; and 72 on the longitudinal phase space bucket, some 
numerical simulations are conducted. We choose the observation point at the middle 
of the RF, where w, = 0 and the beam distribution in the longitudinal phase space is 
upright. The results of the impact of nı and m on longitudinal dynamical aperture 
are shown in Figs. 2.10 and 2.11, respectively. Note that in the plots, we have used 
the longitudinal coordinate z rather than the phase @¢. 

As we can see in Fig. 2.10, the emergence of 7; will make the bucket asymmetric 
in 6, which is as expected as the n = no + 716 is asymmetric in ô. In both directions 
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Fig. 2.10 The impact of 7; on the longitudinal phase space bucket in the longitudinal weak focusing 
regime. Simulation parameters: Arf = | ym, A = 1 x 107-3, Co = 100 m, no = 5 x 1077 
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Fig.2.11 The impact of 72 on the longitudinal phase space bucket in the longitudinal weak focusing 
regime. Simulation parameters: Arf = | ym, A = | x 1073, Co = 100 m, no = 5 x 1077 


(positive or negative), the bucket size shrinks with the increase of 7; and the bucket 
becomes like an upright a-shape, so they are usually referred to as a-buckets. Note 
that we can also classify the bucket to be an RF-bucket or an a-bucket according to 
whether 6 = 0 or 7(6) = 0 at the bucket center, respectively. Such a classification is 
more reasonable from beam dynamics consideration. a-bucket is also a method to 
generate short bunch and there are many interesting beam dynamics issues related 
to such buckets [16]. 

However, the impact of m is different. As can be seen from the simulation results 
presented in Fig. 2.11, the bucket is still symmetric in ô as expected. Besides, when 
a < 0, the stable region of the bucket can be even larger than the case without n2. We 
will see later that in the case of longitudinal strong focusing, such observation can be 
even more notable. Therefore, we can tailor the phase slippage factor as a function 
of energy in the case of longitudinal strong focusing to enlarge the longitudinal 
dynamic aperture. This is very helpful as usually the longitudinal dynamic aperture 
in a longitudinal strong focusing ring is not a trivial issue and needs to be optimized 
to guarantee a sufficient quantum lifetime for example. 
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2.2.2.2 Longitudinal Strong Focusing 


For a longitudinal strong focusing ring, the particle motion is strongly chaotic and 
not integrable, and we cannot get a closed-form Halmitonian for analysis anymore. 
Therefore, we use numerical simulations to study the dynamics directly. Instead of 
a comprehensive investigations, here in this section we aim to give some qualitative 
remarks on the role of the nonlinear phase slippage, i.e., a proper tailoring of the 
nonlinear phase slippage can enlarge the longitudinal dynamic aperture significantly. 

Let us use the schematic layout shown in Fig. 2.6 and parameters choices given 
in Table 2.2 as an example for illustration. Note that the modulation wavelength 
used here is Arp = | pm. We choose to observe the beam at the radiator center, and 


suppose the ring is symmetric with respect to the radiator. The one-turn kick map is 
then 


z=z+ RS, 

6=6 + h/ kgf sin(krfz), 

zzi R26, (2.106) 
6=6 + h/ kgf sin(krfz), 


As we aim to present the main physical picture, here we only consider the nonlin- 
earity of the main ring first, i.e., RO (ô) = —Co (no +ô + mô’). RY and RO in 
principle can also be a function of 5. The simulation results are shown in Figs. 2.12 
and 2.13. 

From Fig. 2.12, we know that, like that in the weak focusing case, nı makes the 
bucket asymmetric in ô and shrinks the bucket size whether 7; is positive or negative. 
From Fig. 2.13, we can see that when a < 0, a proper 72 can help to merge the island 
buckets with the main bucket and broaden the stable region of the phase space, i.e., 
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Fig. 2.12 The impact of 7; on the longitudinal phase space bucket in the longitudinal strong 
focusing regime. Simulation parameters: Arp = 1 pm, h = —50000 m!, RP = 15 um, Co = 
100 m, no = 1 x 1076 
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Fig. 2.13 The impact of 72 on the longitudinal phase space bucket in the longitudinal strong 


focusing regime. Simulation parameters: ARF = 1 pm, h = —50000 m!, RO = 15 um, Co = 
100 m, np = 1 x 1076 


the longitudinal dynamic aperture, significantly. A proper n2 makes the amplitude 
dependent tune shift favorable for the motion to be stable. Note that the fixed points 
of the island buckets may not have period-1 but period-n stability. 
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Chapter 3 A) 
SSMB Transverse-Longitudinal get 
Coupling Dynamics 


After dedicated efforts devoted to the longitudinal dynamics to realize an ultrasmall 
longitudinal emittance and ultrashort bunch length for coherent radiation generation, 
we need to make sure that the coupling arising from transverse dynamics does not 
degrade or even destroy the longitudinal fine structures. Such an argument is based 
on the observation that the transverse beam size in an SSMB ring can be orders 
of magnitude larger than the desired microbunch length. This is the basic motiva- 
tion for us to investigate the transverse-longitudinal coupling (TLC) dynamics. In 
this chapter, we start from the linear TLC and then investigate the nonlinear TLC 
dynamics. For the linear dynamics, first we analyze the passive bunch lengthening 
induced by bending magnets. We then emphasize the fact that TLC can actually be 
actively applied for efficient bunch compression and high harmonic generation when 
the transverse emittance is small. We present three theorems on the application of 
such TLC schemes, with their implications discussed. Further, we have analyzed 
the contribution of modulators to the vertical emittance from quantum excitation, 
to obtain a self-consistent evaluation of the required modulation laser power when 
applying these coupling schemes in a storage ring. The theorems and related anal- 
ysis provide the theoretical basis for the application of TLC in SSMB to lower the 
requirement on the modulation laser power, by taking advantage of the fact that the 
vertical emittance in a planar ring is rather small. Based on the investigations, we 
have presented example parameters sets for the envisioned SSMB storage ring to 
generate high-power EUV and soft X-ray radiation at the end of this dissertation. 
The relation between our TLC analysis and the transverse-longitudinal emittance 
exchange is also briefly discussed. For the nonlinear dynamics, we present the anal- 
ysis and the first experiment proof of a second-order TLC effect on the equilibrium 
beam parameters, which can help to improve the stable beam current and coherent 
radiation power of a ring working in quasi-isochronous regime. Parts of the work 
presented in this chapter have been published in Refs. [1—4]. 
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3.1 Linear Transverse-Longitudinal Coupling Dynamics 


3.1.1 Passive Bunch Lengthening 


In a linear transport line without bending magnets, the transverse and longitudinal 
motions are decoupled in a first-order approximation. However, the situation changes 
when there are bending magnets. Particles with different horizontal (vertical) posi- 
tions and angles will pass through the horizontal (vertical) bending magnets along 
different paths, resulting in differences in the longitudinal coordinate. The transverse 
motion can thus be coupled to the longitudinal dimension. When traversing the bend- 
ing magnets, particles with different energies will also pass along different paths and 
exit with different horizontal (vertical) positions and angles. The longitudinal motion 
can thus also be coupled to the transverse dimension. The physical pictures of the 
linear TLC introduced by the bending magnets are shown in Fig. 3.1. Although this 
passive TLC is a well-understood effect [5-9], here we present a concise analysis 
of this effect with an emphasis on its vital role in microbunching formation and 
transportation for both the transient and steady-state cases. 

We start with a planar x-y uncoupled lattice and assume that the RF cavities are 
placed at dispersion-free locations. We temporarily ignoring the vertical dimension, 
and use the state vector X = (x, x’, z, 5)’. The subscripts 5,6 are used for z, ô 
for consistency with literature. Hereafter, the subscript , in this section is omitted 
unless necessary. As introduced in Sect. 2.1.1, the betatron coordinate, defined by 
Xg = BX, is first used to parametrize the transport matrix in a diagonal form. The 
transport matrix of Xg from sı to sz is then 


Mx (82, 51) 0 
= $ 
Mga(s2, 51) = ( 0 Mz, (52, 81)) ° (3.1) 
Following Courant and Snyder [10], we write Mx, (s2, $1) as 
Mx (92, 91) = AT! (s2)T (s2, s1)A(s1), (3.2) 


x x x 


oL perk , ek, 
< ‘ S 


(b) (c) 


Fig. 3.1 Linear transverse-longitudinal coupling induced by a bending magnet. Particles with 
different horizontal positions (a) and angles (b) pass the horizontal bending magnet along different 
paths, resulting in longitudinal coordinate differences. Particles with different energies (c) also 
pass the horizontal bending magnet along different paths, resulting in horizontal position and angle 
differences 
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with i 
A(si) = Ea ) (3.3) 
TBs VPG) 
and 
_ ( cosy sinyin 
T6251) = (< tk | ; (3.4) 
where 


S2 1 
pi = p — yı = J Toa (3.5) 


51 


is the betatron phase advance from sı to s2. For Mz, (s2, s1), the expression is similar 
to Mx, (s2, 51), but note that if we want to calculate the synchrotron phase advance 
similar to Eq. (3.5), the distance s should be replaced by the effective longitudinal 
drift space, i.e., F = —7(s2, s1)Co defined in Eq. (2.50). If there is no RF cavity 
between sı and s2, we have 


M.,, (52, 81) = (( G n i (3.6) 


The transport matrix of X from sı to sz is then 
M(5p, s1) = Bo! (s2)Mg (s2, s1)B (s1). (3.7) 
After some straightforward algebra, M(s2, s1) can be expressed as 


Ry Riz 0 Dz — Ry Dy — R2 D; 
Ro, R22 0 D} — Ro Dy — R2 D; 


M(52, sı) = Rs, R52 1 F — Rs, Dy — Rs2D} 
0 00 1 
B2 . 
Ri = g% Yiz + &ı sin Yaz], 
Ri =:4/ Bi Bo sin Yiz, (3.8) 


Ry} = — [C1 + aa) sin yiz — (a1 — a2) cos Y2], 


1 
~ Bi Bo 


Ry = J Po Wi2 — œ sin Y2], 
b2 


Rs; = Ry D2 — RD, + D}, 
Rs2 = —Ri2D) + R2 Də — Dj. 
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This matrix can then be used to analyze both the transient and steady-state cases of 
TLC. Note that for a given lattice, F is a function of the initial dispersion and dis- 
persion angle (D,, D{) at sı, although the transfer map M(s2, s1) for the state vector 
X is not, as the transport matrix is fixed once the lattice is given. This dependence is 
a result of the matrix parametrization. 

We consider first the influence of the betatron oscillation on the longitudinal coor- 
dinate. With the help of the Courant-Snyder parametrization, the betatron oscillation 
position and angle at the starting point sı can be expressed as 


xı = y 2J ßı cos Wy, 
xi = —y/2J / Bı (œ1 cos Yı + sin y1), (3.9) 


where J = A (yx? + 2axx’ + px") is the betatron invariant or action of the particle. 
The longitudinal coordinate displacement relative to the ideal particle due to the 
betatron oscillation from sı to s2 is then 


Az = R51x1 + Rs2x| = y 2JHı sin(yı — x1) — V2JA% sin(y72 — x2), (3.10) 


where to obtain the final concise result, D and D’ have been expressed in terms of 
the chromatic H{-function and the chromatic phase x, defined as 


D = YHB cos x, 
D' = —J/H/B (æ cos x + sin x), (3.11) 


where H = y D? + 2aDD! + BD”. If there is no dipole kick between point 1 and 
point 2, H stays constant and x2 — x; = Y2 — Yı, which means Az = 0. Physically, 
this means the contribution of transverse emittance to the bunch length does not 
change during drifting or experiencing quadrupole kicks, as these manipulations 
only affect the beam distribution in the transverse phase space. This argument can 
also be clearly observed in Fig. 2.3c. 

From Eq. (3.10), the root-mean-square (RMS) value of the transient bunch length- 
ening of a longitudinal slice from sı to s2 caused by this linear TLC can be calculated 
to be 


Oaz = Je [m +H — 2V HIH: cos (Ayri — Ax) (3.12) 


The RMS bunch lengthening of an electron beam longitudinal slice after n complete 
revolutions in a ring, due to betatron oscillation, is then 


Oaz = 2y EH |sin(n7v,)| . (3.13) 


The above equations can be used to explain the dependence of the coherent 
synchrotron radiation (CSR) repetition rate on the betatron tune in the bunch slic- 
ing experiment reported in Ref. [8]. A similar approach can be applied to analyze 
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microbunching preservation with beam deflection, for example in FEL multiplexing. 
These equations are also useful for evaluating the influence of the coupling effect 
in the SSMB proof-of-principle experiment [11, 12], which is to be presented in 
Chap. 5. 

If the particle starts with a relative energy deviation of ô, then 


Az = Rs1x1 + Rs2x, + (F — R51 Dı — Rs2D})6 
= y 2JH; sin(Yı — x1) — V2J AH sin(W2 — x2) + Fô. (3.14) 
Note that in Eq. (3.14), the betatron invariant and phase should be calculated accord- 
ing to 
xp =x = D8 = y 2J cosy, 
x, =x’ — D'8 = —/2J/B(acos y + sin Y), 


1 2 1 12 
J= (vx + 2axgx, + Bxe). (3.15) 

For a periodic system, if we observe the particle at the same place n periods later, 
then 


Az = V2JH [sin(w — x) — sin( + 2nz vy — x)] — nnCoô, (3.16) 


where 27 v, is the horizontal betatron phase advance per period. 

We can also obtain the equilibrium second moments in a storage ring by following 
the Courant-Snyder parametrization one step further. The result is the same with 
Eq. (2.18) obtained by SLIM. As can be seen from Eqs. (2.18) and (2.19), if there 
is only passive TLC introduced by bending magnet, the transverse emittance can 
lengthen the bunch at places where H, 4 0, 


(3.17) 


Similarly the energy spread can broaden the beam width at places where D Æ 0, 


Ox = VerBr + €y: D? = y €x bx + o? D2. (3.18) 


To give the readers a more concrete feeling about the bunch lengthening from 
this passive TLC, we have presented in Fig. 2.3c some calculations based on the 
MLS lattice. As can be seen, indeed that the coupling from horizontal emittance can 
contribute significantly, or even dominant the bunch length at places where Hy is 
large. This observation will especially be true in an SSMB ring, where the trans- 
verse size is much larger than the microbunch length. Therefore, the dispersion and 
dispersion angle should be controlled in precision at places where ultrashort bunch 
is desired, for example at the radiator. The bunch lengthening from the transverse 
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Fig. 3.2 Beam current 10 

distributions at places with —H, =0 
different Hy. Bunch length —H: #0 
in an SSMB ring can easily F 1 


be dominated by the 
transverse emittance in 
places where Hy 4 0 


Intensity 


z/ÀL 


emittance will make the current distribution in an SSMB storage ring less sharp and 
more like a coasting beam as places where H, 4 0, as shown in Fig. 3.2. Here we 
make a remark that this coupling effect may be helpful for suppressing unwanted 
CSR and may mitigate the intrabeam scattering (IBS) in SSMB or other applications, 
as extremely short bunches emerge only at dispersion-free locations. 


3.1.2 Coupling for Harmonic Generation and Bunch 
Compression 


The analysis in the above section may lead us to conclude that TLC always lengthens 
the bunch and degrades the microbunching. This, however, is not true, as the above 
analysis is based on the assumption of a planar x-y uncoupled lattice with only 
the passive coupling induced by the bending magnets. In addition to this passive 
coupling, an RF cavity (laser modulator in SSMB) placed at a dispersive location, a 
transverse deflecting RF cavity, etc., are other sources of coupling that can be used 
for subtle manipulation of particle beam in 6D phase space. In fact, one can take 
advantage of TLC for efficient harmonic generation or bunch compression when 
the transverse emittance is small. The reason is that there is some flexibility in 
tailoring the projection of the three eigen emittances of a beam into different physical 
dimensions, although their values cannot be changed in a linear symplectic lattice. 
Here in this and the following sections we investigate the active application of TLC 
for harmonic generation and bunch compression by taking advantage of the fact that 
the vertical emittance of an electron beam in a planar x-y uncoupled ring is rather 
small. 
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Entrance Planar lattice Modulator Planar lattice Radiator 


d: depends solely on ¢, I >> EE ap GS < depends solely on <, 


Fig. 3.3 A schematic layout of applying y-z coupling for bunch compression 


3.1.2.1 Problem Definition 


Let us first define the problem we are trying to solve. We use y-z coupling as an 
example for the analysis, since we aim to exploit the small vertical emittance. The 
case of x-z coupling is similar. Suppose the beam at the entrance of the bunch 
compression section is y-z decoupled, i.e., its second moments matrix is 


€y By; TEyQAyj 0 0 


o [| Hey Oye €EyYyi 0 0 
aa 0 0 Ezbzi —EzQzi |’ Gaa 
0 OEM Ec Yzi 


where «œ, 6 and y are the Courant-Snyder functions [10], the subscript ; means initial, 
and €, and €, are the eigen emittances of the beam corresponding to the vertical and 
longitudinal mode, respectively. For the application of TLC for bunch compression, 
it means that the final bunch length at the radiator depends only on the vertical 
emittance €, and not on the longitudinal one €,. The magnet lattices are all planar 
and x-y decoupled. 

The schematic layout of a TLC-based bunch compression section is shown in 
Fig. 3.3. We divide such a section into three parts, with their transfer matrices given 
by 


r33 r340 d R33 R340 D 
r43 r44 0 d' : ; Ras Rag 0 D' 
M; = , M: = Modulation kick map, M3 = 
i r53 rs4 l 156 : paoi Rs3 Rs4 1 Rs6 
0 001 0 001 
r53 = rad — r33d', rs4 = —r34d' + rad, r33r44 — rara3 = 1, 


Rs3 = Ra3D — R33D', Rs4 = —R34D' + R4D, R33R44 — R34Ra3 = 1, 

(3.20) 
with M; representing “from entrance to modulator”, M3 representing “modulation 
kick” and M; representing “modulator to radiator”. Note that M; and M; are in 
their general thick-lens form as analyzed in last section. The transfer matrix from 
the entrance to the radiator is then 


T = M;M-M.;. (3.21) 


From the problem definition, for o; (Rad) to be independent of €,, we need 
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T55 = 0, 


3.22 
Tss = 0. Gan 


3.1.2.2 Three Theorems on Transverse-Longitudinal Coupling 


Given the above problem definition, we have three theorems which dictate the relation 
between the modulator kick strength with the optical functions at the modulator and 
radiator, respectively. 


Theorem one: If 
1000 
0100 
M, = 0010]? (3.23) 


00h1 


which corresponds to the case of a normal RF or a TEMO0 mode laser modulator, 
then 


h? (Mod)H, (Mod)H, (Rad) > 1. (3.24) 
Theorem two: If 
1000 
0110 
M, = 00101? (3.25) 
tool 


which corresponds to the case of a transverse deflecting RF or a TEMO1 mode laser 
modulator or other schemes for angular modulation, then 


t’ (Mod) £, (Mod)H, (Rad) > 1. (3.26) 
Theorem three: If 
1 0 kO 
0 1 00 
M, = 0010)? (3.27) 
0-k01 


whose physical correspondence is not as straightforward as the previous two cases, 
then 
k? (Mod) yy (Mod)H, (Rad) > 1. (3.28) 


3.1.2.3 Proof 


Here we present the details for the proof of Theorem one. The proof of the other two 
is just similar. From the problem definition, for o; (Rad) to be independent of €,, we 
need 
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Ts5 =hR5 + 1 = 0, 


/ (3.29) 
Ts6 = dRs3 +d R54 + rse(h Rs6 + 1) + R56 =0. 
Note that the harmonic generation schemes in FEL like phase-merging enhanced har- 
monic generation (PEHG) [13, 14] and angular dispersion-induced microbunching 
(ADM) [15] can be viewed as specific examples of the above general relations [2]. 
Under the conditions of Eq. (3.29), we have 


AB 
T= e |) 3 (3.30) 
with 
A= r33 R33 + ra3 R34 + rs3hD r34R33 + r44 R34 + rs4h D 
r33 Ras + ras Ra4 + rs3h D' r34Ra3 + r44 R44 + rs4h D' 
B= hD dR33 + d'R34 + (rssh + 1) D 
~ \AD! dRy + d' Rag + (rsh + 1) D' (3.31) 
C= r33 R53 + ras Rs4 r34Rs3 + r44 R54 l 
ii rs4h 
g Gee rssh + i) i 
The bunch length squared at the modulator and the radiator are 
2 2 
i153 — yirs +r 
ao? (Mod) =€, (Bzi E 2æzirs6 A Vait36) +4 E (By 53 54) 54 
yi 
= €,B,(Mod) + €H, (Mod), (3.32) 
2 2 
ji T: 3: ee: iT: + T. 
o?(Rad) = €, (ByiTs 7 s) + Toa e, H, (Rad). 
yi 


According to Cauchy-Schwarz inequality, we have 


2 2 
[(Byirs3 = ayirs4) + 34] [(Byi 75s — ayi Tsa) + 7, | 
Êyi Êyi 


h? (Mod) Hy (Mod) Hy (Rad) = h? 


r 
z ~ [~ (Byir5s3 — &yirs4) T54 + 154 (Byi T53 — ay; Ts4)]° 
y 
= (Ts3rs4h — Ts4r53h)? = (Ts3 T64 — T3463)" - 
(3.33) 
The equality holds when mil rss rst) = eaa) The symplecticity of T 


requires that TST’ = S, where S = @ 1) and J = oe 5) , so we have 
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T T T T 
(a + BJB’ AJC’ + BJE Jat (3.34) 


CJA’ + EJB’ CJC’ +EJE’ )~ \OJT 


0 O 
hrsgh +1 


As shown in Eq. (3.31), E = ( ) , then EJE” = ( a) Therefore, 


00 
CJC? =J, (3.35) 
which means C is also a symplectic matrix. So we have 
Ts3Te4 — Ts4T63 = det(C) = 1, (3.36) 


where det(C) means the determinant of C. The theorem is thus proven. 


3.1.2.4 Dragt’s Minimum Emittance Theorem 


Theorem one in Eq. (3.24) can also be expressed as 


Ey D Ey 
H, (Mod) yeH, (Rad) ozy (Mod)o; (Rad) ` 


(3.37) 


|n(Mod)| > 
Vé 


Note that in the above formula, o,,(Mod) means the bunch length at the modu- 
lator contributed from the vertical emittance €,. So given a fixed e€, and desired 
o,(Rad), a smaller h(Mod), i.e., a smaller RF gradient or modulation laser power 
(Piaser & |h(Mod)|?), means a larger H, (Mod), thus a longer 0, (Mod), is needed. 
As |h(Mod)|o,(Mod) quantifies the energy spread introduced by the modulation 
kick, we thus also have 

o,(Rad)o3(Rad) > €. (3.38) 


Similarly for Theorem two and three, we have 


€ 

t(Mod)| > x 
|t (Mod)| = s,s (Mod)a,(Rad) ' (3.39) 

and z 
|k(Mod)| > : (3.40) 


oy g (Mod)o; (Rad) ’ 


respectively, and also Eq. (3.38). Note that in the above formulas, the vertical beam 
size or divergence at the modulator contains only the betatron part, i.e., that from the 
vertical emittance €. 

Equation (3.38) is actually a manifestation of the classical uncertainty principle 
[16], which states that 
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2 
X11 22 > Ehin 
233 X44 


e? (3.41) 
X55 266 = € 


IV 


IV 


min’ 

2 

min’ 

in which €min is the minimum one among the three eigen emittances €z 77,7717. In 
our bunch compression case, we assume that €, is the smaller one compared to €z. 
Actually there is a stronger inequality compared to the classical uncertainty principle, 
i.e., the minimum emittance theorem [16], which states that the projected emittance 
cannot be smaller than the minimum one among the three eigen emittances, 


2 — 2 2 

Ex pro = XX — 2i Z Enin’ 

E o = D Ey — Dy = ei (3.42) 
y,pro — “433 144 34 Z Smin’ : 
2 Pn 2 2 

Ez, pro — Us5 266 — 256 Z Ehin 


3.1.3 Normal RF or TEM00 Mode Laser for Coupling 


Now we investigate in more detail about the application of TLC in SSMB for bunch 
compression, using a TEMOO mode laser modulator for the modulation kick. This 
belongs to the category of Theorem one. 


3.1.3.1 Physical Picture 


According to Theorem one, given a vertical emittance €, and modulation kick 
strength h, in principle we can realize as short o, (Rad) as we want by lengthening 
ozy (Mod). In other words, we can lower H, (Rad) by increasing H, (Mod). However, 
such great flexibility of a TLC coupling scheme is not obtained without sacrifice. For 
a premicrobunched beam, and considering that the modulation waveform is actually 
a nonlinear sinusoidal, a bunch lengthening at the modulator will result in bunching 
factor degradation at the radiator. Another key point is that the modulator itself will 
contribute to the vertical emittance through quantum excitation since it is placed at 
a place where Hy # 0. We will elaborate these points more in this section. 

To give the readers a better picture before going into the mathematical details, 
here we summarize in Fig. 3.4 the main information to be presented in this section: 
(1) Compared to bunch compression or harmonic generation scheme in longitudinal 
dimension alone like high-gain harmonic generation (HGHG), TLC schemes like 
PEHG or ADM can reduce the required energy chirp strength, to realize the same 
bunch length compression ratio or harmonic generation number, when the transverse 
emittance is small. (ii) This lowering of energy chirp strength is realized through the 
bunch lengthening from transverse emittance at the modulator, which can degrade 
the bunching factor at the radiator for a pre-microbunched beam due to the nonlinear 
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Fig. 3.4 Application of TLC for bunch compression and harmonic generation, using a TEMO0 
mode laser modulator. Parameters used in this example plot: A; = 1064 nm, ozi = 30 nm and 
ozf = 3 nm for the case of a pre-microbunched beam, os; = 3 x 10-4, Oyi = 2 ym, oy; = l prad. 
The figures show the beam distribution evolution in the longitudinal phase space. Depending on 
the specific lattice scheme, the different colors in the plot correspond to different particle vertical 
positions, angles, or combination of them. The modulation waveforms are shown in the figure as 
the red curves 


nature of the sine modulation. (iii) Addition of the RF or laser harmonics is an 
effective way to mitigate this bunching factor degradation by broadening the linear 
zone of the modulation waveform. 


3.1.3.2 Bunching Factor 


Now we derive the bunching factor degradation at the radiator due to the bunch 
lengthening at the modulator, using ADM as an example. PEHG has a similar result 
as we have proven in the last section the general theorem of bunch lengthening in 
this kind of TLC schemes. Thin-lens kick maps in the last section are again used for 
the analysis, but now with the fact that the modulation waveform is a nonlinear sine 
taken into account. 

Putting in the optimized bunch compression conditions for ADM, namely h Rs6 + 
1 = Oand—d'D + Rs6 = 0, and using the mathematical identity e49" =>” a 
e!” J,,[a], the final bunching factor at the n-th laser harmonic in ADM is 
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2° œ fœ cad inku| -2 ee dz, 
b, = > Jin w f f f dyidy;dzie if TAAAC r 
—oo J—00 J—00 


m=—Oo 
(3.43) 
(a oe) will be non-zero only if m = n, where 


m 
n 


For a coasting beam, (e7!”*LC- 


the bracket (- - -) means the average over all the particles. Therefore, 
bn ,coasting = Jn (n)exp [— (nk, o,(Rad))* /2] . (3.44) 


Note that o; (Rad) = | Dloy; = y €H, (Rad) in this section follows the definition in 
the linear matrix analysis of the previous section, and does not represent the real bunch 
length at the radiator considering the nonlinear modulation waveform. Note also that 
considering the nonlinear sine modulation waveform, the optimal microbunching 
condition for a specific harmonic is slightly different from our simplified linear 
analysis, and J, (n) in Eq. (3.44) should be replaced by J, (—1 Rs6h). In the following 
discussions, we will use the simplified optimal bunch compression conditions, as the 
main physics is the same. 


For a pre-microbunched beam, em [a- Dean) will be non-zero for all m, 
thus 


bp, pre-microbunch = ( 5 Jm (n) exp [— ((n — m)kyo; (Mod))? / a) 3.45) 


m=— 00 


exp |- (nkro; (Rad)}? /2] , 


with o,(Mod) = (d’y; + zi) = af €yH, (Mod) + €,8,(Mod). Note that the bunch 
length o, (Mod) here contains contribution from both €, and €,. 
Now we first investigate two limiting cases. If o,(Mod) = 0, then we have 


oe) 
bn, = ( X Jn o) exp [— (nkro: (Rad))? /2] = exp [— (nkro: (Rad)}? /2] . 
m=— OO 
(3.46) 
This result is the same as that assuming the modulation waveform is linear. Second, if 
o; (Mod) is much longer than the modulation laser wavelength, i.e., ko; (Mod) > 1, 
then the summation terms in Eq. (3.45) will be nonzero only for m = n and we arrive 
at the same result as the coasting beam case Eq. (3.44) as expected. 
Now we conduct a bit more general discussion. Compared to the linear modulation 
case, the reduction factor of the bunching factor Eq. (3.45) is 


Ri = J. Jn (n) exp [= ((n — m)kro:Mod))}? /2]. (3.47) 


m=—oCOo 
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Figure 3.5 shows the flat contour plot for the bunching factor reduction factor | R,| of 
Eq. (3.47) as a function of the harmonic number n and the modulation wavelength- 
normalized bunch length at the modulator kz o; (Mod). As can be seen from the figure, 
the bunch lengthening at the modulator indeed degrades the bunching factor at the 
radiator, due to the nonlinearity nature of sine modulation. The longer this bunch 
lengthening, the more degradation of the bunching factor. The higher the harmonic 
number, the more significant the impact is. The limit of R, with an infinite long 
o,(Mod) is J, (n). Equation (3.47) and Fig. 3.5 is the general result of this bunching 
factor degradation analysis. We emphasize the fact that the discussion of bunching 
factor degradation is more relevant for a pre-microbunhed beam, like that in some 
SSMB scenarios, and is generally not an issue for a coasting beam where the bunch 
duration is much longer than the modulation wavelength, like that in an FEL. 

As the decrease of bunching factor originates from the nonlinearity of the sine 
modulation, we expect that this reduction will be less if we make the modulation 
waveform more like linear, for example by adding a third-harmonic RF or laser 
to broaden the linear zone of the modulation waveform, as also suggested before 


in Refs. [17, 18]. The energy modulation then becomes 6 = 6 + a sin(kyz) + 
TE sin(3kzz). The optimized bunch compression conditions for ADM are now 
(hı + h3)Rs6 + 1 = 0 and —d’D + R56 = 0. The n-th laser harmonic bunching fac- 


tor at the radiator is then 


hı h3 n 
Dn coasting = > Im, (an) Jin; Ger 5) exp |- (nkro; (Rad))? /2] 


mı+3m3=n 
(3.48) 
for a coasting beam, and 
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Fig. 3.6 Left: the bunching factor reduction factor | R,,| of Eq. (3.50) as a function of the harmonic 
number n for kroz (Mod) = 1, with h3 = 0 (red) and h3 = —0.15/ (blue), respectively. Right: the 
bunching factor reduction factor |R| of Eq. (3.50) as a function of the modulation wavelength- 
normalized bunch length at the modulator kz oz (Mod) for n = 79, with h3 = 0 (red) and h3 = 
—0.15h, (blue), respectively. n = 79 corresponds to the case for example a modulation wavelength 
of A, = 1064 nm and a radiation wavelength of Ar = A, /79 = 13.5 nm 


0O OO 


hy hy n 
bpem = D D Jn GG") ms (Gas) 


m,=—0O m3=—00 
exp [- ((n — mı — 3m3)ky,0;(Mod))? /2] exp [- (nkro; (Rad))? /2] 
(3.49) 


for a pre-microbunched beam. Therefore, the reduction factor of the bunching factor 
Eq. (3.49), compared to the linear modulation case, is now 


(oe) CO 


hy h3 5) 
R, = Jm —n | Jin TR E 
2 2 i Ga ) > Ge (3.50) 


m;=—00 M3=—00 


exp [— ((n — mı — 3m3)k,.o,(Mod))* /2]. 


The limit of R, with an infinite long o-(Mod) is Yo, .3sn Jm (zmn) Jn, 


( re a z). It is straightforward to generalize the above derivation and result to the 
case of adding more laser harmonics. 

Now we can use the above formula of R, to do comparison between the cases with 
and without the third-harmonic laser. If h3 = 0, then Eq. (3.50) reduces to Eq. (3.47). 
As can be seen in Fig. 3.6, indeed addition of a third-harmonic laser is effective in 
mitigating the bunching factor degradation arising from the bunch lengthening at the 


modulator. 
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3.1.3.3 Contribution of Modulators to Vertical Emittance and Scaling 
of Required Modulation Laser Power 


We have stated that the main motivation of applying TLC scheme for bunch com- 
pression in SSMB is to lower the requirement on the modulation laser power Pz. 
This is based on the fact that the vertical emittance €, in a planar x-y uncoupled ring 
is rather small. However, since the modulator in this TLC scheme is placed at a dis- 
persive location, i.e., H, (Mod) Æ 0, therefore quantum excitation at the modulator 
will also contribute to €,. With this consideration taken into account, below we try 
to give a self-consistent analysis of the required modulation laser power Pz in these 
TLC schemes. 

To make sure that the TLC-based bunch compression can repeat turn-by-tun in 
a ring, usually two laser modulators are placed upstream and downstream of the 
radiator, respectively, to form a pair. The lattice scheme between these two modulators 
can either be a symmetric one, or areversible seeding one. In both cases, the chromatic 
function H; at two modulators are identical. Assuming the modulator undulator is 
planar, the contribution of these two modulators to €, is then 


Ae, (Mod, QE) = 2 x 


55 arpãy5 f Hy (Mod) 
9/3 av Jo les) 


a (3.51) 

55 ariy? Hy(Mod) 4 

= 2 x - 3 Lu, 
96/3 ov PomMoa 37 
with the vertical damping constant 
alee Eò : C, x 0.2998 Byine[T]E, [GeV] (3.52) 
Ay X x = x U. rin e 
V 2 Eo 2 Y Pring 2 K gs 0 


where C, = 8.85 x 1075 E Pring is the bending radius of dipoles in the ring, and 
PoMoa is the minimum being radius corresponding to the peak magnetic flux density 
Bomoa Of the modulator. Note that in the above analysis we have ignored the contribu- 
tion of the dispersive lattice sections upstream and downstream of the modulators to 
the vertical emittance, as in principle we can minimize their contribution by choos- 
ing weak bending magnets in them. On the other hand, we cannot choose as weak 
modulator as we want since it will also affect the energy modulation efficiency. This 
is the reason why the contribution of modulators to the vertical emittance is of more 
fundamental importance. 

The resonant condition of the laser-electron interaction inside a planar undulator 
is 


2 ius (3.53) 
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with K = Bae. = 0.934 - Bo[T] - A,,[cm] the dimensionless undulator parameter. 
The effective modulation voltage of a laser modulator using a planar undulator is 
related to the laser and undulator parameters according to [19] 


JJIK |4PLZoZ Ey 
v=! l LOOFR tanm! ( l (3.54) 
y ÀL 2Zr 


in which [J J] = Jo(x) — Jı (x) and x = m P; is the modulation laser power, 
Zo = 376.73 Q is the impedance of free space, Zr is the Rayleigh length of the 
modulation laser, L, is the undulator length. The linear energy chirp strength around 
zero-crossing phase is therefore 


V, JJIK [4P,ZZ Lu 
jo BEGET wea kt, (3.55) 
Eo y2mc2 AL 2ZR 


where ky; = 27 /À; is the wavenumber of the modulation laser. 

For simplicity, we set €, = Ae, (Mod, QE), i.e., the vertical emittance is purely 
from the contribution of these two modulators, and assuming that equality holds in 
Theorem one, then the required modulation laser power is 


2 
AL Ey 1 


PL = 
£ — 4ZoZr \ o, (Mod)o, (Rad) eK tan~! (5 )k 
~ y2me2 2Zr L 


2 (3.56) 
= 1 23 55 ape Ry Biya 1 ae 
(LV J]K)* 32°20 48/3 Cy E} Bring 07 (Rad) [tan (zs) 
2Zr 


Now we try to derive more useful scaling laws to offer guidance in our parameters 
choice for a TLC SSMB ring. To maximize the energy modulation, we need Ze = 


0.359 x i. When K > /2,we approximate the resonance condition as Àz © hw 
and [J J] ~ 0.7. Then we have 
3 473 3,8 p3 
P, x Àz Y Bõma 1 m ALY’ Bomoa 1 , (3.57) 
K? Bng +02 (Rad) Bring | 02(Rad) 
The corresponding modulator length scaling is 
Byin Ey 
i a (3.58) 


© are ET 
H, (Mod) Bima 


Putting in the numbers for the constants, we obtain the quantitative expressions of 
the above scalings for practical use 
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P, [kW] ~ 5.67 


’ 


(3.59) 


L,[m] ~ 57 


The above scaling laws are accurate when K > J2. 

Note that €, does not appear explicitly in the scaling of the required laser power. 
It however affects the bunch length at the modulator and therefore the bunching 
factor at the radiator as we have explained. Also it affects the required modulator 
length. In other words, to obtain a desired bunching factor, the smaller €, is, the 
larger o,(Rad) we can use, thus a lower modulation laser power. Generally a shorter 
modulation laser wavelength and lower beam energy is preferred in lowering the 
required laser power. But we need to keep in mind that when the beam energy is 
too low, intrabeam scattering (IBS) could blow up €, [20, 21]. From the scaling, 
a weaker Bomoa Means a smaller modulation laser power will be neadan, But we 
should be aware that the corresponding length of modulator L,, 


of the final chapter, we have presented an example parameters set re a TLC SSMB 
storage ring for high-power EUV and soft X-ray radiation generation, based on the 
investigations presented here. 


3.1.4 Transverse Deflecting RF or TEMOI Mode Laser 
for Coupling 


Now we investigate in more detail about the application of TLC in SSMB for bunch 


compression, using a TEMO1 mode laser modulator for the modulation kick. This 
belongs to the category of Theorem two. 


3.1.4.1 TEM01 Mode Laser Modulator for Bunch Compression 
A laser modulator implementing a TEMO1 mode laser is like a transverse deflecting 


RF cavity in the optical wavelength range. The electric field of a Hermite-Gaussian 
TEMO1 mode laser polarized in the vertical direction is [19] 


2 
ikzziw 1 „kQ A 
E, = Epe (=) apli 2 (0 +y | as 


2 
k 272 fi 
=) exp ee (x + | ae (; = oy’) , 


(3.60) 
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with Q = ee) . The relation between Fo and the laser peak power fora TEMO1 
R lz 
mode laser is given by 
E? ZRÀL 
= 3.61 
L Z, (3.61) 


Note there is a factor of two difference in the above laser power formula compared to 
the case of a TEMOO mode laser. The electron wiggles in a vertical planar undulator 
according to 


K ., 
y@= sin(kyZ), (3.62) 
yk 


u 


and the laser-electron exchanges energy according to 


dW 

— = ev, E, + ev, E. (3.63) 
dt T 

Assuming that the laser beam waist is in the middle of the undulator, and when 
x, y < w(z), which is typically the case in SSMB, we drop exp [i at (x? + y’)| in 
the laser electric field. Since usually Zp >> 4,, we can also drop the contribution 
from E, on energy modulation. Then the integrated modulation voltage induced by 
the laser in the undulator is 


2) (3.64) 


The linear energy chirp with respect to y introduced is then 


L 
Vnoal 2eK[JJ]k, [PiZo 4 
t= = al ER eS A (3.65) 
o y y?mc 1 Lı 
1+ (4) 


Note that the symplecticity of the dynamics requires that the vertical angle of the 
particle after modulation will depend on its initial longitudinal location. This obser- 
vation is also supported by the Panofsky-Wenzel theorem [22] 


Jay 2 2 (=) (3.66) 
as dyo\ y J’ 


where Ay’ and Ay are the electron angular kick and energy change in the laser 
modulator. It is interesting to note that the modulation kick strength depends on the 
ratio between Zr and L,,, instead their absolute values, and the maximal modulation 
is realized when Zr = Ha, 

Now we can do some evaluation based on the formulas. For an example choice of 
parameters, Ey = 400 MeV, àz = 270 nm, K[JJ] = 5, P = 1 MW, Zg = La, 
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we have t = 4 m™!. If e, = 1 pm and the desired bunch length is 3 nm at the 
radiator, i.e., o,(Rad) = 3 nm, then according to Theorem two we have f, (Mod) > 
ate ae = 6.9 x 10° m. The application of TLC with a TEMO1 mode laser to 
compress the bunch length therefore faces the issue of a too large £, at the modulator, 
if the desired bunch length is at nm level, which is needed for example in 13.5 
nm coherent EUV radiation generation. However, if our target wavelength region is 
Ar Z 100 nm, then the idea looks appealing as the required £, is within the reasonable 
reach then. We remind the readers that there could be other more effective angular 
modulation scheme invented such that the issue of large £, at the modulator can be 
solved, even if our target radiation wavelength is still in EUV. 


3.1.4.2 Contribution of Modulators to Vertical Emittance 


As said the advantage of TLC for bunch compression is based on a small vertical 
emittance. Like before let us now investigate the contribution of modulators to €y. 
We remind the readers that for bunch compression using a TEMO1 mode laser, in 
principle, we can place the modulator at a dispersion-free location. Note however, 
if we aim at a complete y-z emittance exchange, the modulator needs to be placed 
at a dispersive location as will be discussed in next section. Here to minimize the 
contribution of modulators to €,, we choose to place the modulator at dispersion-free 
location, which means d = 0 and d’ = 0, then the bunch compression condition is 


T55 = tR54+1=0, 


(3.67) 
T56 = Rs6 = 0. 

Although we have placed the modulator at a dispersion-free location, there is still 
some residual contribution to €, since the transfer matrix of a TEMO1 mode laser 
modulator is intrinsically transverse-longitudinal coupled, and the physical length 
L, and rs54 = 2N,A, of the modulator are nonzero. The thick-lens transfer matrix 
of the laser modulator can be obtained by slicing the laser modulator to tiny slices 
and use the thin-lens kick and drift method to get the total map. If we consider only 
terms to first order of t, rse and L,, then the thick-lens matrix of the TEMO1 laser 
modulator is i in Sates 

u 6 


Ae a0 
56 
O 1lot™ 


Max | ., > 3.68 
2 “sot ceils 1 rs6 ( ) 
t & 0 1 


2 


Note that M; is symplectic to first order of t. The transfer matrix of the state vector 
from the modulator entrance to a distance of s in it is 
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2 t 3 th 
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Here for simplicity we have assumed that the laser is a plane wave such that the 
induced angular modulation strength is proportional to the distance traveled inside 
the modulator. 

Assuming that the one-turn map observed at the entrance of modulator is 


cos ®, + ay sin ®, By sin ®, 0 0 
TO) = —yy sin ®, cos ®, — ay sin ®y 0 0 
0 0 cos ®, + a, sin ®, Bz sin ®, , 
0 0 —yz sin ®, cos ®, — a; sin ®, 


(3.70) 
in which ®, = 27x v, and ®, = 27v,. The eigenvector of the one-turn map corre- 
sponding to the vertical mode at the position s inside the modulator is then 


B+ (i) ba Fe 
i—@y 


1 AV By i 
E; (s) = ME; (0) = —= 2 ; felt, 
J/2 (+) “set By + (+) rsefhy T 
Ly 3 


By 


(3.71) 
Note that to ensure the TLC-based bunch compression can repeat turn-by-tun in 
a ring, usually two laser modulators are placed upstream and downstream of the 
radiator, respectively, to form a pair. According to Chao’s SLIM formalism [23], we 
can calculate the contribution of the two modulators to €, 


Ae, (Mod, QE) = 2 x 


vP 5 Ly E 2 
55 apacy f |E775(s)| ds. (3.72) 
0 


48/3 ay le(s)|? 


When N, > 1, due to the fast oscillating behaviour of sin [2N (+)]. we can 
adopt the approximation 


Ae, (Mod, QE) ~ 2 x 


55 aray? 1 Ja 
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For simplicity, we assume œ, = 0 at the modulator entrance, and usually £, (Mod) > 
La, then 


55 ap&y>? 1 4 r&b 


Ae, (Mod, QE) ~ 2 x 
: 96V3 av Puar 20 


Ly. (3.74) 


Now let us put in some numbers to get a more concrete feeling. For exam- 
ple, if Eo = 400 MeV, pring = 1 m (Bring = 1.33 T), A, = 270 nm, A, = 4 cm, 
K = 3.8, Bo = 1.02 T, N, = 10, rse = 2N, Az = 5.4 wm, Ly = 0.4 m, €, = 1 pm, 
By(Mod) = 100 m, o,(Mod) = ,/e,8,(Mod) = 10 um, o,(Rad) = 2 nm, t= 

! = 9 = 50 m~!, then the contribution of the two mod- 


/ By (Mod)H, (Rad) ay (Mod)o; (Rad) 


ulators to €, is Ae, (Mod, QE) ~ 2.06 fm. So generally, the contribution of the two 
modulators to €, is a small value, if the modulators are placed at dispersion-free 
locations. 


3.1.5 Emittance Exchange 


3.1.5.1 Lattice Condition 


For completeness of the investigation, it might also be helpful to make a short dis- 
cussion on the relation between our TLC analysis and the transverse-longitudinal 
emittance exchange (EEX). For a complete EEX, we need the transfer matrix of the 


form 
0B 
T= k 0) 3 (3.75) 


Therefore, EEX is a special case in the context of our problem definition of TLC- 
based bunch compression, i.e., in EEX the final beam is also y-z decoupled. As 
can be seen from Eq. (3.31), the application of a normal RF or TEMO0 mode laser 
modulator cannot accomplish a complete EEX, as Tg5 = h 4 0. PEHG and ADM 
can thus be viewed as partial EEXs. In contrast, a transverse deflecting RF or TEMO1 
mode laser modulator can be used to obtain a complete EEX. All we need is to add 
another condition to Eq. (3.22), i.e., 


dt+1=0. (3.76) 


After some straightforward algebra, the relations in Eqs. (3.22) and (3.76) can be 
summarized in an elegant form as follows [24] 
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Fig. 3.7 Application of two transverse-longitudinal emittance exchangers to manipulate the bunch 
length in a storage ring 


t = —, 
d 


D = Rad’ + Rzzd, 
D' = Raad' + Rad. 


(3.77) 


Note that the above relations mean that the lattices upstream and downstream the 
transverse deflecting RF are not mirror symmetric with respect to each other [25]. 
Under the conditions given in Eq. (3.77), we have 


0 0 Res d R33 R34 "4 
aes 0 0 Ru ARa — Rag" 
dra3 — r33 Resta drag — 134 Reetaa 0 0 
i ia 0 0 
d d (3.78) 
_| 0 0 Tss T% 
~ | Tss Ts, 0 0 
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3.1.5.2 Two EEXs as an Insertion 


In order to apply EEX to generate short bunch in a storage ring on a turn-by-turn 
basis, another EEX might be needed following the radiator to swap back the €, and 
€, to maintain the ultrasmall vertical emittance €,. If there is only one transverse- 
longitudinal EEX, the ring will then be a transverse-longitudinal Möbius accelerator 
[26], which is also an interesting topic we are not going into in this dissertation. 

Now we consider the application of two y-z EEXs for bunch length manipulation 
in a storage ring as shown in Fig. 3.7. The motivation is still to make use of the fact 
that the vertical emittance €, is rather small in a planar x-y uncoupled ring. The first 
natural idea is to add an inverse EEX unit following the EEX, 


0 O Te —Ts4 
0 0 -Tea Ts 
Ti, -Ts 0 0 i 
-Tis T35 0 0 


T! = (3.79) 


then the total insertion will be an identity matrix and be transparent to the ring. 
The issue of this approach, however, is that we need to design the downstream 
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beamline with an R56 having opposite sign to the upstream beamline, which might 
be a challenging task if we aim at a compact lattice. 

The second natural idea is to implement the mirror symmetry of the upstream 
beamline as the downstream beamline, which is straightforward for the lattice design. 
The transfer matrix of the mirror image is related to that of the original beamline 
according to [27, 28] 


10 00 
0-1 00 
Tinirror UT = U, U = 00 -10 (3.80) 
00 01 
Therefore, 
0 O ~Te4 —Ts4 
T mirror = UT 'U~! = 2 -Tes -Tss (3.81) 


—T4s —T36 0 0 
—Tıs -Tss 0 0 


Note, however, the transfer matrix of the total insertion in this case is generally not 
an identity matrix, 


A’ 0 
M= TinirrorT = ( 0 E) ; (3.82) 
with 
A= — (T53 T64 + Ts4T63) —2T54T64 
—2T53T63 — (T5374 + T5473) } ’ (3.83) 
E= —(T3s T46 + T36 T45) —2T36 Tag l 
—2T35 T45 —(T3s T46 + T36 Tas) j) ` 


A special case of EEX is the phase space exchange (PSX), i.e., the exchange 
happens in the phase space variables apart from a magnification factor. In this case, 
a PSX followed by its mirror can form an identity or a negative identity matrix. 
Case one: 


0 
i 
1 


0 
0 


, T mirror = 


ooos 


0 

0 m 
Om 2 0 
x 0 
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Case two: 
0 0 m 0 o 0 -t 0 
000. 0 0 0 =m 
T= m 0 0 0 » Tmirror = -4 0 0 0 
0-0 0 -m 0 0 
According to Eq. (3.78), for case one, we need 
R34 
= A =0, 
d 
R 
dR + d'R ~ ta 0, 
R 
drs — d'r, a = 0, 
d 
r34 
= 0, 
d 
and 
0 0 0 ¢ 
0 0-2 0 
T= d 
0 0 0 
wÒ 0 0 
~d 
For case two, we need 
rs6R 
d R33 + d' R34 — -A =0, 
_ Ba o 
d ti 
r34 R 
dra — d'ra — ~ = 0, 
d 
Bi 
d kd 
and 
0 0 — 0 
0 0 0 =-& 
T= R34 
_4 9g 0 


75 


» M = ThirrorT = —1. 


(3.85) 


(3.86) 


(3.87) 


(3.88) 


(3.89) 
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3.2 Nonlinear Transverse-Longitudinal Coupling Dynamics 


3.2.1 Average Path Length Dependence on Betatron 
Amplitudes 


After investigating the linear TLC, we will now examine nonlinear coupling. How- 
ever, here we consider only the second-order path lengthening or shortening from 
betatron oscillations, and its impact on equilibrium beam parameters. A general 
discussion of the nonlinear dynamics is beyond the scope of this dissertation. The 
second-order TLC considered here originates from a dependence of the average path 
length on the betatron oscillation amplitudes, which can be expressed by a concise 
formula 


AC = —27 (Ex Jx + y Jy), (3.90) 
where AC is the average path-length deviation relative to the ideal particle, and 
= Tisy and J, y are the horizontal (vertical) chromaticity and betatron invariant, 


respectively. This simple relation is a result of the symplecticity of the Hamiltonian 
dynamics [29-31]. Itis called a second-order coupling because the betatron invariant 
is a second-order term with respect to the transverse position and angle. Note that 
Eq. (3.90) is accurate only for the cases of multiple passes or multiple betatron 
oscillation periods as it is a betatron-phase-averaged result. For the case of a single 
pass with only a few betatron oscillations, there will be an extra term, depending on 
the betatron phase advance, in the path length formula. 

This path length effect has previously been theoretically analyzed by several 
authors in different contexts [29-34]. Due to this effect, particles with different beta- 
tron amplitudes lose synchronization with each other when traversing a lattice with 
nonzero chromaticity. This leads to a stringent requirement on the beam emittance 
for FELs in the X-ray regime (XFELs), as microbunching can be smeared out by 
this effect when the beam is traveling through the undulator [35]. This effect is also 
crucial in non-scaling fixed-field alternating-gradient (FFAG) accelerators for muon 
acceleration [31], as a muon beam typically has a large emittance. Furthermore, the 
natural chromaticities of a linear non-scaling FFAG accelerator are usually not cor- 
rected to achieve a large transverse acceptance. This effect may also have an impact 
on the momentum and dynamic aperture in a storage ring [36, 37], for example, 
due to the Touschek scattering-induced large betatron amplitude or the large natu- 
ral chromaticity in a low-emittance lattice. In this dissertation, we will emphasize 
the importance of this nonlinear TLC effect in precision longitudinal dynamics in a 
storage ring, such as SSMB. 


3.2 Nonlinear Transverse-Longitudinal Coupling Dynamics 77 


3.2.2 Energy Widening and Distortion 


As mentioned, this second-order TLC effect can disperse microbunching in XFELs. 
Methods of overcoming this influence are referred to as “beam conditioning”. Several 
such methods have been proposed since the first publication of Ref. [35]. The basic 
idea of these proposals is to compensate the difference in path length through a 
difference in velocity by establishing a correlation between the betatron amplitude 
and the particle energy. In a storage ring, unlike in a single-pass device, the RF cavity 
will “condition” the beam automatically, causing all particles to synchronize with 
it in an average sense through phase stabilization (bunching). This is accomplished 
by introducing a betatron-amplitude-dependent energy shift to compensate for the 
path-length difference arising from the betatron oscillations, 


Ma, (3.91) 


where œ is the momentum compaction factor of the ring, defined in Eq. (2.1). This 
shift will result in the beam energy widening in a quasi-isochronous ring with nonzero 
chromaticity, because different particles have different betatron invariants [34]. This 
widening will become more significant with the decreasing of the momentum com- 
paction. 

Due to the energy shift, there will also be an amplitude-dependent shift in the 
betatron oscillation center at dispersive locations. The shift direction is determined 
by the signs of a, Ex, y and D, y, and the magnitude of the shift is determined by the 
magnitudes of Jy y, @, Ex, y and D,,. The physical pictures of the betatron center 
shift resulting from this effect are shown in Fig. 3.8. 

When quantum excitation is also taken into account, the total relative energy 
deviation of a particle with respect to the ideal particle is 


Closed orbit 


y 


x Amplitude-dependent \_ á 


| center shift 


Fig. 3.8 Physical picture of the amplitude-dependent shift of betatron oscillation center in the case 
of a positive momentum compaction. Only a horizontal betatron oscillation is considered in this 
illustration 
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TEA Edie, (3.92) 


where dge represents the quantum excitation contribution. Finding a general ana- 
lytical formula for the steady-state distribution of the particles is a complex task 
and, at the same time, not very useful. Simpler expressions can be obtained based 
on reasonable approximations. Since the vertical emittance is usually much smaller 
than the horizontal emittance in a planar uncoupled ring, here, we consider only the 
contribution from the horizontal emittance. When the coupling is not very strong, 
the distributions of J; and dge are still approximately exponential and Gaussian, 
respectively, and are independent of each other, 


de 


l _& 1 = 
(J) = ———-e “o, (Se) = e o, (3.93) 
¥ 20 Ex0 á 1 V 27 050 


where €,9 and oso are the natural horizontal emittance and energy spread. The dis- 
tribution of ô is thus an exponentially modified Gaussian because it is the sum of an 
exponential and a normal random variable, 


AJ *Q25o-28) sen(A) (Ao2, — ô 
metl | ee) ens) | (3.94) 
2 2050 
where A = oe sgn(x) is the sign function and erfc(x) is the complementary error 


function, defined as erfc(x) = 1 — erf(x) = zx J i e-' dt. The direction of long 
non-Gaussian tail of the energy distribution is determined by the signs of œ and &,. 

Because of the dispersion and dispersion angle, the non-Gaussian particle energy 
distribution can also be reflected in the transverse dimension. When this nonlinear 
coupling is considered, the horizontal position and angle of a particle in the storage 
ring are 


2r Ex Jy 
x = 42]; By COS py + Dy | dge + : 


C 
S (3.95) 


x! = —/2J / bx (a, cos Yy + sin yx) + D! (à + = *) ; 


aCo 


It is assumed that the concept of the Courant-Snyder functions is still approximately 
valid in Eq. (3.95). 
With these approximations, the variance of 6 is then 


27 €x0 Fy 4 
of = of, + ( noes) l (3.96) 
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Table 3.1 Parameters of the MLS in the experiment 


Parameter Value Description 

Co 48 m Ring circumference 

Eo 630 MeV Beam energy 

St 500 MHz RF frequency 

Vit <600 kV RF voltage 

Uo 9.1 keV Radiation loss per turn 

Js 1.95 Longitudinal damping partition 

Ts 11.4 ms Longitudinal radiation damping time 
€x0 250 nm Horizontal emittance 

030 4.4 x 1074 Natural energy spread 


By assuming the MLS parameters shown in Table 3.1 and applying a = 1 x 1074 
and £, = 2, one can find that the energy spread contributed by this effect can be as 
significant as its natural value. 
As discussed in Ref. [34], a shift in the energy center corresponds to a shift in the 
synchronous RF phase ¢,, 
Ad, © J; tang, Ad, (3.97) 


where J; is the longitudinal damping partition number and nominally J, ~ 2. There- 
fore, particles with different betatron amplitudes will oscillate around different fixed 
points in the longitudinal dimension, thus lengthening the bunch. The change in the 
synchronous RF phase in a unit of the longitudinal coordinate, Az,, is related to the 
relative change in energy, Ad, according to 


AZs Vs Js tan Øs Ad 1 Aô 
— xX x 
oz  fre/frevlæl oso = Va] oso 


(3.98) 


where v, is the synchrotron tune, frey is the particle revolution frequency in the ring, 
09 and oso are the natural bunch length and energy spread, respectively. 

The critical value of alpha, œ., when the relative change of bunch length and 
energy spread are the same can be calculated to be 


vs J; tan bs Js To tan ds 

De gja a (3.99) 
Fre/frev|Oc| NT fRE/ frevTs 

where ts = 1/az, = a To is the longitudinal radiation damping time. As an exam- 


ple, we use the MLS parameters given in Table 3.1 and consider the application 
of an RF voltage of 500 kV, which corresponds to a synchronous RF phase of 
$s = 0.018 rad. The critical value of alpha is then |a,| ~ 2.1 x 107°, which is about 
four orders of magnitude smaller than the alpha value reachable at the present MLS. 
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Fig. 3.9 Energy widening, bunch lengthening and distortion from a Gaussian distribution induced 
by a nonvanishing horizontal chromaticity. From up to bottom, the particle tracking results for 
distributions of x, x’, z and ô are shown at two dispersive locations in the MLS under three different 
horizontal chromaticities £x. The direction of the long non-Gaussian tail for ô is related to the signs 
of œ and éx, while for x and x’ they are also dependent on Dy and D‘., respectively. The simulation 
was conducted using the code ELEGANT [38] with a beam energy of 630 MeV, an RF voltage of 
500 kV and the application of œ = 1 x 1074. In each simulation, eight particles were tracked for 
5 x 10° turns, corresponding to approximately 73 longitudinal radiation damping times 


Therefore, the relative bunch lengthening resulting from this effect is much less 
significant than the corresponding energy widening at the MLS. 

Several particle tracking simulations were conducted using the MLS lattice with 
the parameters presented in Table 3.1 to confirm the analysis. Two dispersive loca- 
tions, with different signs and magnitudes of D,, were selected as the observation 
points in the simulations. The simulation results are shown in Fig. 3.9. The energy 
widening and distortion from Gaussian behaviors are as expected and, indeed, are 
more significant than the bunch lengthening when a = 1 x 1074. At the two observa- 
tion points, widening and distortion of the particle energy distribution also manifest 
in the transverse dimension through D, and D‘,. The related optic functions at the two 
observation points are also shown in the profiles of x and x’. Note that the directions 
of the long non-Gaussian tails of the profiles and their relations to the signs of £y, 
D, and D'.. We conclude that the simulation results agree well with the analysis and 
physical pictures presented above. 


3.2.3 Experimental Verification 


Here we report the first experimental verification of the energy widening and particle 
distribution distortion from Gaussian due to this second-order TLC effect as analyzed 
above. At the MLS, the Compton-backscattering (CBS) method is applied to measure 
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Fig. 3.10 Measurements of the CBS photons spectra at the 344.28 keV (a) and 778.90 keV (c) 
emission lines of !°*Eu radionuclide to calibrate the channel numbers in terms of keV and the 
HPGe-detector resolution at the CBS cutoff edge (b) in the CBS method of measuring beam energy 
spread 


the electron energy [39]. Nevertheless, within certain limitations, the electron beam 
energy spread can also be evaluated from the CBS photon spectra [40]. The non- 
Gaussian momentum distribution makes the evaluation a bit more involved, but we 
can assume a Gaussian distribution with an equivalent mean energy spread. This is 
a good approximation as long as &, is not too large. 

The experiment is conducted with all 80 RF buckets equally filled. To exclude a 
severe impact from energy widening collective effects, the beam current is decayed 
till the horizontal beam size is not sensitively dependent on it. The average single- 
bunch current is below 12.5 pA (1 electron/1 pA) while doing the CBS measurements. 
To mitigate the influence of a nonlinear momentum compaction, the longitudinal 
chromaticity has been corrected close to zero. The other parameters of the ring in 
the experiment are presented in Table 3.1. 

To get the energy spread based on the CBS method with precision, the HPGe- 
detector used in the measurement should be calibrated in terms of the photon energy 
per channel. This is realized by recording the emission lines from a !>7Eu radionuclide 
simultaneously during the measurement of the CBS photons. Moreover, the width 
of the fitted '°*Eu lines that are close to the edge of the CBS photons have been 
used to determine the detector resolution daet at the photon energy of the CBS cutoff 
edge, Eeage, in our case 707 keV. This is done by a linear interpolation of the width 
of the !?Eu lines at 344.28 and 778.90 keV. The detector resolution daet at 707 keV 
is thus determined to be 0.64(4) keV. The calibration scheme and result is shown in 
Fig. 3.10. 

Figure 3.1 1a shows the typical CBS photon spectra close to the cutoff edge under 
the cases of different £, and a. The adjustment of £y is accomplished by the imple- 
mentation of different chromatic sextupole strengths and the œ by slightly tuning 
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Fig. 3.11 Measurement of the electron beam energy widening brought by the horizontal chro- 
maticity using the CBS method. a The cutoff edges of CBS photon spectra under different &, and 
fs (therefore œ). b Quantitative evaluation of the cut-off edges revealing the energy spread and its 
comparison with theory Eq. (3.96). The error bars in both figures are the one sigma uncertainties 
of the measurements and are due to calibration errors and counting statistics. The data acquisition 
of each spectrum takes 15 min 


quadrupole strengths. In the experiment the RF voltage is kept constant and the syn- 
chrotron frequency f; is proportional to the square root of the magnitude of a. The 
edge in the figure is a convolution of a step function representing the CBS cutoff 
edge with a Gaussian function which attributes to the finite HPGe-detector energy 
resolution and the electron beam energy spread. The fitted line is basically an error 
function from which the energy width of the CBS photons at the edge Geage, and 
therefore the electron beam energy spread o3, can be deduced. It is assumed Gedge 


is given by edge = J A + (2E cages). The second term in the square root is due 
to the electron beam energy spread and is based on the fact that the energy of the 
backscattered photon is proportional to the electron energy squared. 

It can be seen from Fig. 3.1 1a that the edge slope decreases with the magnitude 
increasing of £, and lowering of œ when &, 4 0, which indicates that there is an 
energy widening in the process. Quantitative evaluation of the edges revealing the 
energy spread and its comparison with theory of Eq. (3.96) are shown in Fig. 3.1 1b. 
The energy spread grows significantly with the magnitude decrease of œ when &, = 
—3.14 while it stays almost constant in the case of &, = 0. The agreement between 
measurements and theory is quite satisfactory. This is the first direct experimental 
proof of the impact of this effect on the equilibrium beam parameters in a storage 
ring. 
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Fig. 3.12 Transverse beam intensity distortion from Gaussian at dispersive locations due to a non- 
vanishing horizontal chromaticity measured by the imaging systems installed at QPDO and QPD1. 
Three different £y are applied in both the positive and negative momentum compaction modes with 
a = +8.4 x 1075. There is some residual horizontal-vertical coupling in the positive momentum 
compaction case, which do not influence the principle observation of the non-Gaussian behavior 


As analyzed before, the bunch lengthening due to this second-order coupling 
is much less notable and also due to the limited resolution of the present streak 
camera, we do not measure the bunch lengthening in the experiment. Nevertheless, 
a more comprehensive investigation of this effect can be conducted on the other 
beam characteristics like the transverse intensity distribution. As can be seen from 
Eq. (3.95), particles with different betatron amplitudes oscillate around different 
closed orbits, which is the amplitude dependent center shift [41]. Because of the 
dispersion, the non-Gaussian particle momentum distribution can also be reflected 
to the transverse dimension, which can be observed by the beam imaging systems 
installed at the MLS [42]. 

Figure 3.12 shows the typical transverse beam intensity distribution measured by 
the imaging systems at two dispersive locations, QPDO and QPD1, with different 
values of é, in both the negative and positive momentum compaction modes. The 
relevant optics functions, 6, and Dx, at the two observation points are also shown 
in the figure. Note that D, have different signs and magnitudes at QPDO and QPD1. 
It can be seen that the horizontal beam distribution at these dispersive locations 
becomes asymmetric when éx Æ 0. The long tail direction and the magnitude of 
deviation from Gaussian are determined by the signs and magnitudes of a, Dy, Ex 
and also the value of €,9 and 6,, which fits with the expectations. 

Figure 3.13a demonstrates the typical horizontal beam profiles measured at QPD1 
in the negative momentum compaction mode under three different £, and their good 
agreements with theory. It turns out that both the theoretical and experimental mea- 
sured horizontal coordinate distribution w(x) can be excellently fitted by a skewed 
Gaussian function 


w(x) = : E (1 +ert]a- =") (3.100) 
” Sino 
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Fig. 3.13 Horizontal beam profile distortion from Gaussian by horizontal chromaticty. a Typical 
horizontal beam profile at QPD1 with a = —7 x 1075 under three different £. The closed orbit 
movements of the ideal particle due to the sextupole strengths changes when adjusting £y have 
been compensated in the plot. Cross: beam imaging system measurement results. Dashed line: 
fit of the measurement data by an exponentially modified Gaussian function Eq. (3.100). Solid 
line: theoretical prediction. b Measured and theoretical asymmetry parameter d versus £x at QPDO 
and QPD1 with fs = 5 kHz (a = —7 x 1075); c Measured and theoretical asymmetry parameter 
d versus fs at QPDO and QPD1 with £ = 1.4. All the theoretical curves are obtained based on 
Eggs. (3.93) and (3.95) 


The asymmetry parameter d in Eq. (3.100) is used to quantitatively describe the 
deviation from Gaussian and as a criterion to do comparison between measurements 
and theory. Figure 3.13b and c show the asymmetry parameter d versus &, and fy, 
therefore a, from measurements and theory at QPDO and QPD1. It can be seen that 
the larger the £, and the smaller the a, the more asymmetric the distribution is. Also 
the asymmetry at QPD1 is more significant than that at QPDO as the magnitude of Dy 
at QPD1 is larger while the 6, difference at two places is not much. The agreement 
between measurements and theory confirms that this effect distorts the beam from 
Gaussian in both the longitudinal and transverse dimensions. 

While the energy widening and beam distortion could be a detrimental outcome 
for some applications, it may actually also be beneficial as it can help to stabilize 
collective instabilities. The bunch lengthening on the other hand is much less notable 
compared to the energy widening. So quasi-isochronous ring-based coherent radia- 
tion schemes, like some of the SSMB scenarios, may boost the stable coherent radi- 
ation power by taking advantage of this effect. For example, the stable single-bunch 
current at the MLS can grow for more than one order of magnitude by increasing the 
absolute value of the horizontal chromaticity from zero to a value larger than three, 
with the head-tail and the other collective effects like the longitudinal microwave 
instability properly suppressed. It has been proved at the MLS that the increase of 
THz power due to a higher stable beam current overcompensates the decrease due 
to the slight bunch lengthening of this effect. Therefore, this is now the standard low 
momentum compaction mode at the MLS for the application of Fourier Transform 
Spectroscopy. 

This nonlinear TLC may also be useful in some more applications. For example, 
it can be used for the real-time emittance evaluation in storage rings if the chro- 
maticities, beta function and dispersion are known, which are usually easier to get 
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than measuring the emittance directly. The amplitude dependent center shift can be 
applied to detect beam instabilities which blow up the transverse emittance [43]. A 
strongly asymmetric particle momentum distribution due to this effect cooperating 
with a large momentum compaction lattice can generate a strongly asymmetric dis- 
tributed current, which is favored in some applications such as beam-driven wakefield 
acceleration [44]. 
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Chapter 4 ®) 
SSMB Radiation E 


Having discussed the methods to form and preserve microbunching in the last two 
chapters, now we present the theoretical and numerical study of the average and sta- 
tistical property of coherent radiation from SSMB. Our results show that a kW-level 
average power quasi-continuous-wave EUV radiation can be obtained from an SSMB 
ring, provided that an average current of 1 A and bunch length of 3 nm microbunch 
train can be formed at the radiator which is assumed to be an undulator. Together 
with its narrowband feature, the EUV photon flux can reach 10!5 ~ 10!° phs/s within 
a 0.1 meV energy bandwidth, which is three orders of magnitude higher than that in 
a conventional synchrotron source, allowing sub-meV resolution in angle-resolved 
photoemission spectroscopy (ARPES) and providing new opportunities for funda- 
mental physics research. In the theoretical investigation, we have generalized the 
definition and derivation of the transverse form factor of an electron beam which can 
quantify the impact of its transverse size on coherent radiation. In particular, we have 
shown that the narrowband feature of SSMB radiation is strongly correlated with the 
finite transverse electron beam size. Considering the pointlike nature of electrons 
and quantum nature of radiation, the coherent radiation fluctuates from microbunch 
to microbunch, or for a single microbunch from turn to turn. Some important results 
concerning the statistical property of SSMB radiation are presented, with a brief dis- 
cussion on its potential applications for example the beam diagnostics. The presented 
work is of value for the development of SSMB and better serve the potential syn- 
chrotron radiation users. In addition, it also sheds light on understanding the radiation 
characteristics of free-electron lasers (FELs), coherent harmonic generation (CHG), 
etc. Parts of the work presented in this chapter have been published in Ref. [1]. 
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4.1 Formulation of Radiation from a 3D Rigid Beam 


For simplicity, as the first step we consider only the impacts of particle position x, 
y and z, but ignore the particle angular divergence x’, y’ and energy deviation ô, 
on the radiation. Under this approximation, concise and useful analytical formulas 
of the coherent radiation can be obtained. This approximation is accurate when the 
transverse and longitudinal beam size do not change much inside the radiator, i.e., 
By.y Z Lr and B, Z Rs6,r, where x,y,z are the Courant-Snyder functions of the beam 
in the horizontal, vertical and longitudinal dimensions [2], L, and R56,- are the length 
and momentum compaction of the radiator, respectively. Here in this dissertation, we 
call this approximation the three-dimensional (3D) rigid beam approximation, as the 
beam sizes do not change much during radiation. We will see later in Sect. 4.6 that 
the conditions of rigid beam approximation is generally satisfied in the envisioned 
EUV SSMB. In addition, we will briefly discuss the impact of beam divergence and 
energy spread on coherent radiation in Sect. 4.4. 

Assuming that the vector potential of radiation from the reference particle at the 
observation location is Apoint(@, Ø, t), with @ and @ being the polar and azimuthal 
angles in a spherical coordinate system, respectively, as shown in Fig. 4.1. Under 
far-field approximation, the vector potential of radiation from a 3D rigid electron 
beam containing N, electrons is then 


Fig. 4.1 Coordinate system Zz 
used to calculate the 
undulator radiation 
spectrum. The magnetic field 
is in the y-direction, and the 
electron wiggles in the x-z 
plane 


Observation 


Electron 
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o0 o0 o0 
Abeam(9, y,t) = ne f / / 
—oo J—co J—00 


x sin@cosg+ y sin 8 sing Z 
A point (0.004 ; F zz p(x, y, z)dxdydz, 


(4.1) 


in which £ is the particle velocity normalized by the speed of light in vacuum c, 
and p(x, y, z) is the normalized charge density satisfying [°° SX S2 p(x. y, z) 
dxdydz = 1. Note that we have assumed that the particle motion pattern, therefore 
also the radiation pattern of a single electron, does not depend on x, y and z of 
the particle. In other words, x, y, z of a particle influences only the arrival time of 
the radiation at the observation. This is the reason why their impacts can be treated 
within a single framework. The impacts of x’, y’ and ô are different. Generally, their 
impacts are two-fold. First, they affect the radiation of the single particle itself, i.e., 
the radiation pattern. Second, they affect the electron beam distribution, therefore 
the coherence of different particles, during the radiation process. 
According to the convolution theorem, for a 3D rigid beam, we now have 


Abeam(6, Q, w) = NeApoint(9, P, w)b(0, Q, w), (4.2) 


where 


oo oo oo ` xsinð cosø+ysinð sing , z 
b(0, 9, @) = J l / p(x, y, ze iol i Aleide (4.3) 
—oo J —00 J—00 


Since A(0, g, t) is real, then A (8, y, —w) = A* (0, 9, w). The energy radiated per 
unit solid angle per unit frequency interval is [3] 


d'W 
— (0 = 2|A(6 z 4.4 
Joda. , 9#, ©) = 2|A(O, p, w)| (4.4) 
Therefore, we have 
aw dW 
— (6,9, = N?Ib(0, p, œ|? 6,9, 4.5 
Toda. aT 2 1b, 9, w)| Jodo | Oe (4.5) 


The total radiation energy spectrum of a beam can be obtained by the integration 
with respect to the solid angle 


m 2m dW 
= in d0 do —— (0 
[ sin Í ? Toda. , 9, œ) 


and the total radiation energy of the beam can be calculated by further integration 
with respect to the frequency 


dW 


Ja ; (4.6) 


beam 


beam 
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(4.7) 


beam 


The reason why the lower limit in the above integration is 0, instead of —oo, is that 
there is a factor of 2 in the right hand side of Eq. (4.4). The above formulas can be 
used to numerically calculate the characteristics of radiation from an electron beam, 
once its 3D distribution is given. Note that in the relativistic case, we only need to 
account for 0 several times of 5e as the radiation is very collimated in the forward 
direction. 


4.2 Form Factors 


When the longitudinal and transverse dimensions of the electron beam are decoupled, 
we can factorize b(0, 9, œw) as 


DO, p, w) = b1 (0, 9, w) x bz(w), (4.8) 
where 
oo oo -io(? snfeosptysind sine ) 
bi (0,9, œw) = I f p(x, y)e ž dxdy, (4.9) 
=p — OO 
and Š 
b-(w) = I p(zje mdz. (4.10) 
—0o 


Note that p(x, y) and p(z) are then the projected charge density. b; (œw) is the usual 
bunching factor found in literature and is independent of the observation angle. This 
however is not true for b1 (6, g, œw). For example, in the case of a 3D Gaussian x-y-z 
decoupled beam, 


2 
|b, (6,9, W)? = exp {- (=) [0x sin 0 cosy)? + (oy sin 6 sin %) | ; 
c 


‘ (4.11) 
|b-(w)|? = exp |- (=) a 


where o,,),, are the root-mean-square (RMS) size of the beam in the horizontal, 
vertical and longitudinal dimension, respectively. 

In order to efficiently quantify the impact of the transverse and longitudinal distri- 
butions of an electron beam on the overall radiation energy spectrum, here we define 
the transverse and longitudinal form factors of an electron beam as 
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fo sindo i dg|b_(, Q, w) EX o, p, w) 


FF\(@) = point (4.12) 
Jo sinodo f?” do £40, 9, w) 
point 
and 
FF,(@) = |b-(@)), (4.13) 
respectively. The overall form factor is then 
FF(@) = FF\(@)F F; (w). (4.14) 


The total radiation energy spectrum of a beam is related to that of a single electron 
by 


dW dW 
—| =N2FF(@)—| . (4.15) 
dw 


beam point 


4.2.1 Longitudinal Form Factor 


The longitudinal form factor is the usual bunching factor squared, and have been 
discussed extensively in literature. For example, the longitudinal form factor for a 
single Gaussian microbunch is given Eq. (4.11). When there are multi microbunches 
separated with each other a distance of the modulation laser wavelength Àz like that 
in SSMB, the longitudinal form factor is that of the single bunch multiplied with a 
macro form factor, 


a sin (M24) \" 
F Fm low) = FFs lo) | ———— 3, | > (4.16) 
N; oF) 


where the subscripts mg and sg mean multi bunch and single bunch, respectively, 
and N; is the number of microbunches. This macro form factor of multi bunches is 
a periodic function of the radiation frequency, with a period of the modulation laser 
frequency. The full width a half maximum (FWHM) linewidth around each laser 
harmonics is AWpwum = KA . When N, goes to infinity, this macro form factor will 
become the periodic delta function. Figure 4.2 presents an example plot of the macro 
form factor for three different Np. 

When the radiation wavelength is a high harmonic of the modulation laser wave- 
length, corresponding to these delta function lines in the longitudinal form factor or 
radiation energy spectrum, there will be interference rings in the spatial distribution 
of the coherent radiation from different microbunches. The polar angles of these 
rings, corresponding to the delta function lines in energy spectrum, are determined 
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Fig. 4.2 Macro form factor 
of multi bunches, as a 
function of the number of 
bunches 


by the off-axis resonant condition. Note that when we use the above form factor 
F F;ņyg (œ) to calculate the radiation energy or spectrum from N, microbunches, the 
number of electrons used should be Np Nesp, with Nesp the number of electrons per 
microbunch. We remind the readers that the electron beam energy spread and angu- 
lar divergence will make the linewidth of the delta function lines analyzed above 
become non-zero. For example, the relative bandwidth of the radiation caused by an 
energy spread of o; is 205. 

To make our results more useful, here we also present some analysis applies 
for FEL and CHG. As analyzed in Sect. 2.2.1, the bunching factor for a coasting 
beam-based CHG is 


[0,0] 


w w 1/7% 2 
bz coastin = é6(—— kL) Jn[-=R Alex a [eR ) , 
z,coasting (@) d (- NKL z 856 P| A 5695 
(4.17) 
where 
1, x=0, 
8(x) = 4.18 
) else. ( ) 


Let us now consider the more-often confronted case of a finite bunch length. We 
assume that the initial beam current before microbunching is Gaussian with an RMS 
bunch length of o,. According to the convolution theorem, then 


bz punched (@) = bz coasting (œ) ® Dz Gaussian (a), (4.19) 


. 2 
where ® means convolution and b; Gaussian(@) = exp [-3 (20:) | . Therefore, the 


longitudinal form factor of the beam is now F F; punchea(@) = |b; bunched (@)|?- The 
convolution with b; Gaussian (œ) results in a non-zero bandwidth of each laser harmonic 
line in the longitudinal form factor spectrum 


c/o, — 4/21n2 


(Aowa = 2V2 E = , 
ate J/2 (At) ewum 


(4.20) 
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where At is the electron bunch length in unit of time. Then the relative bandwidth 
of the longitudinal form factor at the H-th laser harmonic can be expressed as 


(=) _ 2721n2 1 O 24/22 A (4.21) 
® J FWHM m — (cAt)rwum/A a  H(cAt)ewum l 
Note that the coherent radiation pulse length is - of the electron bunch length 


J2 
due to the scaling of Peoh « N?, and the above formula means that the coherent 


radiation is Fourier-transform limited. Note also that the absolute width (Aw)rpwum 


is independent of H, while the relative bandwidth (42) an z. 


4.2.2 Transverse From Factor 


Now let us investigate the transverse form factor. Since the transverse form factor 
depends on the radiation process, there is not a universal formula involving only the 
beam distribution. Here for our interest, we focus on the case of undulator radiation. 
We use a planar undulator as an example. The formulation for a helical undulator is 
similar. 

As well established in literature, the planar undulator radiation of a point charge 
in the H-th harmonic is [4] 


aw. 2e?y? 
Zopo] =~" Gé@,y)Fe, 
dodQ A TEC 
point 
sin(x N,€) 2 w w 2 292 

F(e) = ,€= H= 1+K?/2+y76") — H, 

o ( TE ) i œr (0) 2ckuy? ( i í ) 
G(0, p) = Go (0, p) + Gz (0, 9), 

H(£D + 9) cos ) i : 

G, 6,9) = Val © J COS P G,(6 joe Hy@D2 sing 

oO, 9 14 K224 7707 r Sa P= 3 NIF K224 y0) * (4.22) 


1 


ais XO Ja42m—1(He) [Jim (HE) + Im—1(HE)], 


m=— 00 


Dı 


D = X. Jy 42m(Ha)Jm (HS), 


m=—00 
2K y8 cos o K? /4 
a= TE = , 
14 K?/2 + y?0? 1 + K?/2 + y?0? 


in which e is the elementary charge, y is the Lorentz factor, €o is the permittivity of 
free space, œ, (0) is the fundamental resonant angular frequency at the observation 
with a polar angle of 0, ku = A,,/27 is the wavenumber of the undulator, K = eBo 


meks 
0.934 - Bo[T] - 4,,[cm] is the dimensionless undulator parameter, with Bo being the 
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peak magnetic flux density of the undulator and m, being the mass of an electron, 
Jim means the m-th order Bessel function of the first kind. 

Now we try to get some analytical results for the transverse form factor. Instead 
of a general discussion, here we consider only the simplest case of a transverse round 
Gaussian beam, i.e, 


2 
|b (6, 9, w)|? = exp |- (2o, sind) | ; (4.23) 
Cc 


where o is the RMS transverse size of the electron beam. As the radiation is domi- 
nantly in the forward direction in relativistic case, and e~ (Zou sina)’ approaches zero 
with the increase of 6, therefore in Eq. (4.12), the upper limit of 6 in the integration 


can be effectively replaced by infinity, and sin 0 can be replaced by 0 in e~ (Gor sind)” 
Then 


27 m 2 
See d-W. 
Í ap f sin 0dde~ (E21 8in#) Z TE O, 9 wy) 
0 0 dwdQ 


point 


” ao | odo- C0 EYE 4.24 
~ | dof odoo% 0.9. 
f aof oaot 4 0, 0.0) (4.24) 


point 
2 


7 


2n ee) 3 
1 dpG(r, 9) i d(yo)-e 9) F(e), 
0 0 


ITEQC 


where 0; = 1j + K?/2) (4 — 1), wo = œw (0 = 0) = Zro. Here we have 
made use of the fact that there is only one value of 0, i.e., 6,, that contributes signif- 
icantly to the integration over the solid angle © due to the sharpness of F (€) when 
the undulator period number N, > 1, as the spectral width of F (e€) is 1/Ny. 

The transverse form factor corresponding to the H -th harmonic can thus be defined 
as 


B i d(y0)2e~(2219) sinc?(N, 7e) 


FF, (H, œ) = (4.25) 
ree) SE d(yoysinc? (Nz) 
The radiation spectrum of the H-th harmonic is then 
dW dW 
E| =N?FF\(H,o)FF(o)—|__, (4.26) 
dw ie dw : 
eam point 
and the total radiation spectrum of an electron beam is 
dW < dw 
an ZN (4.27) 
dw beam H=1 dw beam 
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Denote 


2 
w w 1 wo, 

=N, —-—H), =N, , = , 4.28 
“(2 ) “2 "o 1 + K2/2 = ( cy ) wae 


then the denominator in Eq. (4.25) is 


Z — Si (2x) + #2 


K2 


(4.29) 


CO 
D(A, w) = f dxsinc? (k, + k2x) = 
0 


where Si(x) = h sin) ay is the sine integral, and the numerator in Eq. (4.25) is 
oo 
N(H, w) = I dxe~*3* sinc? (k1 + K2x) 
0 


#13 


2 K3 
Te [~u sin? (Ke — xai E (2i a) ei( 2ķıi a8)! 
K K2 


= 2 
4k 1 k5 2 


TeK [Fi (200 as) 2i ( as) + Bi ( Dei we), 
K2 K2 K2 
(4.30) 


where Ei(x) = y E “dt is the exponential integral. The transverse form factor is 
then 


N(H, œ) 


FF, (H, w) = DH, o) 


(4.31) 


When w = Hw, then x, = 0, the transverse form factor has a simpler form, 


J dxe™* sinc? (k2x) 


Jy dxsine? (k2x) 
= 2 [0 (3 sin 
= fian 55) + "ICSF! |: 


2 w 2 
K3 o ku% Oo 
SGi,ly.0)= TAR Te = (4.33) 


FF (S) = FF, (H, Hao) = 


(4.32) 


where 


is the diffraction parameter, with L, = N,A, being the length of the undulator. This 
form factor F F; (S) is a universal function and has been obtained before in Ref. [5]. 
Here we have reproduced the result following the general definition of the transverse 
form factor. The variable S is a parameter used to classify the diffraction limit of the 
beam, 

1, S <1, below diffraction limit, 


x45, S 1, above diffraction limit. 


FF\(S) = (4.34) 
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Fig. 4.3 The universal 10° 
function F F, (S) and its 
asymptotic value above 
diffraction limit. The solid 
line comes from Eq. (4.32), 
the dashed line from the 
asymptotic relation above A 
diffraction limit Eq. (4.34) g 10 F 


102- 


This function along with its asymptotic result above diffraction limit are shown in 
Fig. 4.3. 

Note that the decrease of F F, (S) with the increase of o, (S « o?) means that 
the coherent radiation at the frequency w = Haw becomes less when the transverse 
electron beam size becomes larger. This reflects the fact that for a given radiation 
frequency w, there is a range of polar angle 0 that can contribute. For œ = Hao, 
not only 6 = 0, but also 0 very close to O contribute. With the increase of o1, 
the effective bunching factor b(0, gy, œw) at œ = Hwo drops for these non-zero 0 
due to the projected bunch lengthening, therefore the coherent radiation becomes 
less. Another way to appreciate the drop of F F; (S) with the increase of ø, is that 
there is a transverse coherence area whose radius is proportional to ./L,,A9/H with 
ào = 20 a and less particles can cohere with each other when the transverse size 
of the electron beam increases. 

Note that our definition Eq. (4.12) and derivation of the transverse form factor 
Eq. (4.31) is more general than that given in Ref. [5], as it covers other frequencies 
in addition to a single frequency wọ. Therefore, it contains more information than 
Eq. (4.32) as will be presented soon. However, Eq. (4.31) is still not simple enough for 
efficient analytical evaluation to provide physical insight. A further approximation 
is thus introduced, 


i. dxe—** sinc? (K1 + K2x) 
Jy dxsinc? (kı + k2x) 
asa f° dye = sinc?(y) 
en OF Tune 
Je, dysine*(y) (4.35) 


FF\(H,@)= 


aa fo” dye” 2 sinc?(y) 
ee Sn 

Jy dysine*(y) 
oN S(H— 


FF, (S). 
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—Numerical integration —Numerical integration 
—Analytical formula 


—Analytical formula 


0 0.5 1 0 0.5 1 
w/wo w/wo 


Fig. 4.4 The comparison of the transverse form factor, between that calculated from our simplified 
analytical formula Eq. (4.35) and that from the direct numerical integration of Eq. (4.12) for the 
case of H = 1, with N, = 10 (left) and N, = 100 (right), respectively. Other related parameters 
used in the calculation: Eg = 400 MeV, ào = 13.5 nm, A, = 1 cm, K = 1.14, o1 = 5 um 


The condition of applying such simplification is a (w = Hao) X lorS (w = Hao) 
<1, i.e., the beam is below diffraction limit for the on-axis radiation œ = Ha. 
Therefore, the conditions of applying Eq. (4.35) are 


1 Ludo 
N, > lando; < —— 


VHV 2m ` 


If the second condition is not satisfied, the more accurate result Eq. (4.31) should be 
referred. 

As a benchmark of our derivation, here we conduct some calculations of the trans- 
verse form factor based on direct numerical integration of Eq. (4.12), and compare 
them with our simplified analytical formula Eq. (4.35). The parameters used are for 
the envisioned EUV SSMB to be presented in Sect. 4.6, and are given in the figure 
caption. As can be seen in Fig. 4.4, their agreement when N, = 100 is remarkably 
well. Even in the case of N, = 10, the agreement is still satisfactory. There are two 
reasons why the agreement is better in the case of a large N,,. The first is that in the 
derivation we have made use of the sharpness of F (e€), whose width is 1/N,,. The 
second is that S œ N with a given transverse beam size and undulator period length, 
and our simplified analytical formula Eq. (4.35) applies when S (w = Ha) < 1. 

To appreciate the implication of the generalized transverse form factor further, an 
example flat contour plot of the transverse form factor as a function of the radiation 
frequency w and transverse electron beam size o1 is shown in Fig. 4.5. As can be seen, 
a large transverse electron beam size will suppress the off-axis red-shifted coherent 
radiation due to the projected bunch lengthening from the transverse size, thus the 
effective bunching factor degradation, when observed off-axis. Therefore, a large 
transverse electron beam size will make the coherent radiation more collimated in 
the forward direction, and more narrowbanded around the harmonic lines. But note 


(4.36) 
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Fig. 4.5 Flat contour plot of 20 
the transverse form factor 08 
FF (H,w) for H = 1,asa 
function of the radiation 15 0 
frequency w and transverse si 
electron beam size oL, 10 
calculated using our i 0.4 
simplified analytical formula 
Eq. (4.35). Parameters used 0.2 
in the calculation: Ey = 400 
MeV, Ap = 13.5 nm, A, = 1 0 

0 0.5 1 


cm, K = 1.14, N, = 2 x 79 


o, (um) 


A 


w/w 


that not only the red-shifted radiation is suppressed, the radiation strength of each 
harmonic line @ = Hap actually also decreases with the increase of the transverse 
electron beam size, the reason of which we have just explained. 

Now we evaluate the bandwidth and opening angle of the radiation at different 
harmonics due to the transverse form factor. In particular, we are interested in the case 
where the off-axis red-shifted radiation is significantly suppressed by the transverse 
size of the electron beam, which requires that 


NS (w = Hap) > 1, (4.37) 


i.e., o, > ye Vas, Note that to apply Eq. (4.35), we still need the conditions in 
Eq. (4.36). For example, to apply the analytical estimation for the example EUV 
SSMB calculation to be presented in Sect. 4.6, in which A, = 1 cm, Ap = 13.5 nm 
and N, = 2 x 79, we need 1.3 um & o, < 58 pm. The typical transverse electron 
beam size in an EUV SSMB ring is in this range. 

With these conditions satisfied, the value of the exponential factor in Eq. (4.35) 
is more sensitive to the change of w, compared to the universal function F F, (S). 
Therefore, here we consider only the exponential term when w is close to Ha. 
—4N, 7 s(H- 2) 


We want to know the w at which the exponential term gives e Si 


Putting in the definition of S = 7-4, we have we- = ——,——— Hap. Then 
E ee ee 
Fa? kuk 
A@e-! = H wo — We- = Ai d Hw. (4.38) 
L 2 
As o, > ,/ 449, therefore aaa < 1, we have 
Ae- 1 
C) a (4.39) 
Hw |; 2H?ø{ kuko 
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Correspondingly, the opening angle of the coherent radiation due to the transverse 
form factor is 


ye, Awe V2 +K? (4.40) 
7 onl Z —— n 
1 + K?/2 Hoo |; ı 2HyoiJkuko 
It is interesting to note that 
Awe- 1 (4.41) 
x —. : 
H wo ij H? 


As a comparison, the relative bandwidth at the harmonics due to the longitudinal 
form factor is 


Aw 1 
— | & —. (4.42) 
H Wo z H 
Note also that So a | and 6@,-1| are independent of N,, although the approximations 
Wo L 


L 
adopted in the derivation actually involve conditions related to N,. 
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In many cases, the microbunching is formed based on an electron bunch much longer 
than the radiation wavelength, for example in an FEL or CHG. In these cases, the 
linewidth of the longitudinal form factor at the harmonics are usually even narrower 
than that given by the transverse form factor. This also means that the coherent radi- 
ation of a long continuous electron bunch-based microbunching will be dominantly 
in the forward direction, as the bunching factor of the off-axis red-shifted frequency 
is suppressed very fast compared to the on-axis resonant ones. For a more prac- 
tical application, here we derive the coherent radiation power and spectral flux at 
the undulator radiation harmonics in these cases. As we will see, the results can be 
viewed as useful references for SSMB radiation. 

We assume that the long electron bunch, before microbunching, is Gaussain. The 
transverse form factors around the harmonics do not change much, i.e., we assume 


—4N, x S( H- & } . 
e i ( =) x 1 when is close to Ha. Therefore, we only need to take into 


account the Gaussian shape of the longitudinal form factors at the harmonics. The 
RMS bandwidth of the longitudinal form factor for a Gaussian bunch with a length of 


oz is Aw|, = —. Therefore, the coherent radiation energy at the H-th harmonic 
A k oz/ v2 


is 
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= (@ — Ho)? 
A (@ = Hag) x I. exp (eE) dw 


2e? 1+K?/2N, 
= NF FL (Sba n Č GE = 0) F ca x y 7c/0z. 
0 


(0) 


dW. 
Wy = [nerra ) F F; (0) —— 


(4.43) 
For a planar undulator, the o -mode radiation dominates and from Eq. (4.22) we have 


HK/V/2 i 
G0 = 0) = | ——"_—_|_ DJ}, 4.44 
(6 =0) kA [J (4.44) 
2 
in which the denotation JR = [Jsa (Hx)— Jan (Hx)] , with x = E is 


used. Note however, the above expression is meaningful only for an odd H, as the 
on-axis even harmonic radiation is rather weak. The peak power of the odd-H -th 
harmonic coherent radiation is then 


Wu 2 
Pit peak = ~—i— = ZN, Hx (JE, FF (Sb; yI, 
ipek = ee = go NXU IFF Obal (4.45) 


v2 


where Ip = MF < oe is the peak current of the Gaussian bunch before microbunching. 


For a more practical application of the derived formula, we put in the numerical value 
of the constants, and arrive at 


Pu peal KW] = 1.183N, H x (J J Jy F Fi (S)lb; n|" IBIA]. (4.46) 


Note that the above formula applies when the radiation slippage length N,,Ao is 
smaller than the bunch length o;. If not, the above formula will overestimate the 
coherent radiation peak power, as the RMS radiation pulse length is then longer than 


a Note also that given the same bunch charge and form factors, Pr peak % B x ae 


while Wy « —. The reason a shorter bunch radiates more total energy is baci 
that more particles are within the coherence length. 

At a first glance of Eq. (4.45), the coherent radiation power Pon seems to be 
proportional to N„, while an intuitive picture of the longitudinal coherence length 
loon X N, says that the scaling should be Pooh « N2 , as the electron number within 
the coherence length is proportional to N„. This is actually because that F F, (S) is 
also a function of N,. It is interesting to note that 


Pan = x N,, below eee ek (4.47) 
x N2, above diffraction limit, 


which can be obtained from the asymptotic expressions of FF, (S) as shown in 
Eq. (4.34). So for a given transverse beam size, Poon X N? at first when N, is small. 
When N, is large enough such that the electron beam is below diffraction limit, 
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then Pooh X N,. Physically this is because with the increase of N,,, the diffraction 
of the radiation will prevent the radiation from one particle so effectively affect the 
particles far in front of it, as the on-axis field from this particle becomes weaker with 
the increase of the radiation slippage length. 

Our derivation of the coherent radiation power above is for a Gaussian bunch- 
based microbunching. For a coasting or DC beam, we just need to replace Ip in 
Eq. (4.45) by the average current 74, and the peak power is then the average power. 
For a helical undulator, we need to replace K planar/ J2 with Kpeticas and [J J iF with 
1, in the evaluation of the radiation power at fundamental frequency. 

We remind the readers that Eq. (4.45) is for the case of a long continuous bunch- 
based microbunching, for example in FELs and CHG. In some of the SSMB sce- 
narios, the microbunches are cleanly separated from each other according to the 
modulation laser wavelength as will be shown in Fig. 4.9, and usually the radiation 
wavelength is at a high harmonic of the modulation laser. Therefore, there actu- 
ally could be significant red-shifted radiation generated in SSMB as we will see 
in the example calculation in Sect. 4.6. If we put the average current of SSMB in 
Eq. (4.45), what it evaluates is the radiation power whose frequency content is close 
to the on-axis harmonic and will underestimate the overall radiation power. 

After investigating the radiation power, let us now have a look at the spectral flux, 
which is the number of photons per unit time in a given small bandwidth. The spectral 


flux of coherent radiation at the odd-H-th harmonic can be calculated according to 
dWy 
as 


2 dWy Ao 
F (w = How) =| N; FF, (H, w)F F(a) x — | (w= Hao) 
dw ii ha 
1 e? 
= —— ——N,Hx[J J] FF (S)|b; p|? N? (phs/pass/0.1%b.w.), 
1000 Zesch XJ Jg FF (S)|b: u| N; (phs/pass ob.w.) 


(4.48) 
where fi is the reduced Planck’s constant. Again we put in the numerical value of the 
constants, and arrive at 


F (w = Hay) = 4.573 x 1079 N, Hx [J J Vy F FCS) bz, |? N2 (phs/pass/0.1%b.w.). 
(4.49) 


Note that the above spectral flux is for a single pass of the microbunched electron 
beam through the radiator undulator. For the evaluation of the average spectral flux 
in an SSMB storage ring, we need to multiply it with the number of microbunches 
passing a fixed location in one second, namely F rae with F being the filling factor 
of microbunches in the ring, v, being the average longitudinal speed of electron 
wiggling in the modulator undulator, and Àz being the modulation laser wavelength. 
We remind the readers that the above statement means we do not account for the 
radiation overlapping between different microbunches if the radiation slippage length 
is larger than =1_. If there is such radiation overlapping, the flux will be boosted 
further since the electrons in neighboring microbunches can now cohere with each 
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other. To give the readers a more concrete feeling about the high spectral flux in 
SSMB, we just need to multiply the spectral flux of the usual incoherent undulator 
radiation with a factor of N, F F, (S)|b; 4 |2, with N, being the number of electrons 
per microbunch. For example, in the envisioned EUV SSMB to be presented in 
Sect. 4.6, Ne = 2.2 x 104, and FF, (S)|b; H|? can be as large as 0.1. Therefore, the 
EUV spectral flux in an SSMB storage ring can thus be three orders of magnitude 
higher than that in a conventional synchrotron source. 


4.4 Impact of Electron Beam Divergence and Energy 
Spread 


Now we take into account the impact of beam divergence and energy spread on 
the coherent radiation. With an aim to obtain some efficient evaluation, here we 
simplify the analysis by considering only the impact of particle’s x’, y’ and 6 on the 
arrival time of the radiation, not on the radiation pattern. This approximation is valid 


when the beam divergence and energy spread are small enough, such that ox, y < 2 


and o5 < N: Basically we want to get a formula of the effective transverse and 
longitudinal form factors considering the beam size evolution during radiation. 

As an example, here we assume that the beam is a 6D Gaussian one, and round 
in the transverse dimension. Further we assume the beam reaches its minimal in all 
three dimensions at the radiator undulator center, which is desired to get high-power 


radiation, then the effective transverse and longitudinal form factors are 


L 
1 ft I+ le 
FEL E f FF (igs) A 


FR@)=— J È ol +(e razia) P (4.50) 


_ ea Vet (Peat) 
2 2o5N,ào | 
c 


where o1, 09,, 0, and o; are the transverse beam size, divergence, bunch length 
and energy spread at the undulator center, with oog, = €; and 0,05 = €,, and 
erf(x) = a N e`" dt is the error function. Note that VZ ef) < l forx £0. 

Figure 4.6 is an example plot of the effective form factors based on the above 
formulas, and the comparison with the case for a 3D rigid beam. It can be seen 
that given a transverse and longitudinal emittance, there is an optimal transverse 
beam size and bunch length at the radiator center considering the impact of beam 
divergence and energy spread. This is expected, since the beam size or bunch length 
of an over-focused beam grows very fast. 
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Fig. 4.6 An example plot to show the optimization of beam sizes in the middle of the radiator, 
considering the beam divergence and energy spread. Parameters used for calculation: €} = 1 nm, 
€: = 4pm, Ap = 13.5 nm, N, = 2 x 79, Ly = 1.58 m 


4.5 Statistical Property of Radiation 


In the previous sections, we have ignored the quantum discrete nature of radiation. 
Besides, we have derived the coherent radiation property using a smooth distributed 
charge, i.e., we have treated the charge as a continuum fluid. The number of photons 
radiated from a charged particle beam actually fluctuates from turn to turn or bunch 
to bunch if the quantum nature of radiation and the pointlike nature of electrons are 
taken into account. The first mechanism exists even if there is only one electron, and 
the second mechanism is related to the interference of fields radiated by different 
electrons [6]. Using the classical language, the second fluctuation mechanism is from 
the fluctuation of the bunching factor or form factor of electron beam. 

There have been studies on the statistical property of the radiation in FELs [7] and 
also the storage ring-based synchrotron radiation sources [6, 8]. Rich information 
about the electron beam is embedded in the radiation fluctuations, or more generally 
the statistical property of the radiation. For example, the turn-by-turn fluctuation of 
the incoherent undulator radiation can be used to measure the transverse emittance 
of electron beam [8]. The previous treatment, however, usually cares about the cases 
where the bunch length is much longer than the radiation wavelength, i.e., the radia- 
tion is temporally incoherent (in SASE FEL, incoherent at the beginning). In SSMB, 
the bunch length is comparable or shorter than the desired radiation wavelength, and 
the dominant radiation is temporally coherent. Although numerical calculation is 
doable following the general theoretical formulation, an analytical formula for the 
fluctuation in this temporally coherent radiation dominant regime is of value for a 
better understanding of the physics and investigation of its potential applications. 
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4.5.1 Pointlike Nature of Electron 


Here for SSMB we consider first the second mechanism of fluctuation, i.e., the 
radiation power fluctuation arising from the pointlike nature of the radiating electron. 
In this section, to simplify writing, we use the vector notation 

w 


k = — (sin 8 cos ọ, sin 8 sing, 1), 
c P á (4.51) 


r= (x,y,z): 


Then the bunching factor with the pointlike nature of electrons taken into account is 


N, 
1 - —ik-ry 

b(k) = b0, 9, œ) = cA 2 ge, (4.52) 

First we want to evaluate the coherent radiation power fluctuation at a specific 

frequency and observation angle. As the radiation power is proportional to N?|b(k)|?, 
therefore we need to know the fluctuation of |b(k)|?. Since 


Ne N, 
1 aA 1 I 
|b(k) |? = a Yee = aa | bet y eE (4.53) 


€ n=] m=1 € m#én 


we have 
(was?) = 5+ (1- =~) Bao? (4.54) 
Ne Ne , i 
with 
b(k) = b(0, 9, w) = I p(rye kT dr (4.55) 


being the bunching factor we calculated before using a continuum fluid charge dis- 
tribution in Eq. (4.3). As can be seen from Eq. (4.54), when Ne = 1, which corre- 
sponds to the case of a single point charge, then (|b(k)|*) = 1. When N, > 1 and 
N.|b(k)|? <_ 1, which corresponds to the case of incoherent radiation dominance, 
then (|b(k)|?) = ra When N, >> 1 and N,|b(k)|? >> 1, which corresponds to the 


case of coherent radiation dominance, then (|b(k) |?) = |b(k)|?. These results are as 
expected. 
The calculation of (|b(k1) 7 |b (kz) |?) is more involved. More specifically, 


Ne Ne Ne 


Ne 
l œ ; ; 
ko = ee ete la er), (4.56) 


€ n=1 m=1 p=1 q=1 
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Table 4.1 The N+ terms in the quadruple sum of Eq. (4.56) can be placed in 15 different classes, 


as shown in Ref. [9] 


Item number 


Index relations 


Number of terms 


(1) n=m=p=q Ne 

(2) n=m,p=q,n#p Ne(Ne — 1) 

(3) n=m,p#q#n Ne(Ne — 1)(Ne — 2) 
(4) n=p,m=q,n#Am Ne(Ne — 1) 

(5) n=p,m#q#n Ne(Ne — 1)(Ne — 2) 
(6) n=q,m=p,nAm Ne(Ne — 1) 

(7) n=q,mAp#nN Ne(Ne — 1)(Ne — 2) 
(8) n=m=p,n#q Ne(Ne — 1) 

(9) n=m=4q,n#p Ne(Ne — 1) 

(10) n=p=q,n#m Ne(Ne — 1) 

(11) p=q=mnAzm Ne(Ne — 1) 

(12) nxmAzApF| Ne(Ne — 1)(Ne — 2)(Ne — 3) 
(13) p=qnA#zmszp Ne(Ne — 1)(Ne — 2) 
(14) m=q,nA#Ams#zp Ne(Ne — 1)(Ne — 2) 
(15) m=pnAém#q Ne(Ne — 1)(Ne — 2) 


The NA terms in this summation can be placed in 15 different cases, as shown in 


Table 4.1. Corresponding to the 15 cases, we have 


(banta) = 


(4.57) 


N3 [Neb(ky — kı + k2 — k2) 

Ne(Ne — 1)b(kı — k1)b(k2 

Ne(Ne — 1)(Ne — 2)b(k1 — k1)b(k2)b(—k2) 
Ne(Ne — 1)b(ky + kz) b( 

Ne(Ne — 1)(Ne — 2)b(k1 + k2)b(—k1)b(—k2) 
Ne(Ne — 1)b(K1 — k2)b( 

Ne(Ne — 1)(Ne — 2)b(k1 — k2)b(—k1)b(k2) 
Ne(Ne — 1)b(Ki — kı + k2)b(—k2) 

Ne(Ne — 1)b(ky — kı — k2)b(k2) 

Ne(Ne — 1)b (ki + k2 — k2)b(—k1) 

Ne(Ne — 1)b(—ky + k2 

Ne(Ne — 1)(Ne — 2)(Ne 


3)b(k1)b(—k1)b(k2)b(—k2) 


tNe(Ne — 1)(Ne — 2)b(k1)b(—k1)b(k2 — k2) 
tNe(Ne — 1)(Ne — 2)b(k1)b(k2)b(—kı — k2) 


tNe(Ne — 1)(Ne — 2)b(k1)b(—k2)b(—ki + k2)] . 


The above result can be re-organized as 
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([b(k1)|?|b(k2) 7) = y [No + N2(Ne — 1) (|b(k1) I? + [bk2)1°) 
+Ne(Ne — 1)(Ne — 2) (b(ki + k2)b(—k1)b(—k2) + c.c.) 
+Ne(Ne — 1)(Ne — 2) (b(k1 — k2)b(—k1)b(k2) + c.c.) 
+Ne(Ne — 1) (\b(k1 + k2)|? + [bk — k2)1?) 


+N: (Ne — 1)(Ne — 2)(Ne — 3)|b(K1)|71b(k2) 17] , 


(4.58) 
in which c.c. means complex conjugate. 
If ky = kz = k, then 
2 1 2 
(1b) 1$) — (L091?) = 7g [Ne (Ne = I) + 2Ne(Ne = DN. — 2)15(8) |? 
— =2 
+2Ne(Ne — 1)(Ne — 2)Re [B2k)B (—k) (4.59) 


+N.(Ne — 1)|b(2k)|? 
—2N.(N. — 1)(2N. — 3)/B(K) I‘), 


where Re[] means taking the real part of a complex number. 
When N, > 1 and N.|b(k)|? <_1, which is the case for incoherent radiation 
dominance, we have (|b(k) |?) = + and 


(4.60) 


Var [|b(k)|7] 1+0( 1 ) 


(jnanry 


where Var[] means the variance of, and O(x”) means terms of order x” and higher. 
Therefore, the relative fluctuation of incoherent radiation is relatively large. This is 
also the reason why SASE-FEL radiation has a large shot-to-shot power fluctuation. 

When N, > 1 and N,|b(k)|? >> 1, which corresponds to the case of coherent 
radiation dominance like that in SSMB, we have 


Var [Ib 2] 2 (PWP + Re [be] dao ( i ) 
(bf Ne Ib(k)]4 
(4.61) 
The above equation is the main result of our analysis of bunching factor fluctuation 
for the regime of coherent radiation dominance, and to our knowledge is new. The 
formula can be used to evaluate coherent radiation power fluctuation at a specific 
frequency and observation angle. If the transverse electron beam size is zero, or if 
we observe on-axis, then we can just replace b(k) with b; (œ) in the above formula. 
We remind the readers that b(k) in general is a complex number. 
Now we conduct some numerical simulations to confirm our analysis of coherent 
radiation fluctuation. As can be seen from Fig. 4.7, which correspond to the cases of 
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Fig. 4.7 Fluctuation of the longitudinal form factor at 13.5 nm v.s. bunch length with Ne = 2.2 x 
10*. The bunch distribution is assumed to be Gaussian in the left and rectangular in the right, and 


the theoretical fluctuation is calculated according to Eq. (4.61), omitting the term O (az): For each 


parameters choice, 1 x 10+ simulations have been conducted to obtain the fluctuation 


a Gaussian and a rectangular distributed bunch, respectively, the simulation results 
agree well with our theoretical prediction. 

After investigating the expectation and variance of |b(k)|?, one may be curious 
about its more detailed distribution. It can be shown that when N,|b(k)|* >> 1, the 
distribution of |b(k)|? tends asymptotically toward Gaussian. 

As explained before, for a fixed frequency œw, there is a range of polar angle 0 
which can contribute. To evaluate the overall radiation power fluctuation at a specific 
frequency w, we then need to know the fluctuation of the form factor F F (w) which 
involves calculation depending on the specific radiation process. For our interested 
undulator radiation, it seems not easy to get a concise closed-form analytical formula 
to evaluate the total radiation power fluctuation when the beam has a finite transverse 
beam size. So here we refer to numerical calculation to give the readers a more 
concrete feeling about the impact of transverse size on coherent radiation power 
fluctuation. 

For simplicity, we assume that the bunch length is zero and focus on the fluctuation 
of the transverse form factor. As can be seen from the simulation result in Fig. 4.8, the 
larger the transverse beam size, the larger relative fluctuation of the transverse form 
factor. We also observe that in a typical parameters set of the envisioned EUV SSMB, 
the fluctuation of 13.5 nm radiation power due to the finite transverse size is small. 
For example ifo; = 16 um, then the relative fluctuation of the transverse form factor 
as shown in Fig. 4.8 is 0.14%. While the relative fluctuation of the longitudinal form 
factor at 13.5 nm when o, = 3 nm according to Eq. (4.61) is about 2%. Assuming 
that the beam is transverse-longitudinal decoupled, then 


Var[FF(w)] _ VarlFF,] | VarlFF.]  Var[F F, ]Var[ F F] 


(F F (@))? z (FFY (F F;)2 (FF |)2(FF,)2 (4.62) 
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Fig. 4.8 Fluctuation of the transverse form factor at 13.5 nm v.s. transverse beam size with Ne = 
2.2 x 104. The bunch is assumed to have zero length and is round Gaussian in the transverse plane. 
For each parameters choice, 1 x 10° simulations have been conducted to obtain the fluctuation. 
Parameters used for the calculation: Eọ = 400 MeV, A, = 1 cm, K = 1.14, N, = 2 x 79 


Therefore, for the envisioned EUV SSMB, the fluctuation of the longitudinal form 
factor dominates. 

After investigating the power fluctuation at a specific frequency w, now we look 
into the radiation energy fluctuation gathered within a finite frequency bandwidth 
and a finite angle acceptance. We use a filter function of FT (6, p, w) to account for 
the general case of frequency filter, angle acceptance, and detector efficiency. The 
expectation of the gathered photon energy and photon energy squared are 


T 20 foe) 2, 
(W) =n? | sinodo | ay f dwFT (6, 9, w) aW 0, 9, @) (|b, p, œ)1?) 
0 0 0 dwdQ point 
(4.63) 
and 
m 20 oo T 2n oo 
(W°) =N4 / sin 0d0 / dy f dw Í sin 0'd0' | dg’ / dw’ 
0 0 0 0 0 0 
aw aw 
FT 0, , FT 6’, h ! 6, s 6’, A f 
(0,9, @)FTO , p o) ag, g o| imdod Y a, 
(Jb, p, |?) (IDO, p, oP) 920, 8, p, p, w, 0), 
(4.64) 
where 


o, n (BG. p, DPO, g, o) 
BOO Pea a 


whose calculation can follow similar approach of calculating (|b(k1) |? |b (dk2) |?) in 
Eq. (4.57). And the relative fluctuation of the gathered photon energy is 


o3 = —1. (4.66) 
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4.5.2 Quantum Nature of Radiation 


As mentioned, there is another source of fluctuation, i.e., the quantum discrete nature 
of radiation. As a result of the Campbell’s theorem [10], we know that for a Poisson 
photon statistics, the variance of photon number arising from this equals its expec- 
tation value. With both contribution from pointlike nature of electrons and quantum 
nature of radiation taken into account, the relative fluctuation of the radiation power 
or energy at a given frequency and a specific observation angle is 


Var [P(o)] 1 2 (bw)? + Re [Eob —o)] n 1 
(P(@))? — (Nph())Ibeam Ne Dlo) | (az) 
(4.67) 
in which 


(Npn()) Ibeam = [Ne + Ne (Ne — 1)|6(@) 17] (Noh ()) Ipoint 


ays (4.68) 

~ Ne |b(@)|" (Npn()) Ipoint 
is the expected radiated photon number from the electron beam, and (Npn(@)) | point 
is the expected radiated photon number from a single electron. Note that to obtain 
a nonzero expected photon number (Npn(@))|beam, a finite frequency bandwidth is 
needed. Therefore, Eq. (4.67) actually applies to a finite frequency bandwidth close 
to w where b(w) does not change much. 

From Eq. (4.67), it is interesting to note that with the narrowing of the energy 
bandwidth acceptance, i.e., the decrease of (Npn(@)) |beam, the contribution to the 
relative fluctuation from the quantum nature of radiation increases, while the contri- 
bution from the pointlike nature of electron does not change. This reflects the fact 
that one fluctuation is quantum, while the other is classical. 

Note that in our interested case, Ne{Npn(@))|point is usually much larger than 1, 
then the second term in Eq. (4.67) dominants. In other words, the fluctuation due to 
the pointlike nature of electrons dominants. Only when Ne(Npn(@)) [point is Close to 
1, will the first term become significant compared to the second term. 

As the statistical property of the radiation embeds rich information about the 
electron beam, innovative beam diagnostics method can be envisioned by making 
use of this fact. One advantage of using radiation fluctuation in diagnostics is that 
it has a less stringent requirement on the calibration of detectors. Here we propose 
an experiment to measure the sub-ps bunch length accurately at a quasi-isochronous 
storage ring, for example the MLS, at a low beam current, by measuring and ana- 
lyzing the fluctuation of the coherent THz radiation generated from the electron 
bunch. Equation (4.67) or some numerical code based on the analysis presented in 
this section will be the theoretical basis for the experimental proposal. In principle, 
we can also deduce the transverse distribution of electron beam by measuring the 


112 4 SSMB Radiation 


two-dimensional distribution of the radiation fluctuation. More novel beam diagnos- 
tics methods may be invented for SSMB and future light sources by making use of 
the statistical property of radiation. 


4.6 Example Calculations for Envisioned EUV SSMB 


To summarize our investigations on the average and statistical property of SSMB 
radiation, here we present an example calculation for the envisioned EUV SSMB. 
Figure 4.9 is an example plot of the beam current and longitudinal form factor spectra 
of the envisioned EUV SSMB. In the envisioned example, the microbunch length 
is o, © 3 nm at the radiator where 13.5 nm coherent EUV radiation is generated, 
and these 3 nm microbunches are separated from each other with a distance of 


150 T T T 


Fig. 4.9 An example plot of the beam current and longitudinal form factor spectrum of the 
microbunch train at the radiator in the envisioned EUV SSMB. Up: beam current of the 3 nm 
microbunch train separated by the modulation laser wavelength àz = 1064 nm. Bottom: longitu- 
dinal form factor F F; (œw). The exponential decaying envelope corresponds to that of a single 3 
nm Gaussian microbunch, and the green periodic delta function lines correspond to the periodic 
microbunch train in time domain. The desired radiation wavelength is Ag = aL = 13.5 nm 
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àz = 1064 nm = 79 x 13.5 nm, which is the modulation laser wavelength. The 
radiator is assumed to be an undulator. The beam at the radiator can be round or flat 
depending on the lattice scheme, and its transverse size can range from a couple of 
um to a couple of 10 um. Note that a Gaussian-distributed current at the radiator 
is assumed in the plot. This is the case corresponding the usual longitudinal strong 
focusing SSMB. We remind the readers that the current at the radiator in the TLC- 
based bunch compression scheme, i.e., generalized longitudinal strong focusing, is 
actually non-Gaussian considering the nonlinear modulation waveform, as shown in 
Fig. 3.4. As our goal is to give the readers a picture of the radiation characteristics, 
here for simplicity we consider the case of a round beam. We remind the readers 
that the parameters used in this example EUV SSMB radiation calculation is for an 
illustration and is not optimized. 


4.6.1 Average Property 


First we present the result for the average property of the EUV radiation. The calcu- 
lation is based on Eqs. (4.7), (4.11) and (4.22), and the result is shown in Fig. 4.10. 
The upper part of the figure shows the radiation energy spectrum. The lower part 
shows the spatial distributions of the radiation energy. The total radiation power is 
calculated according to 

W 


= e (4.69) 
where W is the total radiation energy loss of each microbunch. For the example radi- 
ator undulator parameters, corresponding to 0, = 5, 10, 20 um, the total radiation 
power are 92 kW, 14 kW, 3.5 kW, respectively. As a reference, the radiation power 
calculated based on Eq. (4.45) for these three transverse beam sizes are 4 kW, 3.6 
kW and 2.6 kW, respectively. The reason Eq. (4.45) gives a smaller value than the 
overall power as explained is that it does not take into account the red-shifted part 
of the radiation. Therefore, Eq. (4.45) can be used to evaluate the lower bound of 
the radiation power from SSMB, once the parameters set of electron beam and radi- 
ator undulator is given. It can be seen that generally, kW-level EUV radiation power 
can be straightforwardly anticipated from a 3 nm microbunch train with an average 
beam current of 1 A. Note that for simplicity, in this example calculation, the filling 
factor of microbunches in the ring is assumed to be 100%, i.e., one microbunch per 
modulation laser wavelength. Then 1 A average current corresponds to the num- 
ber of electrons per microbunch N, = Tahe = 2.2 x 10*, if the modulation laser 
wavelength is àz = 1064 nm. 

Another important observation is that the spectral and spatial distribution of SSMB 
radiation depends strongly on the transverse size of the electron beam. A large trans- 
verse size results in a decrease of the overall radiation power, and also makes the 
radiation more narrowbanded and collimated in the forward direction. This is an 
important observation drawn from our investigation on the generalized transverse 
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Fig. 4.10 An example EUV SSMB radiation calculation with a microbunch length of o, = 3 nm 
and different transverse sizes o}. The upper part shows the energy spectrum. Corresponding to 
o1 = 5, 10, 20 um, the total radiation power are 92 kW, 14 kW, 3.5 kW, respectively. The shaded 
area corresponds to wavelength of 13.5 + ae nm. The bottom part shows spatial distribution of 
radiation energy. From left to right: o1 = 5, 10, 20 um. Parameters used for the calculation: Eg = 
400 MeV, Iavg = 1 A, A, = 1064 nm, A; = ue = 13.5 nm, A, = 1 cm, K = 1.14, N, =2 x 79 


form factor. Using the example parameters, i.e., Eg = 400 MeV, A; = 1064 nm, 
do = ue = 13.5 nm, A, = 1 cm, N, = 2 x 79, K = 1.14, if o, = 10 wm, then 
according to Eqs. (4.39) and (4.40), the relative bandwidth and opening angle due to 
the transverse form factor can be calculated to be 


Awe- 
De I: 
w0 
md (4.70) 
0e-| = 0.21 mrad. 


which is in agreement with the result presented in Fig. 4.10. 

Since N, = 2 x 79 is used in the calculation, which means the radiation slippage 
length is twice the modulation laser wavelength, the energy spectrum and spatial 
distributions presented in Fig. 4.10 is for that of two neighboring microbunches. 
The peaks in the energy spectrum, and the interference rings in the radiation spatial 
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si 
distribution as explained is due to the macro longitudinal form factor ( 


of multiple microbunches, with N, the number of microbunches. There is a one-to- 
one correspondence between the peaks in the energy spectrum and the interference 
rings in radiation energy spatial distribution. The reason for the appearance of such 
peaks and rings is that our radiation wavelength is a high harmonic of the distance 
between the neighboring microbunches. 

As a result of the high-power and narrowband feature of the SSMB radiation, 
a high EUV photon flux of 10!5 ~ 10!° phs/s within a 0.1 meV bandwidth can be 
obtained, if we can realize an EUV power of > 1 kW per 1% b.w. as shown in Fig. 4.10. 
We remind the readers that the radiation waveform of SSMB is actually a CW or 
quasi-CW one, if induction linac is used as the energy compensation system and the 
microbunches occupy the ring with a large filling factor, as assumed in the example 
calculation. This kind of CW or quasi-CW narrowband photon source is favored in 
ARPES to minimize the space charge-induced energy shift, spectral broadening and 
distortion of photoelectrons in a pulsed photon source-based ARPES [11]. Therefore, 
the high photon flux within a narrow bandwith, together with its CW or quasi-CW 
waveform, makes SSMB a promising light source for ultrahigh-resolution ARPES. 
Such a powerful tool may have profound impact on fundamental physics research, for 
example to probe the energy gap distribution and electronic states of superconducting 
materials like the magic-angle graphene [12]. 


4.6.2 Statistical Property 


Now we present the result for the statistical property of the radiation. For the case of a 
Gaussian bunch with o, = 3 nm and N, = 2.2 x 104, from Eq. (4.61) we know that 
the relative fluctuation of the turn-by-turn or microbunch-by-microbunch on-axis 
13.5 nm coherent radiation power will be around 2%. 

Figure 4.11 gives an example plot for the longitudinal form factor spectrum of 
three possible realizations of such a Gaussian microbunch. As can be seen, the spec- 
trum is noisy mainly in the high-frequency or short-wavelength range. Since our 
EUV radiation is mainly at the wavelength close to 13.5 nm, and the longitudinal 
form factor close to this frequency fluctuates together from turn to turn, or bunch to 
bunch. As shown in Fig. 4.8 and discussed before, for the envisioned EUV SSMB, 
the transverse form factor fluctuation is much smaller than that of the longitudinal 
form factor. So the overall radiation power fluctuation is also about 2% as analyzed 
above. This fluctuation is also the fluctuation of microbunch center motion induced 
by the coherent radiation. A small fluctuation as it is, its beam dynamics effects need 
further study. 

Note that this 2% fluctuation of radiation power should have negligible impact for 
the application in EUV lithography, since the revolution frequency of the microbunch 
in the ring is rather high (MHz), let alone if we consider that there is actually a 
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microbunch each modulation laser wavelength and the radiation waveform is CW or 
quasi-CW. 


4.6.3 Discussions 


To resolve possible concerns of readers on the validity of the short bunch length and 
high average current used in the example calculation, here we present a short dis- 
cussion on the related single-particle and collective effects in SSMB. We recognize 
that realizing a steady-state bunch length as short as nanometer level in an electron 
storage ring is non-trivial. Both global and local momentum compaction should be 
minimized to confine the longitudinal beta function, therefore the longitudinal emit- 
tance, in an electron storage ring as analyzed in Sect. 2.1. By invoking this principle 
in the lattice design, a bunch length as short as tens of nanometer can be realized in a 
storage ring, with a momentum compaction factor of 1 x 10~° and modulation laser 
power of 1 MW. 1 MW intra-cavity power is the state-of-art level of present opti- 
cal enhancement cavity technology. Therefore, to realize nanometer bunch length 
at the radiator, we need to compress the bunch further. There are two scenarios 
being actively studied by us, namely the longitudinal strong focusing scheme and 
transverse-longitudinal coupling scheme. The longitudinal strong focusing scheme 
is similar to its transverse counterpart which is the basis of modern particle accel- 
erators [13, 14]. In such a scheme, the longitudinal beta function and therefore the 
bunch length is strongly focused at the radiator, and the synchrotron tune of the beam 
in the ring can be at the level of 1, as analyzed in Sect. 2.1.6 and Ref. [15]. Although 
nanometer bunch length can be realized, this scheme requires a large modulation 
laser power (20 MW level if 270 nm laser is used), causing the optical cavity can 
work only in the pulsed laser mode and the average output radiation power is thus 
limited. To lower the modulation laser power, the transverse-longitudinal coupling 
scheme is thus applied in a clever way by taking advantage of the fact that the ver- 
tical emittance in a planar storage ring is rather small, as analyzed in Sect. 3.1. We 
refer to this turn-by-turn transverse-longitudinal coupling-based bunch compression 
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scheme as the generalized longitudinal strong focusing [16], in which the phase space 
manipulation is 4D or 6D, in contrast to the conventional longitudinal strong focus- 
ing where the phase space manipulation is 2D. This generalized longitudinal strong 
focusing scheme can relax the modulation laser power, but its nonlinear dynamics 
optimization is a challenging task which we are trying to tackle. 

Concerning the high average current, there are two collective effects of special 
importance, namely the intrabeam scattering (IBS) and coherent synchrotron radia- 
tion (CSR). IBS will affect the equilibrium emittance and thus can have an impact 
on the radiation power and also the modulation laser power in the generalized longi- 
tudinal strong focusing scheme. The IBS effect in an SSMB ring thus needs careful 
optimization and the operation beam energy is also mainly determined by IBS. CSR 
is the reason why SSMB can provide powerful radiation. On the other hand, CSR 
is also the effect which sets the upper limit of the stable beam current. In Ref. [17], 
there is some preliminary evaluation of the stable beam current for SSMB based on 
the 1D model of CSR-driven microwave instability. The investigation in this dis- 
sertation implies that the transverse dimension of the electron beam can have large 
impact on the coherent radiation in SSMB. In addition, the bunch lengthening from 
the transverse emittance in an SSMB storage ring can easily dominate the bunch 
length at many dispersive places of the ring, as the transverse size of microbunches 
is much larger than its longitudinal length. This bunch lengthening might be helpful 
in mitigating unwanted CSR. The 3D effect of the coherent radiation is expected to 
be also helpful in improving the stable beam current. With these beneficial arguments 
in mind, we realize that CSR in SSMB still deserves special attention. For example, 
the coherent radiation in the laser modulator could potentially drive single-pass and 
multi-pass collective instabilities in an SSMB storage ring [18, 19]. More in-depth 
study of collective effects in SSMB is ongoing. 
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Chapter 5 A) 
SSMB Proof-of-Principle Experiments geai 


To make SSMB a real option for future photon source facility, a crucial step is to 
experimentally demonstrate its working mechanism. In this chapter, we report the 
first successful experimental demonstration of the SSMB mechanism. Parts of the 
work presented have been published in Refs. [1, 2]. 


5.1 Strategy of the PoP Experiments 


5.1.1 Three Stages of PoP Experiments 


Considering the fact that it is a demanding task to realize SSMB directly in an existing 
machine, part of the reasons we have analyzed in previous chapters, among them the 
most fundamental one is the large quantum diffusion of bunch length in rings not 
optimized for SSMB, the SSMB PoP experiment has been divided into three stages 
as shown in Fig. 5.1. Some brief descriptions of the three stages are as follows. 


e Phase I: a single-shot laser is fired to interact at the undulator with the electron 
beam stored in a quasi-isochronous ring. The modulated electron beam becomes 
microbunched at the same place of modulation after one complete revolution in 
the ring and this microbunching can preserve for several revolutions. By doing 
this experiment, we want to confirm that the optical phases, i.e., the longitudinal 
coordinates, of electrons can be correlated turn-by-turn in a sub-laser-wavelength 
precision. The realization of SSMB relies on this precise turn-by-turn phase cor- 
relation. 

e Phase II: on the basis of Phase I, we replace the single-shot laser with a high- 
repetition phase-locked one to interact with the electrons at the undulator turn 
after turn. In this stage, we want to establish stable microbuckets and sustain the 
microbunching in the microbuckets to reach a quasi steady state. 


© The Author(s) 2024 119 
X. Deng, Theoretical and Experimental Studies on Steady-State Microbunching, 
Springer Theses, https://doi.org/10.1007/978-98 1-99-5800-9_5 


120 5 SSMB Proof-of-Principle Experiments 


Undulator 


Quasi-isochronous storage ring 


e~ beam 
Phase |: a single Phase II: multiple laser Phase Ill: infinite laser 
laser shot, shot-lived shots, quasi-steady-state shots, real-steady-state 
microbunching microbunching microbunching 


: : pen X A 
VV) 000 OO 


Fig. 5.1 Three stages of the SSMB PoP experiments: from single-shot to multiple shots to infinite 
shots laser pulse; from short-lived to quasi-steady-state to real steady-state microbunching 


e Phase III: Phase II is very close to, but is still not, the final SSMB as a true SSMB 
means the balance of excitation and damping. However, the requirement of a true 
SSMB on the magnet lattice is demanding, especially the quantum diffusion of 
longitudinal coordinate z as analyzed in Sect. 2.1. Therefore, this final stage is 
more likely to be realized in a dedicated ring designed for SSMB, which is also 
one of the key ongoing tasks of the SSMB task force [3]. 


Below, we use PoP I, II, II to represent the three stages of the experiment. The key 
words of the three stage experiments are summarized as follows. 


e PoP I: microbunching based on stored electron bunch, turn-by-turn phase correla- 
tion; 

e PoP II: bounded motion in microbuckets, quasi steady state; 

e PoP III: balance of diffusion and damping, real steady state. 


These three stages each have their own significance and are all important for the 
SSMB development. Among them, Phase I is from 0 to 1, and is the most important 
one from conceptual viewpoint. Recently, we have successfully performed the PoP 
I and demonstrated the mechanism of SSMB at the Metrology Light Source (MLS) 
of Physikalisch-Technische Bundesanstalt (PTB) in Berlin [1, 2]. The experiment is 
a collaboration work of Tsinghua, Helmholtz-Zentrum Berlin (HZB) and PTB. 


5.1 Strategy of the PoP Experiments 


Table 5.1 Basic parameters of the MLS lattice 
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Parameter Value Description 

Eo 50-630 MeV Beam energy 

Co 48 m Ring circumference 

ÍRE 500 MHz RF frequency 

VRF <600 kV RF voltage 

n 3 x 107? Phase slippage factor (standard user) 

n —2 x 1075 Phase slippage factor (SSMB experiment) 

Uo 226 eV @250 MeV Radiation loss per turn 

Js 1.95 Longitudinal damping partition 

T 180 ms@250 MeV Longitudinal radiation damping time 

05 1.8 x 10-4@250 MeV | Natural energy spread 

Oz 36 um (120 fs)@250 Zero-current bunch length (SSMB 
MeV experiment) 

Vy 3.18 Horizontal betatron tune 

Vy 2.23 Vertical betatron tune 

Ex 31 nm@250 MeV Horizontal emittance 

ÀL 1064 nm Modulation laser wavelength 

Àu 125 mm Undulator period length 

Nu 32 Number of undulator periods 

Lu 4m Undulator length 

K 2.5 Undulator parameter 


5.1.2 Metrology Light Source Storage Ring 


The MLS is a storage ring optimized for quasi-isochronous operation [4—6], thus an 
appropriate testbed for SSMB physics investigation and PoP experiments. However, 
the partial phase slippage of the MLS is large as the bending angle of each dipole 
is large (7) and the dispersion magnitude inside the dipoles is also large, so it is 
not feasible to realize true SSMB, i.e., PoP III, directly at the MLS. Therefore, the 
SSMB PoP experiment has been divided into three stages as introduced just now, 
and PoP I and II are what we have performed and plan to conduct at the MLS. Some 
basic parameters of the MLS are shown in Table 5.1. The lattice optics of the MLS 
used in the SSMB PoP experiments are shown in Fig. 5.2. 
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Fig. 5.2 The MLS quasi-isochronous magnet lattice used to generate microbunching. The magnet 
lattice and the key are shown at the top. The curves are the model horizontal (red) and vertical 
(blue) -functions and the horizontal dispersion D, (green). Operating parameters of the ring: beam 
energy, Eo = 250 MeV; relative energy spread, os = 1.8 x 1074 (model); horizontal emittance, €x 
= 31 nm (model); horizontal betatron tune, vy = 3.18 (model and measured); vertical betatron tune, 
vy = 2.23 (model and measured); horizontal chromaticity, £y = —0.5 (measured). Note that this 
optics is different from that used in Sect. 2.1.3 for the simulation of partial phase slippage effect. 
(Figure from Ref. [1]) 


5.2 PoP I: Turn-by-Turn Laser-Electron Phase Correlation 


5.2.1 Experimental Setup 


Figure 5.3 shows the schematic setup of the SSMB PoP I experiment. A horizontally 
polarized laser pulse (wavelength, àz = 1064 nm; pulse length, full-width at half- 
maximum, FWHM ~* 10 ns; pulse energy, ~ 50 mJ) is sent into a planar undulator 
(period, à, = 125 mm; total length, L, = 4 m) to co-propagate with the electron 
bunches (energy, Eo = 250 MeV) stored in the MLS storage ring (circumference, Co 
= 48 m). To maximize the laser-electron energy exchange, the undulator gap is cho- 
sen to satisfy the resonance condition Às = Az, where A; = HE ai is the central 
wavelength of the spontaneous undulator radiation, with y « Ey being the Lorentz 
factor and K = a = 0.934 - Bo[T] - A,,[cm] being the dimensionless undulator 
parameter, determined by the undulator period and magnetic flux density. This laser- 
electron interaction induces a sinusoidal energy modulation pattern on the electron 
beam with a period of the laser wavelength. Because particles with different ener- 
gies have slightly different revolution periods, after one revolution in the ring, the 
energy-modulated electrons shift longitudinally with respect to each other, clumping 
towards synchronous phases and forming microbunches. The formed microbunches 
can last several revolutions in the ring. The coherent undulator radiation generated 
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Fig. 5.3 Schematic of the experimental setup. The stored 250 MeV electron bunches are energy- 
modulated by a 1064 nm wavelength laser in an undulator, and become microbunched after one com- 
plete revolution in the 48 m circumference quasi-isochronous storage ring. This formed microbunch- 
ing can then preserve for multiple turns in the ring. Each time the microbunching going through the 
undulator, narrowband coherent radiation will be generated. The undulator radiation is separated 
into the fundamental and second harmonics by dichroic mirrors, and sensitive photodiodes are used 
as the detectors. Narrow band-pass filters can be inserted in front of the photodetectors to pick out 
the narrowband coherent radiation generated from the microbunching. (Figure from Ref. [1]) 


from the microbunches, detected by a high-speed photodetector with a photodiode, 
confirms microbunching. 

The symplectic longitudinal dynamics of the the above experiment processes can 
be modeled by 


ôi = ĉo + Asin(kz Zo), 


(5.1) 
Z1 = Z0 — nCod1, 
for the first revolution with laser modulation, and 
Ôm = Ôm ; 
H (5.2) 


Zm+1 = Zm — NCodm-+1; 
for the later revolutions. This demonstration proves that the longitudinal dynamics 
described by the one-turn map Eq. (5.1) can be extrapolated from the RF wavelength 


(metre scale) to laser wavelength (micrometre scale) for a stored electron beam, thus 
validating the SSMB microbunching mechanism. 


5.2.2 Physical Analysis of Microbunching Formation 


5.2.2.1 Storage Ring 


Operation energy The above models Eqs. (5.1) and (5.2), however, do not consider 
the non-symplectic, transverse-longitudinal coupling and nonlinear lattice dynamics, 
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which all could lead to degradation of the microbunching. It turns out that the first 
non-symplectic dynamics we need to account for is the synchrotron radiation. As 
we know, when a relativistic electron is subjected to an acceleration normal to its 
velocity exerted by a bending magnet, it radiates electromagnetic energy. This radi- 
ation is characterized by the quantum nature of the photon emission process. The 
photon energy and emission place or time are both stochastic, giving rise to changes 
on particle energy (instantly) and the longitudinal coordinate z (non-instantly), as 
analyzed before in Sect. 2.1. Of special interest in the experiment is the root-mean- 
square (RMS) quantum diffusion of z in one turn d,. According to Eq. (2.41), we 
have 


2. 
d; = (22) — R} = | (F2)(N) (a). (5.3) 


with F (s2, s1) = —N(s2, s1)Co. For the MLS quasi-isochronous magnet lattice used 
in the PoP experiment as shown in Fig. 5.2, d, is as large as 260 nm at its standard 
operation energy of 630 MeV, deteriorating the sub-micrometre microbunching con- 
siderably. Therefore, the beam energy needs to be lowered to mitigate this diffusion, 


as ,/ (N) ( “| œ y*>. At the same time, a lower beam energy gives a smaller energy 


spread and is also beneficial for microbunching, as the smearing from the natural 
uncorrelated energy spread becomes smaller. Nevertheless, the beam energy cannot 
be too low, otherwise the beam parameters and lifetime could be profoundly affected 
by scattering among particles [7—9]. An electron beam energy of 250 MeV is adopted 
in the experiment to balance these issues. At Eg = 250 MeV, we have d, = 26 nm 
and o; = 1.8 x 1074. 


Phase slippage factor Because the laser wavelength is much smaller than that of 
an RF wave, the phase slippage factor 7 needs to be ultrasmall. That is, the ring 
should be quasi-isochronous to allow turn-by-turn stabilization of the electron optical 
phases, i.e., the longitudinal coordinates, for particles with different energies. More 
quantitatively, the RMS spread of z in one turn that arises from the uncorrelated 
electron energy spread should be adequately smaller than the laser wavelength, 


Azes = |nCoo5| < 41/27. (5.4) 


To fulfill this requirement, the phase slippage factor of the MLS was lowered to 
nx —2 x 1075, which is three orders of magnitude smaller than its standard value 
of 3 x 107”. By implementing these parameters, Azgs = |7Coo3| © 0.17 pm (about 
0.6 fs), enabling the formation and preservation of sub-micrometre microbunching. 

Such a quasi-isochronous magnet lattice is achieved by tailoring the horizontal 
dispersion functions D, around the ring so that a particle with non-ideal energy trav- 
els part of the ring inwards and part of the ring outwards compared to the reference 
orbit, thus having a revolution period nearly the same as that of the ideal particle. 
The tailoring of Dx is accomplished by adjusting the (de)focusing strengths of the 
quadrupole magnet. The operation of the MLS as a quasi-isochronous ring also 
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benefits from the optimization of the sextupole and octupole nonlinear magnet 
schemes to control the higher-order terms of the phase slippage [5, 6], which affect 
both the equilibrium beam distribution in the longitudinal phase space before the laser 
modulation and the succeeding microbunching evolution as analyzed in Sect. 2.2.1. 

In the experiment, the value of the small phase slippage factor is quantified by 
measuring the electron orbit offsets while slightly adjusting the RF frequency up and 
down at a beam position monitor (BPM), where D, is large. From the offsets and 
D,, the phase slippage-dependent electron energy shifts caused by the RF frequency 
adjustment can be derived, as well as the phase slippage factor. The D, value at 
the BPM is acquired using the same method in a reversible way, that is, based on a 
known phase slippage factor. This is done at a larger phase slippage factor, at which 
its value can be determined from its relation to the synchrotron oscillation frequency 
of the electron beam as given in Eq. (2.28). The synchrotron oscillation frequency 
can be measured accurately at a phase slippage factor such as —5 x 1074, and the 
model confirmed that the relative change of D, at the highly dispersive BPM is small 
(< 4%) when the phase slippage factor is reduced from —5 x 1074 to the desired 
—2 x 107 by marginally changing the quadrupole magnet strengths. 


Bunching factor With the electron beam evolved according to Eq. (5.1) for one 
revolution in the ring, the bunching factor at the n-th laser harmonic as analyzed in 
Sect. 2.2.1 is 


(5.5) 


kinCoos)2 
bn = Ja (nkr nCoA)exp [0p] l 


2 


where J,, is the n-th order Bessel function of the first kind. The coherent radiation 
power at the n-th harmonic is proportional to the bunching factor squared, P,, coh X 
|b, |?. With the dynamics in the following turns modeled by Eq. (5.2), for the m-th 
revolution, we just need to replace the nCo in Eq. (5.5) with mnCo. 

The maximum reachable bunching factor becomes larger with the decrease of n, 
given that the optimal A can always be realized. n = —2 x 1075 is approximately 
the present lowest reachable value at the MLS, so here below we use this 7 for the 
analysis. As we will explain soon, our signal detection at first focuses on the sec- 
ond harmonic of the undulator radiation, since it is easier to mitigate the impact of 
the modulation laser on signal detection compared to the fundamental frequency. 
As shown in Fig. 5.4, given n = —2 x 1075, a modulation strength of A ~ 1.505 
results in the maximum bunching at the second harmonic. Correspondingly, given 
A = 1.50;3, the optimal n for the fundamental frequency and second harmonic bunch- 
ing is a bit smaller than 2 x 1075. Note that the optimized conditions for the fun- 
damental frequency and second harmonic bunching are different. Figure 5.5 shows 
the bunching factor evolution with respect to the revolution number, with A = 1.505 
and |n| = 2 x 1075. As can be seen, the second harmonic bunching can last only one 
turn, while the fundamental frequency bunching can last about three turns or even 
more if 7 becomes smaller, as also can be seen in the right part of Fig. 5.4. These 
expectations have also been confirmed in the experiment as will be presented soon. 
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Fig. 5.4 Left: impact of the energy modulation strength (A « ./ PL) on the bunching factor |b, | at 
the fundamental and second harmonic, with n = —2 x 1075. Right: impact of the phase slippage 
factor ņ on the bunching factor |b, | at the fundamental and second harmonic, with A = 1.505 


Fig. 5.5 The evolution of 0.4 
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Chromatic H, function and chromaticity £. Apart from the longitudinal beam 
dynamics, the coupling of the particle betatron oscillation to the longitudinal dimen- 
sion is also critical. The reason is based on the fact that the horizontal beam width 
at the undulator is about 600 um (model value), three orders of magnitude larger 
than the sub-micrometre longitudinal structures that we aim to produce. Because the 
vertical emittance is much smaller than the horizontal one in a planar x-y uncoupled 
storage ring, in the following we consider only the impact of the horizontal betatron 
oscillation. 

According to Eq. (3.13), for a periodic system, the RMS bunch lengthening of 
an electron beam longitudinal slice after m complete revolutions in the ring, due to 
betatron oscillation, is 


Azam = 2V eH, |sin(mzv,)| . (5.6) 
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Fig. 5.6 Influence of H, at 
the undulator on the 
bunching factor one turn 
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With this bunch lengthening taken into account, the bunching factor at the n-th laser 
harmonic after m revolutions will be 


nkimnCoos)? + (nkt: V&H, |sin(m7 vy 
ban = nko Ades | (nk,mnCoos)? + ( - « |sin( 1) l 


(5.7) 
The relative bunching factor reduction due to the non-zero H, at the undulator can 
thus be defined as 


(5.8) 


2 


Ram (Hy) = exp l Can renan] 


Putting in €, =31 nm (model) and v, = 3.18 (model and measured), we need H, < 
0.8 um at the undulator to have Azg,ı < àz/2x. Figure 5.6 shows the bunching 
factor reduction at the fundamental frequency and the second harmonic one turn 
after the laser modulation as a function of the H, at the undulator. As can be seen, 
the second-harmonic bunching is even more sensitive to the H, at the undulator. 
This stringent condition on H, (note that Hy at other places of the ring is typically 
>0.1 m) is satisfied by fine-tuning the quadrupole magnet (de)focusing strengths to 
correct the dispersion D, and dispersion angle D’. at the undulator to the level of 
millimetre and 0.1 mrad, respectively (see Fig. 5.2). 

In addition to the linear-order oscillating bunch lengthening, as discussed in 
Sect. 3.2, the betatron oscillation also produces an average path lengthening or short- 
ening (second-order effect) described by the formula 


ACg = —27 J Ex, (5.9) 


with ACg being the average change of the particle recirculation path length, and 
&, = dv,/dé being the horizontal chromaticity of the ring. Because different parti- 
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cles have different betatron oscillation amplitudes (actions), this effect results in a 
loss of synchronization between particles and degrades microbunching. Moreover, 
it broadens the equilibrium energy spread and distorts the beam from the Gaussian 
form before the laser modulation as investigated in Sect. 3.2, which also affects the 
microbunching. Therefore, the horizontal chromaticity should be small, to moderate 
its detrimental outcome, and simultaneously sufficient to suppress collective effects 
such as the head-tail instability [10]. As a consequence, a small negative chromaticity 
is used in the experiment. 


5.2.2.2 Modulation Laser 


Long-pulse laser A long-pulse laser (FWHM ~ 10 ns) has been used to simplify 
the experiment by avoiding a dedicated laser-electron synchronization system, as 
the shot-to-shot laser timing jitter is fjitter < 1 ns (RMS). According to Eq. (5.5) and 
Fig. 5.4, for a given phase slippage factor and harmonic number, there is an optimal 
laser-induced energy modulation amplitude A (A œ ./P;, with Pz the laser power) 
that gives the maximum bunching factor. The laser used in the experiment (Beamtech 
Optronics Dawa-200) has multiple longitudinal modes, and its temporal profile has 
several peaks and fluctuates considerably from shot to shot (see Fig. 5.7). Therefore, 
the laser-induced electron energy modulation amplitudes are different from shot to 
shot and from bunch to bunch. When the modulation amplitude matches the phase 
slippage factor, the energy-modulated electrons are properly focused at synchronous 
phases, which gives optimal microbunching. For some of the shots, the laser intensity 
is higher or lower than the optimal value, and the electrons are then over-focused or 
under-focused, giving weaker microbunching and less coherent radiation. As we will 
see soon, this explains the large shot-to-shot fluctuation of the coherent amplified 
signals shown in Fig. 5.16c, e. 


Power and Rayleigh length The electric field of a TEM00 mode Gaussian laser 
beam is [11] 


Pa ease 1 k 
Ey = E,gel2-ie'+ido —exp E LO (2 + | ; 


1+i Ze 2 (5.10) 
E, ~% -—E,x, 
with Zr = uaa the Rayleigh length, wo the beam waist radius, and Q = EET 
Z 
The relation between Exo and the laser peak power is given by i 
P, = ExoZeke (5.11) 
ee Agee 


in which Zp = 376.73 Q is the impedance of free space. The electron wiggles in an 
undulator according to 
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Fig. 5.7 Fluctuating temporal profiles of the multilongitudinal-mode laser. a, Temporal profiles of 
two example consecutive laser shots (red and blue) and the averaged waveform of 200 consecutive 
laser shots (black). b, Statistical distribution of the laser power at t = 0 ns in a for 10* consecutive 
laser shots, where the red curve is a gamma distribution fit. Laser: compact Nd:YAG Q-switched 
laser (Beamtech Optronics Dawa-200). Detector: ultrafast photodetectors (Alphas UPS-40-UVIR- 
D; rise time < 40 ps). Measurement system: digital oscilloscope (Teledyne LeCroy WM825Zi-B; 
bandwidth 25 GHz; sample rate 80 billion samples per second). (Figure from Ref. [1]) 


K . 
x(z) = — sin(k,z), (5.12) 
yk 


u 
and the laser-electron exchange energy according to 


a = ev, E, + ev, E. (5.13) 
Assuming that the laser beam waist is in the middle of the undulator, and when 
x, y < w(z), which is the case for SSMB PoP I, we drop the exp [5e (x? + y’)] 
in the laser electric field. Further, when Zr >> 4,, we can also drop the contribution 
from E; on the energy modulation. The integrated modulation voltage induced by 
the laser in the planar undulator, normalized by the electron beam energy, is then 


[11] 
e Vinod _ = 4P,ZoZr tan-! Ly (5.14) 
Eo y mce ÀL 2ZR 


in which [JJ] = Jo(x) — J1 (xX) and x = zÉ- 
As can be seen from the above formula, when + is kept constant, then Vmod & 
VZr X 4 Lu. Inour cases L,, is fixed, then as res in Fig. 5.8, to maximize Vinod We 


need Ze = = 0. 359 ~ 3. On the other hand, the modulation strength does not depend 
on ZR Seine when Zp is larger than this optimal value. To make the laser beam 
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waist larger than the electron beam and thus induce the same energy modulation 
on different electrons, a larger Rayleigh length might be in favored in our case. 
For example if Zr = 2L,, then in order to induce an energy modulation depth of 
A = 1.505, we can calculate that the laser power required is 430 kW. Considering the 
non-ideal conditions and the fact that the laser used contains higher-order Gaussian 
modes, one order of magnitude higher laser power might be required in the actual 
case. 


5.2.2.3 Microbunching Simulation 


Based on the above parameters, we have conducted the simulation of microbunching 
formation in the storage ring. The beam current and bunching factor one turn after 
the laser modulation is shown in Fig. 5.9. As can be seen from the comparison with 
Fig. 5.5, the simulation agrees with theory well. 


5.2.3 Microbunching Radiation Calculation 


Now we evaluate what radiation we can obtain from the formed microbunching. 
The numerical calculation of incoherent and coherent undulation radiation shown 
in this section are obtained using SPECTRA [13]. The beam energy and undulator 
parameters used are those in our SSMB proof-of-principle experiment, i.e., Eo = 
250 MeV, Az; = ào = 1064 nm, A, = 125 mm, K = 2.5, N, = 32. The results are 
also used to compare with the theoretical formulas presented in Chap. 4. Note that the 
numerical calculation of coherent radiation with a 3D charge distribution is usually 
time-consuming. This is also one of the motivations for us in Chap. 4 to develop the 
simplified analytical formulas with the main physics accounted for. 
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Fig. 5.9 Example current profile and bunching factor spectrum one turn after the laser modulation 
in SSMB PoP I, obtained from ELEGANT [12] tracking. Parameters used: op = 100 fs, €y = 31 
nm, €y = aie: 1 x 10° particles are simulated, meaning 160 fC for a one-to-one correspondence 


The left part of Fig. 5.10 shows the incoherent undulator radiation flux of 10° 
electrons (0.16 pC) versus the opening angle of a circular aperture placed in the 
forward direction of electron traveling. As can be seen, with the increase of the 
aperture opening angle, the red-shifted part of the radiation grows. For the total 
flux, there are sharp spikes near the odd harmonics and no clear spikes near the 
even harmonics. This is due to the fact that there is no on-axis radiation at the even 
harmonics. Also note that with the change of the aperture opening angle, there are 
jumps in the flux at a specific frequency nwo. This is due to the fact that the red- 
shifted radiation of higher harmonics m > n can contribute to the flux at œ = nwo 
when the aperture is large enough. 

Now we calculate the coherent radiation of the laser modulation-induced 
microbunched beam, using an RMS bunch length of 100 fs (o; = 30 um). An exam- 
ple beam current and bunching factor spectrum of the laser modulation-induced 
microbunched beam are shown in Fig. 5.9. We remind the readers that the bunch 
length in the real machine is typically longer than that used in the calculation here. 
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Fig. 5.10 Left: incoherent radiation photon flux of 10° electrons within a circular aperture placed 
in the forward direction, whose largest opening angles are Omax. Right: coherent radiation photon 
flux of 10° electrons within the circular aperture whose largest opening angle Omax = 8 mrad. The 
beam current and bunching factor spectrum used in the calculation are shown in Fig. 5.9. Other 
related parameters: Eg = 250 MeV, Ao = 1064 nm, A, = 125 mm, K = 2.5, N, = 32 


Therefore, the coherent radiation is even more narrowbanded than presented in the 
example calculation in this section. 

At first, we ignore the influence of electron beam’s transverse dimension, i.e., a 
thread beam is assumed. The coherent radiation spectrum of 10° electrons is shown 
in the right part of Fig. 5.10. As can be seen, there is narrowband coherent radiation 
only at the modulation laser harmonics, which fits with expectation due to the fact that 
there is only notable bunching factor at the laser harmonics. More closer look of the 
first two harmonics versus the aperture opening angle are shown in the upper part of 
Fig. 5.11. As can be seen, indeed the relative bandwidth of the coherent radiation at the 
fundamental frequency and second harmonic are 1% and 0.5%, respectively, agreeing 
with the values calculated from Eq. (4.21). In addition, the flux of the fundamental 
mode H = | atthe fundamental frequency agrees reasonably well with that according 
to Eq. (4.48), i.e., Fi(@ = Œœ, 0, = 0 um) = 3.4 x 107 (photons/pass/0.1% b.w.). 
In other words, the amplification factor of the flux at œ = Hap is indeed N2? b; ul? 
when o, = 0 pm. Also note that the jumps of the flux with the change of aperture 
opening angle as we commented just now. 

Now we investigate the impact of transverse electron beam sizes on the coher- 
ent radiation. As can be seen from the middle and bottom parts of Fig. 5.11, which 
correspond to a transverse electron beam size of 100 um and 400 um, respec- 
tively, and the comparison with the upper part, the transverse sizes of the electron 
beam suppress the coherent radiation. And the calculated fluxes at the fundamen- 
tal frequency agrees well with those predicted by Eq. (4.49), i.e., Fi (@ = wo, 0, = 
100 um) = 3.1 x 107 (photons/pass/0.1% b.w.) and Fi (@ = wo, o1 = 400 um) = 
1.6 x 10’ (photons/pass/0.1% b.w.). The suppression from transverse beam size is 
even more significant at the higher harmonics and the suppression factors agree with 
those predicted according to the transverse form factor Eq. (4.32). Also note that 
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Fig. 5.11 Coherent radiation photon flux of 10° electrons versus opening angle Omax of the cir- 
cular aperture for the first two harmonics, with o} = 0 wm (up), o1 = 100 um (middle) and 
oj; = 400 um (bottom), respectively. The beam current and bunching factor spectrum used in 
the calculation are shown in Fig. 5.9. Other related parameters: Ey = 250 MeV, Ao = 1064 nm, 


àu = 125 mm, K = 2.5, N, = 32 
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different from that of incoherent radiation, when o = 100 um or 400 um, there is 
no visible jump of the flux with the aperture opening angle Omax grown from | mrad 
to 8 mrad. This is because that the off-axis red-shifted coherent radiation of higher 
modes are suppressed now. 

We have also confirmed our derivation of the coherent radiation power by com- 
paring it with simulation. As shown in Fig. 5.12, the calculated peak powers of 
coherent radiation with different transverse electron beam sizes also agree well 
with the theoretical predictions from Eq. (4.46), i.e., Pi peak(@, = 0 um) = 363 W, 
Pi peak(@1 = 100 pm) = 332 W, Pi peak(o, = 400 um) = 168 W. 

From the calculation and analysis, we know that the coherent radiation from 
the formed microbunching is mainly at the fundamental frequency and second har- 
monic of the modulation laser, and in the forward direction. The coherent radiation 
is narrowbanded, and stronger than the incoherent radiation. These calculations and 
observation are the basis for our signal detection scheme. 


5.2.4 Signal Detection and Evaluation 


After investigating the microbunching formation beam dynamics and radiation char- 
acteristics of the formed microbunching in the above sections, now we consider how 
we can measure and evaluate the signals. 


Measurement and evaluation of bunch charge The bunch-by-bunch charge (cur- 
rent) in the experiment is measured by a single-bunch current monitor [14], which 
analyses the electron beam-induced RF signals from a set of four stripline elec- 
trodes (3 GHz bandwidth). To minimize the influence of neighbouring bunches on 
the signal, the pulse response of the electrodes is reshaped by a 500-MHz low-pass 
filter. The current calibration of the monitor is conducted using a parametric current 
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Fig. 5.13 Evaluation of bunch charge based on the stripe line measurement. Blue dots are the 
measurement results with the systematic offset subtracted and the red curve is a fit by the sum of 
two exponential functions, Q(t) = Qiexp(—t/t1) + Q2exp(—t/t2), performed at different time 
intervals, with the fit results connected by a smoothed line. (Figure from Ref. [1]) 


transformer [4] at higher current, and the linearity of the system at lower current 
is checked with the signal of the photodiode illuminated by synchrotron radiation. 
During the current decay in the experiment, one data point of the result given by the 
monitor is saved every second for each individual bunch. The averaged measurement 
of ten unfilled bunches preceding the homogeneous filled bunches (10 ns time gap 
in between) is used as the systematic offset. To smooth the measurement noise and 
at the same time account for the change of the beam lifetime, the time evolution of 
the offset-removed data points is then fitted by the sum of two exponential func- 
tions, Q(t) = Q,exp(—t/t1) + Q2exp(—t/t2), at different time intervals, with the 
fit results connected smoothly. One example evaluation of the bunch charge measure- 
ment result is presented in Fig. 5.13. Based on the evaluated data, we obtain a linear 
bunch-charge dependence of the broadband incoherent signal that is detected by the 
photodetector without the 3-nm-bandwidth band-pass filter, as shown in Fig. 5.14, 
confirming the reliability of the bunch-charge measurement and evaluation method. 
If the ring works in single-bunch mode, there is a more direct and accurate method 
of bunch charge measurement based on the measurement of synchrotron radiation 
strength using a photodiode. Both methods have been used in our experiments. 


Detection and evaluation of undulator radiation The long-pulse laser (FWHM 
= 10 ns) is used to simplify the experiment by avoiding a dedicated laser-electron 
synchronization system, given that the shot-to-shot laser timing jitter is fitter < 1 ns 
(RMS). However, the photodetector (Femto HSPR-X-I-1G4-SI; rise/fall time, 250 
ps) becomes saturated and even damaged by the powerful laser (Beamtech Optronics 
Dawa-200) if it is placed in the path of the laser. To address this issue, the undula- 
tor radiation is separated into the fundamental and second harmonics by appropri- 
ate dichroic mirrors (Thorlabs Harmonic Beamsplitters HBSY21/22), as shown in 


Fig. 5.3, and the signal detection at first focuses on the second harmonic with the 
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Fig. 5.14 Linear dependence of the broadband incoherent undulator radiation on the bunch charge. 
a, Results corresponding to individual laser shots; the shading (light red) represents 30 of the 
detection noise. b, The result after 200-consecutive-laser-shot averaging. The blue dots are the 
experimental data of a bunch not modulated by the laser and the red curves are linear fits. (Figure 
from Ref. [1]) 


wavelength centred at 532 nm. The photodetector output voltage, which is propor- 
tional to the radiation power, is then measured by a digital oscilloscope (Tektronix 
MS0O64:6-BW-4000; bandwidth, 4 GHz; sample rate, 25 billion samples per sec- 
ond). Later we will also present the result of the 1064 nm radiation by implementing 
Pockels cells to block the modulation laser and let pass the radiation in the following 
turns. 

An example radiation waveform of the second harmonic is shown in Fig. 5.16. To 
avoid the impact of the signal waveform offset caused by stray laser light, the data 
analysis takes the peak-to-peak value of the photodetector output voltage as a mea- 
sure of the radiation power. The coherent radiation power that corresponds to each 
individual laser shot, obtained during a time interval with a decaying beam current, is 
presented in Fig. 5.15a, where the modest contribution on the measured quantity from 
the small amount of incoherent radiation transmitted through the 3-nm-bandwidth 
band-pass filter has been eliminated. As can be seen, the coherent signal fluctuates 
considerably from shot to shot. This is attributable to the shot-to-shot fluctuation 
of the laser intensity profile (see Fig. 5.7) and the measurement noise. Despite the 
fluctuation, quadratic functions fit reasonably to the lower and upper bounds of the 
data points, which correspond to the cases of minimum and maximum bunching 
factors induced by the fluctuating laser, respectively. When performing the fits, we 
took into account that the measured quantity is the real radiation signal convoluted 
with the detection noise. The impact of this noise, obtained by analysing the mea- 
surement result of the unfilled bunches, on the bounds of the measured data points is 
visualized as shading in Fig. 5.15a for the coherent signal and in Fig. 5.14a for the 
incoherent signal. To smooth this shot-to-shot fluctuation, a 200-consecutive-laser- 
shot averaging is conducted and the results are presented in Figs. 5.14b and 5.15b 
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Fig. 5.15 Quadratic dependence of the narrowband coherent undulator radiation generated from 
microbunching on the bunch charge. a, Results corresponding to individual laser shots; the shading 
(light red and grey) represents 30 of the detection noise. b, The result after 200-consecutive-laser- 
shot averaging; the plot is the same as Fig. 5.17 and is presented again here for comparison with a 
and with the incoherent signal in Fig. 5.14. The blue dots represent the experimental data and the 
red curves are quadratic fits. (Figure from Ref. [1]) 


for the narrowband coherent and broadband incoherent signals, where a quadratic 
and a linear fit have been performed, respectively. 


5.2.5 Experimental Results 


5.2.5.1 Second Harmonic Radiation 


Figure 5.16 shows the typical measurement results of the second-harmonic undulator 
radiation emitted from a homogeneous stored bunch train, with a charge of about 
1 pC per bunch and a time spacing of 2 ns, supplied by the 500 MHz RF cavity at 
the MLS. The spikes in the waveforms are the signals of different bunches. The left 
and right panels show the results corresponding to 2 and 40 consecutive laser shots, 
respectively. To smooth the measurement noise and signal fluctuation, the waveforms 
in the right panels have been averaged. Figure 5.16a, b shows the signals one turn 
before the laser shot, which correspond to the incoherent radiation and reflect the 
homogeneous bunch filling pattern. Figure 5.16c, d shows the radiation one turn after 
the laser shot, from the same bunches as those in Fig. 5.16a, b. The five larger spikes 
at the centre correspond to the bunches modulated by the laser. The enhanced signals 
of these five spikes indicate the formation of microbunching and the generation of 
coherent radiation from the laser-modulated bunches. 

The laser used in the experiment has multiple longitudinal modes, and its temporal 
profile has several peaks and fluctuates considerably from shot to shot (see Fig. 5.7). 
Therefore, the laser-induced electron energy modulation amplitudes are different 
from shot to shot and from bunch to bunch. When the modulation amplitude matches 
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Fig. 5.16 Waveforms of the undulator radiation produced from a homogeneous stored bunch train. 
a, b, Radiation one turn before the laser shot. The photodetector output voltage is proportional to 
the radiation power. c, d, Radiation one turn after the laser shot, from the same bunches as those of a 
and b, where the central five bunches are modulated by the laser pulse. The offset and general slight 
decreasing trend of the waveforms are due to the photodetector being saturated by stray light from 
the modulation laser one revolution before and not having completely recovered. e, f, Radiation 
one turn after the laser shot, obtained with a narrow band-pass filter (centre wavelength, 532 nm; 
bandwidth, 3 nm FWHM) placed in front of the photodetector, with bunch filling and charge similar 
to those in a to d. (Figure from Ref. [1]) 


the phase slippage factor, the energy-modulated electrons are properly focused at 
synchronous phases, which gives optimal microbunching. For some of the shots, 
the laser intensity is higher or lower than the optimal value, and the electrons are 
then over-focused or under-focused, giving weaker microbunching and less coherent 
radiation. This explains the shot-to-shot fluctuation of the coherent amplified signals 
shown in Fig. 5.16c, e. 

As analyzed before, the microbunching coherent radiation is much narrowbanded 
compared to the incoherent radiation. To confirm that the amplified radiation is due to 
microbunching, we tested this narrowband feature of the coherent radiation. A band- 
pass filter (Thorlabs FL532-3; centre wavelength, 532 nm; bandwidth, 3 nm FWHM) 
was inserted in front of the detector. The radiation one turn after the laser shot is 
shown in Fig. 5.16e, f, which was obtained with a bunch filling and charge similar to 
that of Fig. 5.16c, d. From the comparison between Fig. 5.16e and c (Fig. 5.16f and 
d), we can see that the broadband incoherent signals are nearly completely blocked 
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by the filter, whereas the amplified part is not affected much, confirming that the 
amplification is the narrowband coherent radiation generated by the microbunches. 

Finally, we investigated the dependence of the coherent radiation on the bunch 
charge. To mitigate collective effects such as intrabeam scattering and head-tail insta- 
bility, which could change the electron beam parameters, this investigation was con- 
ducted at low beam current, and the coherent signal was optimized by fine-tuning the 
machine to ensure a sufficient signal-to-noise ratio. Because the longitudinal radia- 
tion damping time in the experiment was 180 ms, we operated the laser at 1.25 Hz 
repetition rate to ensure that the electron bunches had time to recover their equilib- 
rium parameters before each laser shot. The 3-nm-bandwidth band-pass filter was 
inserted to block the incoherent radiation, and the coherent signal corresponding to 
each individual laser shot was saved, with the beam current decaying naturally until 
the signal was at the detection noise level. The measurement results of the bunch 
closest to the laser temporal centre (t = 0 ns in Fig. 5.16) are used for quantitative 
analysis as introduced above. To lessen the impact of the laser temporal profile fluctu- 
ation and measurement noise, a 200-consecutive-laser-shot averaging is performed to 
obtain the data point for each bunch charge. The coherent undulator radiation power 
versus the single-bunch charge is shown in Fig. 5.17, where a quadratic function fits 
well to the experiment data. The quadratic bunch charge dependence, together with 
the narrowband feature of the coherent radiation, demonstrates unequivocally the 
formation of microbunching. 
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Fig. 5.17 Quadratic dependence of the coherent undulator radiation generated from microbunching 
on the bunch charge. The blue dots represent the experimental data and the red curve is a quadratic 
fit. Each data point represents the averaged result of 200 consecutive laser shots. The error bars 
denote the standard deviation of the averaged results when the averaging time window shifts for 
+100 consecutive laser shots from the corresponding data point. (Figure from Ref. [1]) 
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5.2.5.2 Fundamental Frequency Radiation 


The above deciding experimental results were obtained in the year of 2020, and 
have been published in Ref. [1]. After that, there are two main upgrades on the 
experimental setup. First, the multi-longitudinal-mode laser has been replaced by a 
single-longitudinal-mode one (Amplitude Surelite I-10). Second, Pockel Cells have 
been installed along the signal detection optical path to block the modulation laser 
and let pass the radiation in later revolutions, thus allowing the detection of the 
fundamental frequency radiation [15]. As shown in our analysis, we expect that the 
coherent radiation at the fundamental frequency is much stronger than that at the 
second harmonic and the microbunching can last multiple turns. These expectations 
have been confirmed in our following experimental investigations at the MLS [2]. 

Figure 5.18 shows the typical experimental results of the multi-turn coherent radi- 
ation at the fundamental frequency. A narrowband-pass filter has been inserted to 
select the narrowband coherent radiation. Figure 5.19 is the more quantitative data 
analysis of the signal of the first three turns after each laser shots. The bunch charge 
has now also been obtained in a more accurate way by using the incoherent syn- 
chrotron radiation signal measured by a photodiode. Several important observations 
are in order concerning the experiment results: 


e First, signals of all three turns have shown nice quadratic bunch charge fits, con- 
firming again the formation of microbunching and coherent radiation generation 
from it. 

e Second, we mentioned in the above section and also in Ref. [1] that the huge shot- 
to-shot coherent radiation signal fluctuation obtained before the upgrades is mainly 
due to the laser profile fluctuation arising from its multi-longitudinal-mode nature. 
This argument has also been confirmed by the results shown in Fig. 5.19. From 
the comparison of Figs. 5.15a and 5.19, we can see that the shot-to-shot coherent 
signal with the present single-longitudinal-mode laser is much more stable. 
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Fig. 5.18 Raw data of the multi-turn microbunching preservation experiment result. The signal is 
for the fundamental-mode undulator radiation, i.e., 1064 nm. (Refer to J. Feikes’ talk in IPAC2021 
and Ref. [2] for more details) 
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Fig. 5.19 The bunch-charge scaling of the coherent undulator radiation signals of the first three 
turns after the laser shots. Up: each data point corresponds to the radiation signal strength after 
each single laser shot. Bottom: average of the radiation signal from 40 consecutive laser shots. The 
saturation level of the detector is about 2 V 


e Third, the signal deviates from the quadratic scaling at high current and starts 
to saturates about 1.3 pC, which we believe is due to the influence of collective 
effects. 


More in-depth investigations on the multi-turn microbunching is still ongoing and 
will be reported in the future [2]. 


5.2.6 Summary 


In conclusion, we have demonstrated the mechanism of SSMB for the first time 
in a real machine. This demonstration represents the first milestone towards the 
implementation of an SSMB-based high-repetition, high-power photon source. As 
great as the experimental results are, to avoid confusion, here we make clear that 
here we do not report an actual demonstration of SSMB, but rather a demonstration 
of the mechanism by which SSMB will eventually be attained. First, the formation 
of microbunching after one complete revolution of a laser-modulated bunch in a 
quasi-isochronous ring and the maintenance of microbunching for multiple turns 
demonstrate the viability of a turn-by-turn electron optical phase correlation with a 
precision of sub-laser wavelength. Second, this microbunching is produced on the 
stored electron bunch, the equilibrium parameters and distribution of which before the 
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laser modulation are defined by the same storage ring as a whole. The combination 
of these two crucial factors establishes a closed loop to support the realization of 
SSMB, provided that a phase-locked laser interacts with the electrons turn by turn. 


5.3 PoP II: Quasi-steady-State Microbunching 


On the basis of the PoP I, the next step is to replace the single-shot laser by a high- 
repetition phase-locked one to interact with the electrons turn-by-turn. By doing so, 
we want to form stable microbuckets to constrain the microbunching in it to reach 
a quasi steady state. This is the SSMB PoP II as introduced in the beginning of this 
chapter. 


5.3.1 Phase-Mixing in Buckets 


To reach a quasi steady state, the particles need to do synchrotron oscillations to reach 
phase mixing in the microbuckets, as a result of longitudinal amplitude dependent 
tune spread of the electron beam. Here we present a remarkable feature of phase 
mixing or filamentation in RF or optical buckets. As we will see soon, given an initial 
DC mono-energetic beam, there will be an equilibrium phase space distribution after 
phase mixing in the bucket. We find that in this final steady state, the beam current 
distribution has little dependence on the bucket height. This feature is favorable for 
the SSMB PoP II, as the requirement on the modulation laser power can then be 
much relaxed compared to PoP I. This effect is also of relevance to the injection 
process of the final real SSMB storage ring. 

As the phase mixing is a rather fast process compared to radiation damping, we 
consider only the symplectic dynamics in this section for simplicity. The symplectic 
longitudinal dynamics of a particle in a storage ring with a single RF system, in 
SSMB a laser modulator, can be modeled by the well-known “standard map” [16] 


Tn41 = In +K sin On 


(5.15) 
On+1 = On + In41) 
in which 
VRF 
0 = kgrz, I = Rsekęrô, K = E, Boker, (5.16) 
0 
with R55 = —nCo. Note that K in this section is not the undulator parameter. Equa- 


tion (5.15) can be described with the pendulum Hamiltonian driven by a periodic 
perturbation 
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Íy = 
H(I,0,1) = 59° + K cos XC cos(2znt). (5.17) 


n=— 00 


The dynamics is given by a sequence of free propagations interleaved with periodic 
kicks. For K 4 0, the dynamics is non-integrable and chaotic. But for a K much 
smaller than 1, which is the case for usual storage rings working in the longitudi- 
nal weak focusing regime, the motion is close to integrable and the differences in 
Eq. (5.15) can be approximately replaced by differentiation, and the Hamiltonian 
Eq. (5.17) can be replaced by a pendulum Hamiltonain 


1 
H = 51 +K cosð. (5.18) 


The separatrix of the pendulum bucket is H = K with a bucket half-height of 2/K. 
Or in unit of ô, the bucket half-height is 
61 


_ 24 VerRsoker 1 


2 | Rsokre| Beskrr’ 


(5.19) 


where f,s is the longitudinal beta function at the RF as analyzed in Sect. 2.1.2. The 
synchrotron tune is 
VK 


Vs & —sgn(K)- (5.20) 


When K is large, the strongly chaotic dynamics can also be used for interesting 
applications, for example applying the bucket purification to generate short bunches 
as proposed in Ref. [17]. 

Figure 5.20 shows a simulation result of the evolution in the longitudinal phase 
space of a mono-energetic DC beam after injection into RF or optical buckets 
described by Eq. (5.15). We have chosen to observe the beam in the middle of the 
RF kick so the beam distribution in the longitudinal phase space is upright. As can 
be seen, there is a steady-state beam distribution due to phase mixing in the bucket. 
Note that the bucket center is at@ = 2 when K > 0.If K < 0, then the bucket center 
will be at 0 = 0. 

As the longitudinal form factor, thus the coherent radiation power, depends more 
directly on the beam current (namely the longitudinal coordinate z of the electrons) 
rather than the energy spread, now we try to get an analytical formula for the steady- 
state beam current. For convenience, we shift the bucket center to the origin, which 
means 0 — x — 0, or a sign change of K. What we want to know is the steady- 
state distribution of 0, i.e., f(0,t —> oo). In action-angle (¢, J) phase space, the 
distribution function evolves according to 


S@, J,t) = f(@—- w(J)t, J, 0). (5.21) 
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Fig. 5.20 Phase mixing (filamentation or decoherence) of a mono-energetic particle beam trapped 
in RF or optical bucket with K = 0.01. Up: particle distribution in longitudinal phase space, with 
red curves being the separatrices. Bottom: the corresponding beam current distribution 


Fig. 5.21 A plot to help 
better understand Eq. (5.23) 


When there is a tune dependence w(J) on J, then in the limit of t + oo, the steady- 
state distribution depends only on the initial distribution of action J as a result of 
phase mixing, 


2n 


f(¢, J,t > 0) = L fle, J, t = 0)d¢. (5.22) 
27 0 


The final angle for each action J will uniformly distributed in [0, 277). As shown in 
Fig 5.21, after reaching the steady state, the percentage of the particles with 0 < x 
for0 <x < zis 


P@<x)=~+4 f (1 _ ees >) 1 iB. (5.23) 
x Jy m m 


in which Z (x, 8) represents the 7-coordinate of a point on the (ġ, J) phase space 
trajectory traversing (£, 0) with a 6-coordinate of x. 
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After getting P(@ < x), the current distribution can then be calculated according 


to 
aP 
f@)= Pr (5.24) 
x 


x=0 


However, ġ (x, I (x, B)) has a complex form, and it is hard to get a simple analytical 
expression for f (0). Here we simplify the discussion by approximating all the phase 
space trajectories in the bucket by ellipses to arrive at an analytical formula for f (0). 
For an ellipse phase space trajectory, we have 


(x, I(x, B)) = arccos a (5.25) 


Note that the result in Eq. (5.25) has no dependence on K. For a real RF or optical 
bucket, there is a dependence of ġ (x, I (x, 6)) on K, but the dependence is weak, 
especially for trajectories close to the bucket center. So we expect our approximated 
Eq. (5.25) is valid to a large extent. Substituting Eq. (5.25) into Eqs. (5.23) and (5.24), 
we have 


n+ /n?—@ 
0 


oP 
f@) = ay 


X 


E i 1 TEE 
=f 9\2 B aoe 
1-(§) 
(5.26) 


Note that our simplified theoretically current distribution f (0) is independent of K, 
which means the steady-state current distribution is independent of the bucket height. 

Figure 5.22 shows the simulation result of the steady-state current distribution 
under different K, i.e., different bucket heights, and simultaneously our simplified 
theoretical distribution Eq. (5.26). As can be seen that indeed the steady-state current 
distribution has little dependence on the bucket height, and our simplified analysis 
is quite accurate. Note that the origin is not shifted in the plot. 

The analysis reveals a remarkable feature of phase mixing in RF or optical bucket, 
i.e., the final steady-state current distribution after a mono-energetic beam getting 
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Fig. 5.23 The steady-state 8 <<a 
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trapped by RF or optical buckets has little dependence on the bucket height. This 
is helpful for our Quasi-SSMB experiment since it means the requirement on the 
modulation laser power is not that demanding. A bucket height several times of the 
natural energy spread is sufficient. 

The above result is based on a constant RF voltage in the phase mixing process, 
it can be anticipated that more particles will be bunched closer to the bucket center 
phase when we increase K after the beam reach its steady-state distribution after 
phase mixing. Similar steps to the above section can be invoked for calculating the 
new steady-state current distribution. A transformation of the action when K changes 
is all that needed. Figure 5.23 shows the simulation result of the steady-state current 
distribution by increasing K in two consecutive steps from 0.001 to 0.03. As can be 
seen, the current are more concentrated to the center after the increase of K. 

A discrete change of K can boost bunching as shown in Fig. 5.23. However, 
it is not without sacrifice, as the filamentation process will result in longitudinal 
emittance growth. This emittance increase is unwanted in some cases. As well-studied 
in RF gymnastics [18], an adiabatic change of RF voltage or lattice parameters can 
manipulate the bunch length while preserving the longitudinal emittance. Similar 
ideas can also be applied to boost microbunching with little emittance growth. A 
simulation of trapping of microbunch with K linearly ramped from 1 x 10~° to 1 x 
107? is shown in Fig. 5.24. Note the drastic difference between Figs. 5.24 and 5.20. 
The spirit of adiabatic buncher [19, 20] is the same with adiabatic trapping introduced 
here, for enhancing microbunching while preserving longitudinal emittance which 
is useful for FEL and inverse FEL. The adiabatic trapping mechanism can also be 
applied in the beam injection of the SSMB or other storage rings whose momentum 
aperture is of concern. It is interesting to note the connection of adiabatic trapping 
with the microbunching process in a high-gain FEL [21-23]. 
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Fig. 5.24 Trapping of particles with a linear increase of K from 1 x 1076 to 1 x 107? in 10+ 
turns. Up: particle distribution in longitudinal phase space, with red curves being the separatrices. 
Bottom: the corresponding beam current distribution 


5.3.2 Experimental Parameters Choice 


As the quantum diffusion of z is large for the MLS lattice (26 nm RMS per turn 
at 250 MeV corresponding to optics in Fig. 5.2), it is not feasible to realize true 
SSMB inside a 1064 nm wavelength microbucket at the MLS. Therefore, the goal 
of SSMB PoP II is to accomplish microbunching for 100 to 1000 consecutive turns 
to reach a quasi steady state. Based on the beam physics and noises analysis, the 
tentative experimental parameters choice is as shown in Table 5.2. We have conducted 
numerical simulations based on the parameters set, from which we observe that the 
typical evolution of electrons in PoP II experiment can be divided into several stages. 


e I: with the modulation laser turned on, the bunching factor reaches the maximum 
after about one quarter of the synchrotron oscillation period; 

e II: after several synchrotron oscillation periods, the whole microbuckets are filled 
with particles like that shown in the right part of Fig. 5.20 as a result of phase 
mixing. Bunching factor after reaching this quasi-steady state will be stable if 
there is no quantum excitation or other diffusion effects; 

e III: due to quantum excitation and various diffusion effects, energy spread starts 
to increase and particles continue to leak out the microbuckets and begin to hit on 
the vacuum pipe and become lost. Bunching factor in this stage decreases; 

e IV: after a while, all the particles are lost in the end. 


To make the experiment more realistic, each time we fire the laser, we only want 
to accomplish the above stages I and II, but avoid III and IV, i.e., to avoid particle 
loss, otherwise it will be too time-consuming to do the experiment. This is based on 


148 5 SSMB Proof-of-Principle Experiments 


Table 5.2 Tentative parameters of the Quasi-SSMB experiment to be conducted at the MLS 


Parameter Value Description 

Eo 250 MeV Beam energy 

Co 48 m Ring circumference 

ÍRE 500 MHz RF frequency 

n In| < 2 x 1075 Phase slippage factor 

T 180 ms@250 MeV Longitudinal radiation damping time 

os 1.8 x 1074@250 MeV | Natural energy spread 

Ex 31 nm@250 MeV Horizontal emittance 

Àu 125 mm Undulator period length 

Nu 32 Number of undulator periods 

Ly 4m Undulator length 

K 2.5 Undulator parameter 

AL 1064 nm Modulation laser wavelength 

ZR ~ tu Rayleigh length 

Poeak >10kW Modulation laser peak power 

é1 >1.505 Microbucket half-height 

AR 1064 nm Radiation wavelength 

by >0.01 Bunching factor (1064 nm) in 
quasi-steady state 


the fact that preparing the beam and storage ring state is time-consuming, while the 
particle can be lost in milli seconds with the laser keep firing. This is why we aim for 
preserving mircorbunching for 10° turns, instead of 10° turns or a longer time. The 
experiment is under preparation and more progress will be reported in the future. 
One thing worth mentioning is that the second-harmonic bunching in the quasi-steady 
state is negligible at the MLS, therefore, fundamental frequency radiation detection 
is needed in SSMB PoP II. 
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Chapter 6 A) 
Summary get 


In this final chapter, we give a brief summary of the dissertation, and present some 
useful results for practitioners. 


6.1 Summary of the Dissertation 


The contribution of this dissertation consists of three parts: in Chaps. 2 and 3, we 
answer the question of how to realize SSMB; in Chap. 4, we investigate what radi- 
ation characteristics can be obtained from the formed SSMB; and in Chap. 5, we 
experimentally demonstrate the working mechanism of SSMB in a real machine for 
the first time. 

In Chap. 2, to account for the impact of local phase slippage factors on beam 
dynamics in a quasi-isochronous electron storage ring, we have developed and 
applied the Courant-Snyder formalism in longitudinal dimension to derive new for- 
mulae of bunch length, energy spread and longitudinal emittance beyond the classical 
scaling laws. The method of optimizing the global and local phase slippages simul- 
taneously to minimize the longitudinal £ function at the bending magnets has been 
proposed based on the analysis, to generate an ultrashort bunch length and ultrasmall 
longitudinal emittance, as required by SSMB. Further, we have derived the scaling 
law of the theoretical minimum bunch length and longitudinal emittance with respect 
to the bending radius and angle of the bending magnet. The use of transverse gra- 
dient bends for minimizing the longitudinal emittance has also been investigated. 
The application of multiple RF cavities, or laser modulators in an SSMB storage 
ring, for longitudinal strong focusing has been discussed using the same formalism 
with important observations made. Considering the momentum compaction of a laser 
modulator, its thick-lens linear and nonlinear maps have been derived and simulated 
for a more accurate modeling of beam dynamics in it. We have also studied the 
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application of the higher-order terms of phase slippage for high-harmonic bunching 
and longitudinal dynamic aperture optimization. Based on the investigations in this 
chapter, we have presented in Table 6.1 an example parameters set of a longitudinal 
weak focusing SSMB storage ring for high-power infrared radiation generation. 

In Chap. 3, we have presented a concise analysis of the bending magnet-induced 
passive bunch lengthening from transverse emittance of the particle beam. After that, 
we have generalized the analysis and proved three theorems on the active applications 
of transverse-longitudinal coupling (TLC) for efficient harmonic generation or bunch 
length compression. These theorems dictate the relation between the modulation lick 
strength and the lattice optical functions at the modulator and radiator, respectively. 
Further, we have analyzed the contribution of modulators to the vertical emittance 
from quantum excitation, to obtain a self-consistent evaluation of the required mod- 
ulation laser power in applying these TLC schemes in a storage ring. These theorems 
and related analysis provide the theoretical basis for the application of TLC in SSMB 
to lower the requirement on the modulation laser power, by taking advantage of the 
fact that the vertical emittance in a planar ring is rather small. The relation between 
our TLC analysis and the transverse-longitudinal emittance exchange is also briefly 
discussed. In addition to the investigation on linear TLC dynamics, we have also 
reported the first experimental validation of particle energy widening and distortion 
by anonlinear TLC effect in a quasi-isochronous ring, which originates from an aver- 
age path-length dependence on the betatron oscillation amplitudes. The result could 
be important for quasi-isochronous rings, SSMB, nonscaling fixed-field alternate 
gradient accelerators, etc., where very small phase slippage factor or large chro- 
maticity is required. Based on the investigations in this chapter, we have presented 
in Table 6.2 an example parameters set of a transverse-longitudinal coupling SSMB 
storage ring for high-power EUV and soft X-ray radiation generation. 

In Chap. 4, we have presented theoretical and numerical studies of the average 
and statistical property of the coherent radiation from SSMB. Our results show that 
kW-level average power of 13.5 nm-wavelength EUV radiation can be obtained from 
an SSMB ring, provided that an average current of | A and bunch length of 3 nm 
microbunch train can be formed at the radiator. Such a high-power EUV source 
is a promising candidate to fulfill the urgent need of semiconductor industry for 
EUV lithography. Together with its narrowband feature, the EUV photon flux can 
reach 10° ~ 10!° phs/s within a 0.1 meV energy bandwidth, which is appealing 
for fundamental condensed matter physics research. In the theoretical investigation, 
we have generalized the definition and derivation of the transverse form factor of an 
electron beam which can quantify the impact of its transverse size on the coherent 
radiation. In particular, we have shown that the narrowband feature of SSMB radiation 
is strongly correlated with the finite transverse electron beam size. Considering the 
pointlike nature of electrons and quantum nature of radiation, the coherent radiation 
fluctuates from microbunch to microbunch, or for a single microbunch from turn to 
turn. Some important results concerning the statistical property of SSMB radiation 
have been presented, with a brief discussion on its potential applications for example 
the beam diagnostics. The presented work is of value for the development of SSMB 
and better serve the potential synchrotron radiation users. In addition, it also sheds 
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light on understanding the radiation characteristics of free-electron lasers, coherent 
harmonic generation, etc. 

In Chap. 5, we have reported the first demonstration of the mechanism of SSMB 
at the Metrology Light Source in Berlin. We have shown that electron bunches stored 
in a quasi-isochronous ring can yield sub-micrometre microbunching and narrow- 
band coherent radiation, one complete revolution after energy modulation induced 
by a 1064 nm wavelength laser, and this microbunching can preserve for multiple 
turns. These results verify that the optical phases, i.e, the longitudinal coordinates, 
of electrons can be correlated turn by turn in a storage ring at a precision of sub-laser 
wavelengths. On the basis of this phase correlation, we expect that SSMB will be 
realized by applying a phase-locked laser that interacts with the electrons turn by 
turn. This demonstration represents the first milestone towards the implementation 
of an SSMB-based high-power, high-repetition photon source. 


6.2 Useful Formulas and Example Parameters for SSMB 
Storage Rings 


To make our investigations more useful for practitioners, especially concerning the 
parameters choice for an SSMB storage ring, here we present some important for- 
mulas. Generally we group our formulas into two categories, i.e., a longitudinal 
weak focusing storage ring for a desired radiation wavelength Xr = 100 nm, and 
a transverse-longitudinal coupling storage ring for a desired radiation wavelength 
1 nm SAr < 100 nm. In each category, we have presented an example parameters 
set for the corresponding SSMB storage ring. 


6.2.1 Longitudinal Weak Focusing SSMB 


The relation of bending radius p and magnetic flux density B of the bending magnet 


ss 1 B[T] 
— = 0.2998 —__—_, (6.1) 
p Eo[GeV] 


with E the electron energy. 
Assuming that the storage ring consists of isomagnets, then the radiation loss of 
an electron per turn is 
Et 
Regs (6.2) 


Pring 


with C, = 8.85 x 1075 eVi? Pring the bending radius of bending magnets in the ring. 
The horizontal, vertical and longitudinal radiation damping constants for a planar 
uncoupled ring are 
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g (-—2n) Dy ds 
where D = a with n = =F p the field gradient index and D, is the hor- 


izontal dispersion. Nominally for a planar uncoupled ring using bending magnets 
with no transverse gradient, we have D < 1. 
The horizontal, vertical and longitudinal radiation damping times are 
Co / C 


+ 
QH,V.L 


(6.4) 


THV.L >= 


with Co the ring circumference and c the speed of light in vacuum. 
The natural energy spread of electron beam in a longitudinal weak focusing ring 


is 
Cy 
oss = E (6.5) 
Js p 


with C4 = = > = 3.8319 x 10-3 m, a, = + = 386 fm is the reduced Compton 
wavelength of electron, J; = 2 + D is the foneitidinal damping partition number, 
y is the Lorentz factor. 

The natural bunch length at the laser modulator is 


0-5 = Oss fzS, (6.6) 


where £5 is the longitudinal beta function at the laser modulator to be given soon. 
The effective modulation voltage of a laser modulator using a planar undulator 


is [1] 
JJIK [4P,ZoZ La 

ğal l AAN E an= ( ): (6.7) 
y ÀL 2Zr 


in which [J J] = Jo(x) — Jj (x) and x = en J, is the n-th order Bessel function 
of the first kind, K = <2- = 0.934. Bo[T] -A,[cm] is the undulator parameter, 


determined by the peak Ma Atie flux density Bo and period i, of the undulator, Pz 
is the modulation laser power, Zo = 376.73 Q is the impedance of free space, Zp is 
the Rayleigh length of the laser, L,, is the undulator length. 

The linear energy chirp strength around zero-crossing phase is related to the laser 


and modulator undulator parameters according to 
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V JIJIK |4PLZoZ L, 
pa RAA EER l kı, (6.8) 
Eo y2mc? AL 2ZR 


where ky; = 27 /À; is the wavenumber of the modulation laser. 
Linear stability of the longitudinal motion requires 


0 < hnCo < 4, (6.9) 


where 77 is the phase slippage factor of the ring. 

Considering the fact that the modulation waveform is sinusoidal and the longi- 
tudinal dynamics is more accurately modeled by a “standard kick map”, to avoid 
strong chaotic dynamics, an empirical criterion is 


0<hnCy S 0.1. (6.10) 
In a longitudinal weak focusing ring (v; < 1), the synchrotron tune is 


Ji 
pe Ne (6.11) 


In a longitudinal weak focusing ring, the longitudinal beta function at the laser 


modulator is 
nCo 
bs ©] (6.12) 


A 2 
ôi = ; 
4 PzskL 


The micro-bucket half-height is 


(6.13) 


If there is a single RF or laser modulator in the ring, and J; = 2, then the theoretical 
minimum bunch length and longitudinal emittance in a longitudinal weak focusing 
ring with respect to the bending radius p and angle 6 of each bending magnet are 


o; minum] © 4.93p? [m] Eo[GeV]6° [rad], 


> ; (6.14) 
€; minlnm] © 8.44E2[GeV]6? [rad]. 


Scaling law of the horizontal emittance in an SSMB storage ring is 


e [nm] ~ —366.5E2[GeV]9° [rad] a Px + £ cot Px (6.15) 
ay aA ae 9 2 10 ay 
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with ®, the horizontal betatron phase advance per cell which usually lies in (x, 27). 
The above scaling is derived by considering only the contribution of main cells, and 
ignoring that from the matching section. 

Coherent undulator radiation power at the odd- H -th harmonic from a transversely- 
round electron beam is 


Pit peakKW] = 1.183N, H X [J J Jy F F1 (S)lbz n| IRIA], (6.16) 
2 
where N, is the number of undulator periods, JI = [ys (Hx) - Jus (Hx)] ; 


with x = ae and the transverse form factor is 
FF, (5) a tan! : +51 Qs) (6.17) 
= — | tan — n| ——— ]], ; 
g x 25 (2S)? +1 
with S$ = ae and o; the RMS transverse electron beam size, b; y is the bunching 


factor at the H-th harmonic, and Tp is the peak current. 
The relative fluctuation of coherent radiation power considering the pointlike 
nature of electrons is 


(bk)? Ne |b(k)|4 


Var [|b(K)|?] 2 [b(k)? + Re [B2k)b (—k) | ; +o( 3 


(6.18) 
where b(k) is the bunching factor at the wavevector k, and Ne is the number of 
electrons. 

Based on the above formulas, here we present an example parameters set in 
Table 6.1 of a longitudinal weak focusing SSMB storage ring, aimed for high-power 
infrared radiation generation. As can be seen, such a compact SSMB storage ring 
can be used for power amplification of the injected seed laser. The requirement on 
the stored laser power is easy to realize in practice. All the other parameters are also 
within practical range. A sharp reader may notice that the microbucket half-height 
is only twice the natural energy spread of the electron beam. Therefore, in addition 
to these shallow microbuckets, we need a larger bucket, for example a barrier bucket 
formed by an induction linac, to constrain the particles in the ring to ensure a large 
enough beam lifetime. 


6.2.2 Transverse-Longitudinal Coupling SSMB 


For a transverse-longitudinal coupling (TLC) based SSMB, or a generalized longi- 
tudinal strong focusing SSMB [2], using TEMO00 mode laser modulator for energy 
modulation, we have the following important formulas. 
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Table 6.1 Example parameters set of a longitudinal weak focusing SSMB storage ring for infrared 


radiation generation 


Parameter Value Description 

Eo 250 MeV Beam energy 

Co 50m Circumeference 

n 4x 1076 Phase slippage factor 

Pring 0.6 m Bending radius of dipoles in the ring 

Bring 1.39 T Bending field in the ring 

0 7 Bending angle of each dipole 

OSS 2.76 x 1074 Natural energy spread 

Oz lim 86 nm Theoretical lower bunch length limit 

ÀL 1064 nm Modulation laser wavelength 

h 500 m7! Energy chirp strength 

ozs 175 nm Natural bunch length 

ô 1 5.36 x 1074 Microbuket half-height 

uMod 5cm Modulator undulator period 

BoMod 0.92 T Modulator peak magnetic field 

LuMod lm Modulator undulator length 

PL(ZR = Lu) 24 kW Modulation laser power 

g 5 x 103 Optical enhancement cavity gain 

Pin 4.8 W Injection laser power 

ÀR =S ÀL 1064 nm Radiation wavelength 

by 0.59 Bunching factor 

oL 100 um Transverse electron beam size at the 
radiator 

ÀuRad 5 cm Radiator undulator period 

Borad 0.92 T Radiator peak magnetic field 

LuRad 2m Radiator length 

Pr 1kW @ Ip = 0.55 A Radiation peak/average power 


Relation between energy chirp strength and optical functions at the modulator 


and radiator 


h’ (Mod)H, (Mod)H, (Rad) > 1, (6.19) 


where H, is a chromatic function quantifying the contribution of vertical emittance 


to bunch length. 
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Put the above relation in another way, 


Ey 


h> : (6.20) 
ozy (Mod)ozy (Rad) 


Bunching factor at the n-th laser harmonic in TLC SSMB at the radiator 


[0,6] 
b= ( Y Jin (n) exp [- ((n — m)kLoz (Mod))? z) exp [- (nkro; (Rad))? /2| 


m>=—Oo 
(6.21) 
where o, (Mod) = V €z: (Mod) + €y Hy (Mod) and o; (Rad) = „/€y H, (Rad) are the 
linear bunch length at the modulator and radiator, respectively. 
Contribution of two modulators to €, from quantum excitation 


55 apa2y> H, (Mod) 4 


Ae, (Mod) = 2 x 
i 96/3 æv Pma 3T 


Lis (6.22) 


where œp = a is the fine-structure constant. 
Assuming €, = Ae,(Mod), which means the vertical emittance is solely from 
the two modulators, then the required modulation laser power and modulator length 


scaling are 
8 7 
p, [kW] ~ 5.6722 lnm) Eo [GeV] BomoalT] 
ge (Rad)[nm] Bring [T] 


Ly [m] x57 Bring[T]Ey La , 
H, (Mod)[um] BeyjoalT] 


$ 


(6.23) 


where Bomoa is the peak magnetic flux density of the modulator undulator, B,ing is 
the magnetic flux density of bending magnets in the ring. The above scaling laws are 
accurate when K, > /2. For the more general case, refer to Eq. (3.56). 

Based on the presented formulas, here we present an example parameters set in 
Table 6.2 of a TLC SSMB storage ring, aimed for high-power EUV and soft X-ray 
radiation. It can be seen that as long as we can realize a coasting beam of 1.5 A 
average current, and an optical cavity stored power of 3100 kW, we can realize 1 
kW average power 13.5 nm EUV and 6.75 nm soft X-ray radiation. All the other 
parameters applied should be realizable, including the small €, considering IBS. 
Even if we can only realize an average beam current of 1 A or less, we can take 
advantage of the fact that Poon « I 3 to realize an average radiation power of kW 
level, by decreasing the filling factor of electron beam in the ring but increasing the 
peak current as long as the value is below the collective instability threshold.! Since 
there is no requirement on the longitudinal emittance for a coasting beam, thus no 
requirement on the fine control of longitudinal 6 function, the circumference of this 
ring has great flexibility, which means the ring can be very compact, for example a 


! Private communication with Alex Chao. 
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Table 6.2 Example parameters set of a transverse-longitudinal coupling SSMB storage ring for 


EUV and soft X-ray radiation generation 


Parameter Value Description 
Eo 800 MeV Beam energy 
Co ~100 m Circumference 
Bring 1.33 T Bending field in the ring 
Pring 2m Bending radius in the ring 
oss 4.85 x 1074 Natural energy spread 
€y 2 pm Vertical emittance 
AL 270 nm Modulation laser wavelength 
oj 10 um Transverse electron beam size at the 
radiator 
EUV (13.5 nm) oz (Rad) 2 nm Linear bunch length at the radiator 
ozy (Mod) 1.85 um Bunch lengthening from ey at the 
modulator 
h 541 m7! Energy chirp strength 
AuMod 0.5m Modulator undulator period 
BoMod 0.039 T Modulator peak magnetic flux density 
LuMod 1.5m Modulator length 
PL(ZR = 3H) 141 kw Modulation laser power 
AR= wa 13.5 nm Radiation wavelength 
b20 0.11 Bunching factor 
AuRad 2cm Radiator undulator period 
BorRad 1.15T Radiator peak magnetic flux density 
LuRad 3.2m Radiator length 
Pr 1kW @Ip=1.5A_ | Radiation peak/average power 
Soft X-ray (6.75 nm) | oz (Rad) 1 nm Linear bunch length at the radiator 
ozy (Mod) 1.85 um Bunch lengthening from éy at the 
modulator 
h 1082 m7! Energy chirp strength 
AuMod 0.5m Modulator undulator period 
BoMod 0.039 T Modulator peak magnetic flux density 
LuMod 1.5m Modulator length 
PL(ZR= fu) 564 kW Modulation laser power 
AR= A 6.75 nm Radiation wavelength 
b40 0.085 Bunching factor 
AuRad 1.5 cm Radiator undulator period 
Borad 111T Radiator peak magnetic flux density 
LuRad 24m Radiator length 
PR 1kW @ Ip =2.2A Radiation peak/average power 
Soft X-ray (2.7 nm) | o;(Rad) 0.5 nm Linear bunch length at the radiator 
ozy (Mod) 2.81 um Bunch lengthening from €y at the 
modulator 
h 1423 m7! Energy chirp strength 
AuMod 0.6 m Modulator undulator period 
BoMod 0.028 T Modulator peak magnetic flux density 
LuMod 1.8m Modulator length 


(continued) 
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Table 6.2 (continued) 


Parameter Value Description 

Pi (ZR = tu) 1 MW Modulation laser power 

AR= ih 2.7 nm Radiation wavelength 

b100 0.049 Bunching factor 

AuRad 1 cm Radiator undulator period 

BorRad 0.86 T Radiator peak magnetic flux density 
LuRad 2m Radiator length 

PR 1kW @Ip=5A Radiation peak/average power 


circumference of 100 m should be feasible. This compact high-power EUV radiation 
source is promising to fulfill the urgent need of EUV lithography for high volume 
manufacture, and also serve the future lithography like Blue-X which invokes 6.x nm- 
wavelength light source. Such an SSMB-based high-power soft X-ray photon source 
could be of great value for fundamental science like high-resolution angle-resolved 
photoemission spectroscopy and can also bridge the water window gap. 
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